UNCLASSIFIED 


AD  NUMBER 


AD819970 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  28  JUN 
1963.  Other  requests  shall  be  referred  to 
Army  Air  Mobility  Reseacrch  and 
Development  Lab.,  Fort  Eustis,  VA. 


AUTHORITY 


USAAMRDL  ltr,  4  May  1971 


THIS  PAGE  IS  UNCLASSIFIED 


Statements  and  opinions  contained  herein  are  those  of  the  authors 
and  are  not  to  be  construed  as  reflecting  the  view*  of  the  Department  of 
the  Army  or  the  cosponsors  of  this  Symposium- -The  U.  S.  Army 
Transportation  Research  Command  or  Cornell  Aeronautical  Laboratory,  Inc. 


FOREWORD 


Rapid  advances  made  in  helicopter  and  V/STOL  aircraft  development 
during  the  past  few  years  have  created  many  new  and  unique  design  require¬ 
ments.,  In  the  field  of  dynamic  loads,  in  particular,  an  increasing  number 
of  problems  have  required  fundamental  research  in  mechanics,  aeromechanics, 
and  physics  to  acquire  the  understanding  necessary  for  the  development  of 
satisfactory  method*  of  solution.  Specific  research  has  resulted  in  significant 
advances  in  the  "state  of  the  art"  of  prediction  techniques  for  dynamic  loads 
on  helicopter  and  V/STOL  aircraft.  The  objectives  of  this  Symposium  are  to 
aid  the  dissemination  of  the  latest  research  results  and  to  provide  a  timely 
opportunity  for  open  discussion  among  specialists  in  the  field  of  current  and 
future  "dynamic  load  problems  associated  with  helicopters  and  V/STOL 
aircraft". 

In  keeping  with  these  objectives,  the  first  two  days  of  the  Symposium 
are  devoted  to  the  presentation  and  discussion  of  specific  technical  papers 
on  various  aspects  of  dynamic  load  problems  which  run  the  entire  gamut 
from  basic  research  to  design  methods.  In  addition,  n.  third  day  is  provided 
for  a  review  of  the  current  and  future  dynamic  load  problems  by  panels 
composed  of  outstanding  representatives  from  the  helicopter  and  V/STOL 
aircraft  industries. 

The  proceedings  of  this  Symposium  have  been  bound  in  three  volumes  — • 
one  volume  for  each  day  of  the  general  meetings,  and  one  for  the  panel  meetings. 
We  are  indebted  to  the  authors  for  preparing  their  final  manuscripts  in  a  form 
that  could  be  directly  reproduced.  This  material  was  published  herein  as 
provided  by  the  authors  and  was  neither  checked  nor  edited  by  CAL  or  TRECOM. 
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DYNAMIC  LOADS  IN  A  RIGID  ROTOR 
FROM  FULL-SCALE  WIND-TUNNEL  TESTS 

J.F.  Johnston,  J.E.  Rhodes,  and  R.H.  Cotton 
Lockhood,  California 


DYNAMIC  LOADS  IN  A  RIGID  ROTOR  FROM 
FULL  SCALE  WIND  TUNNEL  TESTS 


J.  F.  Johnston,  Manager,  Flight  Dynamics  Department 
J.  E„  Rhodes,  Group  Engineer,  Helicopter  Dynamic  Loads 
R.  H.  Cotton,  Res  ear ch  Specialist,  Flight  Analyses 

INTRODUCTION 

In  setting  up  the  high  speed  rigid  rotor  research  program  on  the 
XH-51A,  it  was  soon  recognized  that  a  wind-tunnel  test  program  prior 
to  flight  would  expedite  and  increase  the  safety  of  the  flight  test  program. 
With  the  aid  of  BuWeps,  TRECOM,  and  NASA  personnel,  it  was  decided 
that  a  full-scale  test  using  as  much  as  possible  of  the  flight  hardware 
would  be  most  advantageous.  In  addition  to  verifying  and/or  improving 
the  loads,  control,  and  performance  prediction  methods  for  the  rigid 
rotor,  such  a  test  would  also  provide  a  good  checkout  of  the  flight 
components.  Accordingly,  these  tests  were  conducted  in  the  NASA  Ames 
40  x  80  foot  wind-tunnel  in  two  phases  during  February  and  August,'  1962. 
This  paper  discusses  the  dynamic  stability  and  dynamic  loads  aspects  of 
those  te  sts. 

Additional  material  relative  to  the  XH-51A  rigid  rotor  research 
helicopter  is  available  in  References  1,  2,  and  3.  Reference  1  covers 
some  of  the  basic  design  considerations,  particularly  in  the  dynamics 
area.  Reference  2  the  loads  prediction  methods,  and  Reference  3  is  an 
overall  progress  report  on  the  research  program  up  to  May  1963. 


BACKGROUND 


During  the  CL-475  program  a  whirl  tower  test  method  for  the  rotor  system 
had  been  evolved  which  contained  fuselage  pitch  and  roll  representation,  Figure  1. 

In  essence  the  rotor  was  mounted  on  an  inertia  frame  which  represented  the 
fuselage.  This  simulated  fuselage  contained  all  the  rotor  controls  and  a  pilot's 
seat  from  which  the  rig  could  be  controlled  in  a  manner  similar  to  that  on  the 
actual  flying  machine.  Power  to  the  rotor  was  fed  from  an  external  source.  The 
inertia  frame  was  gimbaled  to  a  fixed  frame  rigidly  attached  to  the  deck.  The 
pilot  on  the  whirl  tower  was  thus  able  to  apply  given  lifts,  pitch  and  roll  about  the 
gimbal.  While  this  proved  an  excellent  tool  for  checking  rotor  stability  and  control 
characteristics  in  hover  the  essential  element  of  forward  velocity  was  missing. 

The  rotor  blade  was  fully  instrumented  and  it  was  proposed  to  go  up  in 
speed  on  the  vehicle  in  small  increments  carefully  monitoring  loads  and  stability 
trends  at  each  step. 

The  problem  presented  in  this  approach  was  the  fact  that  one  of  the  more 
critical  areas,  transition,  was  approached  almost  instantaneously  in  the  15  m?h  to 
35  mph  range  (Figure  2  shows  the  typical  cyclic  load  plot  vs.  speed).  To  provide 
coverage  in  this  rSnge  an  Electra  turboprop  aircraft  was  rim  in  front  of  the  whirl 
tower.  This  gave  up  to  a  40  mph  slipstream  which  although  somewhat  turbulent 
provided  the  necessary  assurance.  The  pilot  was  able  to  fly  the  rig  in  this 
turbulence  without  difficulty.  It  was  from  these  early  experiments  that  the  potential 
of  ground  testing  in  the  wind  tunnel  became  apparent. 
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NASA  AMES  40  X  80  WIND 


TUNNEL  TEST  PROGRAM 


CL-475  experience  had  indicated  one  very  important  factor.  In-plane  and 
higher- frequency  loading  levels  and  stability  characteristics  depended  on  the  rotor- 
body  interrelationship.  The  more  dynamic  analysis  would  progress  the  greater 
the  realization  that  the  rotor  was  not  an  entity  on  its  own.  The  dynamic  system 
was  interdependent  (Reference  1).  With  the  help  and  guidance  of  NASA,  TRECOM 
an!  BUWEPS  personnel  an  ambitious  program  for  testing  the  XH-51A  rotor -body 
system  was  evolved.  The  ground  rules  were  as  follows^ 

(a)  The  rotor -body  dynamic  relationship  should  be  simulated  to  the 
degree  to  make  load  level  measurements  meaningful. 

(b)  The  control  system  including  the  prime  stability  items  should  be 
represented. 

(c)  The  complete  system  would  be  flown  with  only  translation  movement 
restricted.  The  control  system  was  to  be  similar  to  that  of  the  vehicle. 

(d)  Provisions  should  be  made  for  varying  prime  control  characteristics. 

(e)  Complete  load  measuring  instrumentation  would  be  provided, 

(f)  The  rig  would  be  designed  to  allow  its  complete  functional  check  out 
before  installation  at  AMES. 

The  succese  of  the  program  was  dependent  on  three  prime  items. 

(1)  Aerodynamic  simulation 

(2)  Dynamic  simulation 

(3)  Ability  of  the  control  personnel  to  be  able  to  fly  the  rig  from  a  remote 
station. 

The  only  question  with  respect  to  aerodynamic  simulation  was  the  effect  of  the  tunnel 
walls  and  posoible  free  stream  turbulence.  This  problem  can  best  be  understood  by 
comparison  of  the  item  to  be  tested  to  the  tunnel  size.  Figure  3  shows  the  rig 
mounted  in  the  tunnel.  In  retrospect,  aerodynamic  simulation  proved  to  be  excellent 
except  at  very  low  speeds  and  hover.  In  hover  the  rotor  tunnel  wall  effects  created 
considerable  turbulence.  This  problem  was  overcome  by  tilting  the  ivU»*  In  » 
10-degree  nose  down  attitude  to  provide  down  stream  flow. 


Dynamic  simulation,  at  first  very  ambitious,  had  to  be  modified  for 
practical  reasons  and  for  some  of  special  stability  problems  associated  with  the 
teat  rig.  These  will  be  discussed  later. 


DESCRIPTION  OF  TEST  RIG 


Figure  4  shows  the  final  test  rig.  This  rig  could  be  mounted  either  on  a 
support  structure  for  whirl  tower  testing  or  on  the  AMES  tripod  electric  driven 
shaft  system.  The  pitch,  roll  and  collective  control  systems  were  operated  by 
electric  actuators;  the  control  of  these  actuators  was  outside  the  tunnel  in  a 
control  station.  The  mount  system  freedom  was  of  prime  importance  for  dynamic 
simulation.  Figure  5  shows  the  initial  objectives  and  those  finally  used. 

A  fuselage  having  the  external  shape  of  the  XH - 51 A  was  mounted  on  the 
inertia  frame,  the  object  being  to  obtain  as  much  information  on  aero  charac¬ 
teristics  as  possible  from  the  tunnel  tests. 


INSTRUMENTATION 


The  blade  was  instrumented  to  obtain  *he  blade  flap  moments  at  seven 
sp&nwise  stations  and  tl;e  chord  or  in-plane  moments  at  five  stations.  At  any  one 
time  the  complete  apanwise  bendings  on  one  blade  were  obtained  along  with  the 
root  bending  moments  on  the  other  two  blades.  All  blades  were  strain  gauged  so 
thay  could  be  hooked  up  at  any  time.  Figure  6  shows  the  locations  of  the  loads 
measurement*  on  the  hub  and  blades. 


TEST  PROGRAM 


nwrtMuatM  C  CV*C  ~  thw 
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DYNAMIC  STABILITY  CONSIDERATIONS 


Early  consideration  of  whirl-tower  and  wind-tunnel  testing  indicated  two 
distinct  areas  in  which  the  stability  of  the  gyro -controlled  rigid  rotor  would  be 
compromised  by  the  conditions  of  test.  These  problems  were  confirmed  during 
the  testing  and  a  third  was  found.  Fortunately,  solutions  were  available  for 
each  problem  that  would  not  compromise  obtaining  valid  rotor  dynamic  loads  data. 

Whirl-Tower  Mode  -  The  first  problem  area  was  dubbed  the  "whirl -tower 
mode".  It  arose  from  the  use  of  a  remote  power  source  to  drive  the  rotor.  A 
universal  joint  coincident  with  the  body  gimbal  center  was  provided  in  the  rotor 
drive  shaft  to  permit  the  body  to  pitch  and  roll. 

If  the  body  and  rotor  are  moved  in  roll  (for  example)  the  resulting  moment 
system  due  to  the  drive  torque  is  shown  in  Figure  8.  ! 

The  displaced  rotor  shaft  absorbs  the  lower  shaft  drive  torque  partly  as 
torque  in  the  displaced  plane  of  rotation  and  partly  as  a  moment  bending  the  upper 
shaft  forward  -  a  pitching  moment  if  the  rotor  and  shaft  are  displaced  in  roll, 
or  a  rolling  moment  if  the  rotor  is  pitched.  Since  the  rotor  is  flexible,  the  upper 
Shaft  Will  tend  to  pitch  in  the  direction  of  the  unbalanced  moment.  This  will 
slowly  precess  the  gyro  in  pitch,  due  to  friction  and  other  moments  which  are 
always  at  work  tending  to  erect  the  control  gyro  normal  to  the  rotor  shaft.  Thus 
the  gyro  and  rotor  will  have  an  angular  velocity  in  phase  with  ^he  applied  moment. 
If  the  body  is  free  in  both  pitch  and  rol*  about  the  gimbal  s,  the  system  is 
spirally  unstable  in  an  advancing  mode,  the  rate  of  spiral  divergence  being 
proportional  to  the  rotor  torque  and  the  characteristic  erection  time  of  the  gyro. 

Inasmuch  as  this  "whirl-tower  mode"  did  not  show  up  in  the  whirl-tower 
tests  of  the  rotor  system  for  the  lower -powered  CL -47  5,  it  was  hoped  that  the 
higher -powered  XH-51A  system  would  also  be  relatively  free  of  it.  However, 
the  spiral  instability  occurred  during  early  high-iitt,  mgh-power  teste,  tie  cove  ry 
was  made  by  caging  the  body.  Further  difficulty  was  avoided  by  stiffening  the 
body  mounting  in  roll,  Benv4ng.it> free  to  pitch.  The  roll  control  inputs  were  there¬ 
after  based  on  roll  momenta  rather  than  roll  attitude  idicafcion. 
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It  should  be  reiterated  that  the  ’’whirl-tower  mode"  is  associated  with  a 
remote  drive,  and  cannot  occur  in  free  flight. 

Wind-Tunnel  Pitch-Uj?  -  The  second  area  in  which  the  difference  of  test 
conditions  from  free  flight  introduced  stability  problems  was  the  pitch  stability 
in  the  wind  tunnel  at  forward  speed. 

In  the  wind  tunnel,  pitch  attitude  is  synonymous  with  angle  of  attack, 

ninr-  ii  i  n  Mi  i 

whereas  in  free  flight  the  freedom  in  plunge  makes  it  necessary  to  maintain  a 
pitch  rate  in  order  to  maintain  an  angle  of  attack  other  than  that  for  steady 
flight. 


Pitch  stability  in  the  Lockheed  rigid  rotor  system  is  obtained  from  the  fact 
that  the  control  gyro  tends  to  remain  fixed  in  space,  so  that  a  change  in  rotor 
attitude  produces  a  corrective  feathering  input.  The  gyro  will  however,  process 
slowly  to  erect  nprnal  to  the  shaft,  slowly  removing  the  corrective  input.  In 
free  flight  this  erection  rate  is  small  compared  to  the  rate  of  pitch  required  to 
maintain  an  out-of-trim  angle  of  attack,  with  the  result  that  the  vehicle  is 
stable  with  respect  to  angle  of  attack.  In  the  wind  tunnel,  h-'wever,  this  erection 
rate  produces  a  slow  divergence.  It  was  assun^ed  that  the  operator  controlling 
rotor  attitude  could  handle  -this  slow  divergence  satisfactorily.  The  control 
rates  available  were  slow  enough,  however,  to  cause  the  operator  considerable 
trouble.  During  a  75-knot  run  it  became  necessary  to  cage  the  system  due  to 

a  pitch -up.  The  remainder  of  the  tests  at  and  above  75  knots  were  run  with  the 

.....  ,  • 

fcyStdtil  caged,  although  it  later  became  necessary  to  reduce  the  caging  stiffness 
in  pitch  in  order  to  provide  the  solution  to  the  third  potential  dynamic  instability 
associated  with  the  wind  tunnel  test  conditions. 

Regressive  Whirl  Mode  -  The  characteristic  wobble  of  a  gyroscope  will 
be  in  an  advancing  or  regressing  mode  depending  on  whether  the  forces  acting 
on  it  are  destabilising  or  stabilising,  respectively.  A  toy  gyroscope  spinning 
atop  a  table  (force  of  gravity  destabilising)  will  wobble  in  the  same  direction  as 
it*  "srhtrr*;  the  uwl«w  a  wiring  {force  ox 

gravity  stabilising)  will  wobble  in  a  direction  opposite  its  rotation  after  a  disturbance. 

In  changing  by  stages  from  a  system  free  in  both  roll  and  pitch  to  a  system 
restrained  about  both  aso»»,  the  characteristic  dynamic  motion  of  the  rotor  had 


6 


changed  from  an  advancing  mode  characteristic  of  free  flight  to  a  regressive  mode. 
The  effect  of  this  change  began  to  appear  in  testing  at  100  knots  at  a  10-degree 
nose  down  attitude,  when  a  beat  frequency  began  to  appear  in  the  traces  of  the 
blade  chord  bending.  .A  typical  trace  is  shown  in  Figure  9.  Inasmuch  as  the  beat 
frequency  was  very  close  to  the  regressive  chord  bending  frequency  as  seen  in 
the  non-rotating  system,  it  was  concluded  that  the  trace  indicated  an  incipient 
dynamic  instability.  The  boundaries  for  these  occurrences,  in  terms  of  rotor 
attitude  versus  speed,  are  .shown  in  the  Io-yer  part  of  Figure  9. 

Fortunately,  dynamic  analyses  of  the  wind  tunnel  test  condition#  were  in 
process  at  the  same  time.  These  analyses  showed  the  possibility  of  such  an 
instability.  More  importantly,  they  showed  how  to  cure  the  instability.  The 
results  are  given  in  Figure  10,  showing  the  system  stability  as  a  function  of  the 
body  pitching  frequency  against  the  wind  tunnel  balance  tripod.  The  system  is 
stable  for  all  body  pitch  frequencies  below  1,  6  cps,  and  unstable  for  all  higher 
frequencies.  The  cut-off  frequency  of  1.  6  cps  is  the  regressive  chord  bending 
frequency  as  seen  by  the  body.  At  support  frequencies  below  1.  6  cps,  the  body 
would  pitch  out  of  phase  with  the  chord  bending,  and  the  regressive  whirl  mode 
pattern  would  be  broken  up.  The  support  stiffness  was  then  modified  by  placing 
rubber  blocks  on  the  balance  tripod  platform  against  which  the  body  was  pulled 
by  the  caging  cylinders,  and  adjusting  the  spacing  of  the  rubber  blocks  until  the 
body  pitch  frequency  was  reduced  to  1.  5  cps.  This  minimum  margin  of  0. 1  cps 
was  selected  in  order  to  maintain  the  maximum  possible  pitch  stiffness 
{approximately  150,000  ft.  lb.  per  radian),  and  because  the  physical  relation- 

i  ■ 

chips  indicated  by  the  dynamic  analysis  made  good  sense.  From  this  point 
testing  proceeded  up  to  the  maximum  speed  with  no  further  dynamic  problems. 

Adequacy  of  Simulation  -  From  these  experiences  it  can  be  concluded  that 
the  low-frequency  dynamics  of  the  rigid-rotor  helicopter  cannot  be  simulated 
with  a  restrained  model  in  the  wind  tunnel.  Th®  simulation  improve  a  for  frequencie 
Above  the  body  restraint  frequencies,  inasmuch  as  the  phase  of  the  body  response 
becomes  similar  to  free  flight.  The  relative  amplitudes  of  body-rotor  interactions 
are  normally  satisfactorily  simulated  at  frequencies  above  twice  the  body 
restraint  frequency. 

With  respect  to  the  rotor  dynamic  loads  for  steady  flight  conditions  - 
level  flight  at  any  simulated  gross  weight,  drag,  or  longitudinal  or  lateral  e,  g. 
displacement  -  thfe  simulation  is  very  good.  The  IF  vertical  bending  accurately 
reflects  She  blade  loads  due  to  steady  body  momenta.  The  IF  bending 
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requires  no  body  reaction  and  is  thus  independent  of  the  conditions  of  restraint. 

The  2P  and  4P  vertical  and  in-plane  blade  loads  in  a  3 -  blade  rotor  are  reflected  as 
3P  to  the  body,  which  is  well  above  the  body  restraint  frequency.  Normally  these 
blade  loads  will  be  well  simulated  if  the  body-rotor  springing  is  simulated,  although 
the  body  responses  are  more  open  to  question. 

The  3F  blade  loads  are  collective  in  a  3-blade  rotor.  The  vertical  bending  in 
most  important,  inasmuch  as  the  aerodynamic  input  is  appreciable  and  the  blade 
second  bending  frequency  is  near  3P.  Analyses  indicated  that  as  long  as  the  body 
appeared  massive  to  the  rotor,  the  degree  of  massiveness  had  little  effect  on  the 
blade  3P  responses  in  the  rvnge  investigated,  which  was  from  one -quarter  to  full 
effective  body  mass.  Vortical  springiness  between  the  rotor  shaft  and  the  body 
may,  however,  have  a  strong  effect  on  the  blade  loads.  Proper  simulation  is 
important  if  this  connection  is  relatively  flexible. 

The  3P  in-plane  responses  were  not  properly  simulated  because  of  the 
high  inertia  of  the  drive  mechanism,  but  these  have  proved  negligible  both  in 
flight  and  in  the  wind  tunnel. 


ROTOR  DYNAMIC  LOADS 


The  methods  used  in  loads  prediction  work  was  based  on  Reference  2.  The 
level  flight  loads  are  dependent  on: 

1.  Body  moment  due  to  center  of  gravity  offset  and  aerodynamic 
pitching  moment. 

2.  '  Vehicle  attitude  and  speed  which  dictate  collective  and  cyclic 

blade  angles. 

In  the  A mas  tunnel  the  attitude  of  the  vehicle  was  set  by  tilting  the  tripod 
tiwuui.  Tile  aiitiuua  was  based  on  an  assumed  arpg  tor  tne  vehicle.  A  check  on 
this  drag  from  tunnel  balance  measurement  was  difficult  due  to  the  large  tare 
drag  of  the  tripod.  In  order  to  provide  a  complete  coverage  the  vehicle  attitude 
was  varied  at  each  speed  point  up  to  120  knots. 

In  comparing  predicted,  in-flight,  and  Ames  data  it  should  be  realised  that 
bhikitiing  ss  vehicle  condition  which  is  identical  for  all  three  in  term*  of  body 

lift,  speed,  drag,  -etc*  la  not  always  possible.  The  load  levels  therefore. 


can  only  bo  thought  of  as  general  trends  as  a  function  of  speed.  At  Ames,  for 
example,  the  operators  had  to  iqrim  out  both  pitch  and  roll  moments  using  an 
indicator  linked  to  the  load  cells.  At  the  instant  of  taking  the  record,  there  was 
often  quite  &  high  level  of  body  moment  left  untrimmed. 

In  spite  of  these  variables  the  plot  of  cyclic  in-plane  moments  at  the  hub 
versus  speed.  Figure  11,  shows  good  agreement  with  predicted  values.  These 
valuers  are  predominantly  one  per  revolution  and  first  mode. 

Good  agreement  with  the  one  per  revolution  first  mode  flapping  and 
steady  flapping  was  obtained  when  full  consideration  was  taken  of  the  body 
moment  and  lift  for  the  particular  condition.  Figure  12  Jhcws  the  IP  flapping 
moment  distribution  comparison  between  measured  and  predicted  values. 


Loads  prediction  for  the  2P,  3P,  and  higher  harmonics  was  less  sophisti¬ 
cated  at  the  time  the  wind  tunnel  tests  were  made,  inasmuch  as  the  levels 


extrapolated  from’cfe-475  experience  indicated  Ikit  feiey  were  of  less  importance 


to  rotor  fatigue  than  the  effects  of  c.g.  offset  and  fuselage  aerodynamic  moments. 
The  early  analyses  did  not  include  considerations  of  vortex  shedding  at  low 
speed  or  of  the  energy  admittance  of  the  harmonic  airloads  into  the  higher  modes 
of  blade  bending. 


Figure  13  shows  comparisons  of  the  hub  2P  and  3P  flapping  moments  vs. 
speed.  The  2P  moments  measured  in  the  tunnel  were  generally  higher  than 
predicted.  These  moments  were  effectively  removed  in  the  flight  article  by 
providing  a  shaft  mounting  frequency  below  3P.  This  resulted  in  a  17%  reduction 
of  the  pitch  stiffness  of  the  body  relative  to  the  rotor  plane.  This  elimination  of 
the  2P  blade  loading  is  an  example  of  the  importance  of  the  body-rotor  dynamic 
interaction^,  Referencel. 


The  3P  hub  flapping  moments  were  found  to  be  of  consequence  throughout 
the  flight  regime.  The  xlight  results  are  in  general  agreement  with  the  wind- 
tunnel  up  to  100  Knots,  and  indicate  a  higher  rate  ox  increase  above  that  speed. 

Body -rotor  dynamic  interaction  is  prime  importance  in  the  2P  and  3F 
dynamic  response.  The  XH-51A  Rotor  has  a  second  mode  fi&ppihg  frequency 
very  near  3?  as  shown  in  Figure  14,  The  exact  phasing  and  combined  frequency 
is  dependent  on  the  mass  springing  below  the  rotor.  Figure  IS  which  is  extracted 
from  Reference  1  shows  this  relationship. 
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A  reasonable  tolerance  on  the  higher  frequency  loadings  is  acceptable,  ai  j 
they  contribute  less  than  50%  of  the  total  cyclic  moments  causing  fatigue  damage. 
Figure  16  shows  the  total  cyclic  flapping  moment  at  the  hub  versus  speed  measured 
at  Ames,  and  the  contribution  of  the  3F  component. 

FULL  SCftLE  VS.  MODEL  TESTING 


The  prime  advantage  of  full-scale  wind-tunnel  rotor  testing  as  compared 
with  testing  at  model  scale  is  that  actual  flight  hardware  is  tested,  minimizing 
the  difficulties  in  interpretation  from  wind  tunnel  to  flight  results. 

Aside  from  the  obvious  advantages  of  testing  at  true  Reynolds  and  Mach 
numbers,  the  full  scale  hardware  has  the  further  advantage  of  including  the  true 
effects  of  the  non-linearities  due  to  friction  and  dynamics  of  the  various  mechanisms. 
Mechanical  malfunction  is  often  the  trigger  for  catastrophic  failure.  A  good 
overall  dynamic  design  can  be  ruined  by  poor  mechanical  detail.  Full  scale  testing 
is  the  only  way  to  cover  this  important  area.  The  Ames  test  uncovered  two 
important  problems  in  this  category. 

The  first  was  the  failure  of  a  blade  due  to  the  separation  of  a  poorly  bonded 
leading  edge  balance  as  described  in  Reference  3.  This  failure  led  to  complete 
redesign  of  the  blade  to  provide  redundancy  and  inspectability  in  the  blade -$ub 
system.  The  discovery  of  this  weakness  alone  justified  the  full  scale  wind  tunnel 
program. 

A  second  less  catastrophic  discovery  involved  the  feathering  bearings.  In 
surveying  the  records,  the  feathering  angle  trace  ehow^d^.flat  top  to  the  sine 
curve.  Examination,  of  the  bearings  showed  that  they  tended  to  seize  duerto  locking 
of  the  needles  in  the  roller  bearings.  The  design  was  immediately  improved  to 
correct  the  situation. 

The  bearing  sticking  would  have  shown  up  in  the  flight  vehicle  as  poor  and 
erratic  dynamic  characteristics.  The  true  3,  might  have  been  difficult  to 

trace  in  flight  tests  because  of  the  reduced  instrumentation  possible  in  the  flight 
vehicle. 
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The  prime  disadvantage  of  full-scale  rotor  testing  Xi«s  in  the  large  rotational 
energy  in  the  rotor,  making  the  tests  always  somewhat  hazardous.  For  this  reason, 
no  testing  in  known  hazardous  areas,  such  as  fully  developed  retreating  blade  stall, 
was  planned.  The  exploration  of  hazardous  areas  is  better  done  at  model  scale  where 
the  danger  to  test  personnel  can  be  orders  of  magnitude  less,  and  the  cost  of  the 
many  replacement  parts  likely  to  be  required  is  correspondingly  less. 


CONCLUDING  REMARKS 


In  reviewing  this  experience  of  testing  the  XH-51A  rigid  rotor  up  to  160  knots 
in  the  NAS/.  Axhes  40  x  80  foot  wind  tunnel,  it  was  co:. -',luded  that  the  prime  objectives 
had  been  achieved.  The  hardware  had  been  wrung  out  and  important  mechanical 
deficiencies  had  been  found  and  corrected.  Proof  and  further  understanding  of 
the  fundamental  dynamic  behavior  of  the  system,  within  the  limits  imposed  by  the 
necessary  wind  tunnel  restraints,  had  been  obtained.  The  blade  dynamic  loads 
associated  with  steady  flight  hand  been  determined. 

Should  all  new  rotor  systems  be  tested  in  this  manner?  fphere'is  probably 
no  general  answer  to  this  question.  Full-scale  wind  tunnel  testing  is  not  a  substitute 
for  the  necessary  development  flight  testing;  it  is  rather  a  means  of  increasing  the 
safety  of  the  flight  work  and  of  getting  qm  early  lead  on  some  of  the  problems  likely 
to  be  encountered  in  flight.  It  is  more  valuable  and  more  costly  than  model 
testing  and  thus  will  command  more  attention  and  higher -grade  personnel.  Inasmuch 
as  flight  hardware  is  involved,  details  that  might  be  passed  off  as  "bugs  in  the 
model"  cannot  be  ignored  in  full-scale  testing. 

In  short,  each  situation  must  be  judged  on  its  merits.  In  this  particular 
case,  full-scale  testing  had  paid  off  hands  oGnedy* 
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SUMMARY 


Presented  results  of  an  -aeroeiastic  investigation  to  explore  effects 
of  flight,  rotor,  aiid  bld'de  parable  ter  s  op  blade  stresses  and  performance  of 
articulated  dnd  rigid  rotor  systems  ih  high-speed  flight.  Variations  in  blade 
planform  dnd  twist,  fdtor  lift  and  brppulsive  force,  blade  stiffness  and  mass 
distribution,  are  investigated  for  helicopters  and  compound  helicopters  at 
several  groSs  weights  for  speeds  from  150  to  300  knots.  The  analysis  takes 
into  account  the  fully  coupled  flsttwise -edgewise- torsional  response  of  the 
rotor  blades.  Work  was  performed  under  U,  S.  Army  TRECOM  Contract 
DA  44- 177 -TC- $31. 


INTRODUCTION 

Design  of  helicopters  and  compound  air  vehicles  for  higher  airspeeds 
requires  greater  knowledge  of  rotary  wing  dynamics  and  performance  than  is 
presently  available  from  published  analyses  and  investigations.  For  fuller 
understanding  of  advanced  flight  regimes,  a  research  program  is  required, 
which  will  explore  effects  of  forward  speed,  blade  planform,  blade  twist,  rotor 
lift  and  propulsive  force,  blade  stiffness  and  mass  distribution  on  resulting 
blade  stresses  and  performance  of  articulated  rigid  rotor  systems. 

To  provide  this  information,  Sikorsky  Aircraft  last  year  undertook  a 
study  for  the  U.  S. ,  Army  Transportation  Research  Command.  This  paper 
gives  prelimihary  results  of  that  study.  The  work  was  achieved  by  a  joint 
effort  of  the  company’s  Aerodynamic,  Blade  Design,  and  Dynamics  Sections. 
Three  well  developed  analytical  programs  were  used  in  the  program.  These 
were:  (X)  an  Advanced  Performance  Analysis,  which  was  used  to  determine 
the  basic  rotor  parameters  to  be  studied;  (2)  a  Horvay-type  Blade  Design 
Analysis,  which  has  been  responsible  for  Sikorsky’s  present-day  successful 
blade  designs;  and  (3)  a  Coupled  Blade  Flatwise-Edgewise-Torsionai  Aero- 
elastic  Analysis,  which  yields  the  full  spectrum  of  blade  dynamic  information. 
The  aeroelasrir  analysis  ie  general  and  weii  smtea  tor  studies  of  a  wide  range 
of  helicopters  and  VTOL-type  aircraft  with  rotbr  blades  and  propellers  of  all 
types.  As  will  be  shown  in  die  paper,  good  correlation  has  been  achieved  with 
the  method  for  both  rigid  and  articulated  rotor  systems. 
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1.  DESCRIPTION  OP  AIRCRAFT  CONSIDERED 

The  aeroelaatic  research  program  is  being  carried  out  under  U,  S.  Army 
TRECOM  Contract  DA  44-177-TC-831.  Included  in  the  overall  TREGOM  study 
are  eight  basic  aircraft,  four  in  the  12,000-lb,  class  and  four  in  the  33,000-lb. 
class.  There  are  two  150-knot  helicopters,  two  200-knot  helicopters,  two 
200-knot  wingless  compounds,  and  two  300-knot  winged  compounds.  Four 
planform  variations  are  considered  in  each  case,  and  blades  are  taken  as  both 
rigid  and  articulated.  Planform  variations  are  given  in  Figure  (1).  A  total  of, 
470  flight  conditions  will  be  explored  in  the  contract. 

Results  presented  in  this  paper  are  taken  from  initial  investigations  of 
the  program.  Considered  are  the  150-knot  helicopters  at  12, 000  and  33, 000 
lb.  gross  weights.  Also  shown  are  some  results  for  the  12, 000-lb.  winged 
compound.  Substantiation  for  the  method  of  analysis  is  based  upon  correlation 
studies  with  flight  test  data  from  the  Lockheed  CL- 475  rigid  rotor  helicopter, 
and  Sikorsky's  S-58  and  S-61  articulated  helicopters.  Basic  data  for  the  five 
aircraft  considered  here  are  presented  in  the  following  table: 

TABLE  I 


BASIC  DATA  FOR  AIRCRAFT  INVESTIGATED 


Lockheed 
CL- 475 

.Sikorsky 

S-58 

Sikorsky 

S-61 

12000- lb. 
Helicopter 

33000-lb. 

Helicopter 

12000rlb, 

Compound 

Blade  Radius 
(ft.) 

16 

28 

31 

31 

36 

28 

Offset  (inches) 

IX.* 

12 

12. 625 

12. 625 

24 

12.63 

Number  of 
Blades 

3 

4 

5 

5 

6 

i  .5 

Chord  (inches) 

12.0 

16.4 

18.  25 

18.  25 

23.65 

18.25 

Design  Taper . 

1:1 

1:1 

1:1 

1:1 

1:1 

i  1:1 

Design  Twist 
(degrees) 

0 

-8 

-8 

(-8  artic.  (-6  artic.  (-4  artic. 
(  0  rigid  (  0  rigid  (-4  rigid 

Airfoil  Section 

0012 

0012 

.0012 

,0012 

0012  i 

0012 

Solidity 

.  0596 

.062 

. 07804 

.  07S04 

.1046 

.056 

Tip  Speed  (&>s) 

504 

620 

660 

660 

696 

493cruisi 

Rotor  Speed 
(rpm) 

300 

216 

203 

203 

185 

168 

Gross  Weight 
(lbs.) 

1970 

11300 

17000 

12000 

33000 

12Q00 

Percent  G.  W. 
Carried  by 

Rrrtw** 

100 

ion 
v  - 

*8  nn 

r  aw 

1  AA 

« 

i  *  AA 

|  AW 

i  .. 

’Pitch  Bearing  Center  Line 


Tabulated  in  the  Appendix  is  k  set  of  representative  data  for  blade 
stiffness,  and  mass  distribution  taken  from  this  program.  Rigid  and  articu¬ 
lated  blade  data  are  shown  for  the  150-knot  33, 000-lb.  helicopter.  Also  given 
are  fuselage  lift  and  drag  characteristics  for  this  aircraft. 


2.  BLADE  DESIGN  CRITERIA 


Individual  blade  designs  were  determined  as  follows.  Initially,  an 
Advanced  Performance  Analysis  was  used  to  obtain  basic  rotor  parameters. 
These  include  tip  speed,  rotor  diameter,  blade  chord,  number  of  blades,  and 
design  twist,  Based  on  these  parameters,  the  Blade  Design  Section  then 
developed  structural  character! stiqs  for  the  articulated  1:1  planform  blade  for 
each  aircraft  studied.  Blades  were  designed  by  standard  Sikorsky  procedures 
and  met  basic  load  requirements.  Characteristics  of  separate  planform 
blades  were  then  extrapolated  from  1:1  nlanfo-m  designs  by  the  following 
criteria. 


1.  At  corresponding  blade  radial  stations  cross-sectional  area 
was  held  equivalent.' .  By  this  means  blade  weights  and  centrifu¬ 
gal  stresses  were  kept  equal.  The  only  deviation  permitted 
was  when  extrapolated  wall  thickness  became  small.  Here,  to 
avoid  local  spar  buckling  a  minimum  value  of  0. 15  inches,  was 
set. 

2.  Section  properties  were  extrapolated  based  on  the  ratio  of  spar 
wall  thickness  to  chord.  Knowing  chord  and  section  area,  this 
value  could  be  determined,  then  used  to  calculate  flatwise, 
chordwise,  and  torsional  stiffness  assuming  a  typical  D-spar 
blade. 

For  rigid  rotor  systems  the  blade  root  had  to  be  specifically  designed  to 
meet  Lockheed's  criteria  as  set  forth  in  Reference  (1): 

1.  For  planforms  of  aspect  ratio  greater  than  18  a  soft  inplane 
system  was  used.  This  required  the  first  inplane  bending 
frequency  to  be  between  ,  6X1  and  .  in.  The  inpiane  frequency 
was  achieved  by  reducing  root  stiffnees  by  removing  material 
and  attributing  15%  of  the  frequency  decrease  to  bearing 
flexibility. 

2.  For  planforms  of  lees  than  18  a  stiff  inplane  system  was  used. 
Here,  blades  were  tuned  so  that  the  first  inplane  bending  mode 
was  near  1.  4DU 


3,  Flatwise  flexibility  of  rigid  blades  was  designed  sq  that  the 
second  flatwise  sending  mode  was  sufficiently  removed  from 
n/rev. 


4.  Prelag  and  precone  angles  were  selected  based  on  zero  steady 
bending  moment  at  the  blade  root  in  cruise  flight  condition. 


3.  description  OF  AEROELASTIC  ANALYSIS 

The  method  of  aeroelastic  analysis  is  based  upon  superposition  of 
separate  harmonics  of  blade  forced  response,  which  results  from  response  of 
the  blade  to  individual  harmonics  of  airloads.  A  detailed  description  of  the 
method  has  been  given  in  References^)  and  (3).  A  brief  description  follows. 

a.  Calculation  of  Aerodynamic  Loads 


For  the  helicopter  to  be  analyzed,  the  gross  weight,  drag,  speed, 
and  rotor  rpm  must  be  given.  Also  for  the  blades,  steady-state  two- 
dimensional  airfoil  characteristics,  structural  stiffnesses,  mass  distri¬ 
bution,  twist,  and  root  retention  must  be  specified.  The  rotor  disk  is 
considered  to  be  moving  at  the  proper  forward  tilt  to  provide  enough 
propulsive  force  to  overcome  the  net  drag  of  the  aircraft.  It  must  also 
support  the  aircraft,  and  there  must  be  sufficient  cyclic  pitch  to  keep 
the  rotor  in  equilibrium.  Certain  simplifying  assumptions  are  made  to 
initiate  the  calculation;  such  as  the  approximate  coning  angle,  an 
estimate  of  the  rotor  drag,  and  art  estimate  of  the  radial  position  of  the 
resultant  thrust  vector.  These  approximations  do  not  affect  the  final 
accuracy,  for  if  they  are  too  far  in  error,  this  can  be  remedied  by  a 
second  or  third  iteration. 

For  a  high-speed  condition,  constant  inflow  is  taken.  The  blade  is 
subdivided  into  twenty -four  elements.  For  each  of  36  ten-degree 
azimuth  intervals,  the  blade  is  considered  set  at  two  blade  angles. 

These  angles  bracket  the  expected  blade  angles.  Blade-element  aero¬ 
dynamic  lifts  are  then  computed,  from  which  the  moment  of  the  thrust 
about  the  flapping  hinge  is  calculated  as  a  function  of  blade  angle  and 
azimuth  position.  The  cyclic  pitch  necessary  to  maintain  the  rotor 
system  in  equilibrium  is  then  calculated  by  an  iteration  to  enforce  the 
condition  that  the  first  harmonic  thrust  moment  about  the  flapping  hinge 
is  zero. 

For  calculation  of  aerodynamic  loads,  two-dimensional  0012  airfoil 
data  is  used.  Compressibility  effects  are  taken  into  account  by  using 
separate  Cl,  Cq,  and  Gm  versus  curves  for  angles  of  attack  from 
0  to  30  degrees  for  Mach  numbers  up  to  0. 95  in  5%  increments.  To 
define  stall  regions  above  angies-of-attack  of  30  degrees,  single  Cl»  ' 
Gy,  and  Cm  versus  curves  arc  taken  from  data  in  Reference  (8). 

Final  determination  of  cyclic  pitch  yields  angle-of -attack  distribu¬ 
tion,  rotor  drag,  power  required,  location  of  resultant  thrust  vector, 
thrust  moments,  and  provides  pitching  moments,  and  resolved  thrusts 
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and  drags  on  twenty-four  biade  elemen  j  for  ten- degree  azimuth 
intervals.  A  harmonic  analysis  is  performed  on  this  loading  and  the 
steady  plus  the  first  seven  harmonics  of  blade  element  loading  (pitching 
moments,  thrusts,  and  drags)  are  obtained  in  complex  form. 

b.  Blade  Dynamic  Response 

Individual  harmonics  of  airloads  are  next  introduced  into  the 
coupled  biade  dynamic  analysis.  This  method  is  based  upon  an  extension 
of  Myklestad's  analysis  for  rotating  beams.  There  is  provision  for  up 
to  twenty-four  flatwise,  edgewise,  and  torsional  degrees -of- freedom 
with  coupling  due  to  twist.  Equations  are  in  complex  form  to  allow  for 
aerodynamic  damping  and  phasing  of  aerodynamic  loads.  Boundary 
conditions  at  the  tip  of  the  biade  require  that  shears,  moments,  and 
torque  are  zero.  At  the  root  of  the  blade  boundary  conditions  are 
applied  consistent  with  blade  root  restraint,  whether  articulated,  teeter¬ 
ing,  or  rigid.  There  is  provision  for  a  lag  damper.  Also,  with  torsion 
added,  provision  has  her  made  for  control  system  flexibility  and 
damping. 

Analysis  for  flexible  blade  dynamics  with  torsional  coupling  is 
done  in  three  steps.  In  the  first  step,  the  blade  is  treated  as  infinitely 
rigid,  and  rigid  body  coning,  lagging,  and  twisting  calculated.  In  the 
second  step  using  an  iterative  relaxation -type  method,  steady  coupled 
flatwise-edgewise-torsional  bending  is  computed.  The  iteration  is 
initiated  by  using  the  blade  slopes  and  deflections  determined  in  the  first 
step.  In  the  third  step,  the  first  through  the  seventh  harmonics  of  blade 
coupled  responses  are  determined.  Steady  response  values  from  the 
second  step  are  used  so  that  torsional  coupling  of  vibratory  forces  times 
steady  deflections  and  steady  forces  times  vibratory  deflections  can  be 
included. 

Total  forced  response  of  the  blade  is  next  determined  by  super¬ 
position  of  separate  harmonics  of  blade  dynamic  response.  This  yields 
the  azimuthwise  distribution  of  moments,  torques,  deflections,  twist, 
and  stresses  at  each  of  twenty™  four  biade  stations  for  ten-degree 
azimuth  intervals. 


4.  CORRELATION  WITH  FLIGHT  TEST  DATA 

Using  constant  inflow,  the  method  of  analysis  described  shows  good 
agreement  with  measurements  in  predicting  performance.  This  is  illustrated 
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favorably  for  the  S- 61  at  two  gross  weight  conditions.  The  computed  values 
include  power  required  for  the  tail  rotor  and  accessory  drives.  Coefficients 
of  drag  for  the  blades  were  incremented  by  a  Ac  0  «  .002  to  account  for 

additional  roughness  of  the  actual  blade  over  a  highly  polished  wind-tunnel 


specimen.  Note  that  good  agreement  is  achieved  even  at  low  airspeeds  where 
effects  of  variable  inflow  have  been  shown  to  be  large.  Comparit  >ns  also 
indicate  that  the  method  described  provides  performance  results  consistent 
with  standard  performance  calculations  in  use. 

Figures  (3)  and  (4)  compare  calculated  blade  vibratory  bending  moments 
and  stresses  with  flight  test  values  for  an  articulated  rotor  helicopter.  Given 
in  Figure  (3)  is  a  comparison  of  the  radial  distribution  of  one-half  peak-to- 
peak  bending  moments  against  flight  test  data  for  the  S-58  (H-34)  helicopter. 
Data  is  shown  for  110  knots.  The  test  helicopter  was  instrumented  by  Sikorsky 
Aircraft  under  U,  S.  Army  TRECOM  contract.  Flight  tests  were  conducted  by 
NASA  at  Langley  Field,  and  preliminary  results  of  tests  were  released  in 
Reference  (5). 

Effect  of  forward  speed  on  vibratory  blade  stresses  and  results,  of 
correlation  are  shown  in  Figure  (4).  Here,  flight  test  data  taken  at  Sikorsky 
Aircraft  is  compared  with  calculated  values  for  the  S-58  at  speeds  above  100 
knots.  Each  test  point  shown  on  the  plot  represents  an  averaging  of  data  taken 
in  three  separate  flights.  Data  was  obtained  from  Reference  (d). 

Correlation  studies  for  a  rigid  rotor  helicopter  were  based  on  data 
furnished  by  the  Lockheed- California  Company  for  the  CL-475  helicopter 
(Reference  (1)).  Data  included  detailed  information  on  blade  etiffness  and 
mass  distribution.  Also  provided  was  a  blade  resonance  diagram,  a  curve 
showing  blade  static  deflection,  and  flight- measured  vibratory  bending 
moments  along  the  blade  at  a  number  of  airspeeds. 

Results  of  earlier  studies  (Reference  (2)  and  (3))  have  shown  good 
correlation  on  predicting  one- half  peak- to- peak  stresses  at  higher  airspeeds 
using  constant  induced  velocity.  Thus,  for  correlation  'imposes  the  highest 
airspeed  was  selected  for  which  cotnplete  CL-475  data  had  been  furnished. 

This  was  100  mph.  The  blade  was  subdivided  into  twenty  segments  in  the  fully 
coupled  flatwise-chordwise-torsiorial  analysis,  allowing  sixty  degrees-of~ 
freedom.  First  step  in  the  analysis  was  to  run  the  calculation  as  a  non¬ 
rotating  beam  to  check  analytical  stiffness  and  weight  characteristics  against 
Lockheed's  data  for  static  deflection.  Good  agreement  was  achieved  as  shown 
by  Figure  (5). 

Second  step  in  the  rigid  rotor  correlation  program  was  to  determine 
flexibility  of  spindle  bearings  at  the  blade  root.  This  was  achieved  by  varying 
the  flatwise  and  chordwise  root  springs  in  the  analysis  until  there  was  good 
agreement  with  Lockheed's  blade  resonance  diagram  as  shown  in  Figure  (6). 
From  this,  effective  springs  of  1. 75  million  inch  lbs.  per  radian  and  4. 00 
million  inch  lbs.  per  radian  were  determined  for  respecdve  flatwise  and  edge- 
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flatwise  modes  shown  by  Figure  (6)  attributed  to  motion  of  the  rotor  head.  Test 
results  had  been  obtained  by  shaking  the  hub  with  non-rotating  blades  where 
motion  of  the  hub  would  tend  to  increase  the  natural  frequency. 
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The  third  step  in  the  correlation  was  to  run  the  complete  aeroeiaetic 
analysis  at  100  mph  using  constant  inflow.  Here,  the  rotor  was  trimmed, 
cyclic  pitch  determined  and  aerodynamic  loads  calculated.  Harmonics  of 
aerodynamic  loads  were  then  applied,  and  individual  harmonics  of  blade 
response  superimposed.  The  radial  distribution  of  calculated  blade  vibratory 
bending  moments  are  shown  to  compare  favorably  with  Lockheed's  flight- 
measured  values  in  Figure  (7). 


5.  RESULTS  OF  ANALYSIS 


a.  Variation  of  Blade  Root  Stiffness 


Proper  treatment  of  blade  root  flexibility  is  of  prime  importance 
in  design  of  rigid  rotor  systems.  Rigid  rotor  designers  point  out  that 
the  term  "rigid"  is  a  misnomer  and  that  "hingeless”  is  a  preferred 
definition.  In  truth,  the  root  region  of  these  blades  is  far  from  being  a 
mathematical  cantilever.  Designed  flexibility  contributes  to  a  root 
restraint,  which  is  between  a  theoretical  rigid  and  hinged  condition. 

Effects  of  root  flexibility  on  vibratory  bending  moments  and 
stresses  are  shown  in  Figures  (8),  (9),  (10),  and  (11).  Presented  are 
results  of  calculations  for  the  Lockheed  CL-  475  helicopter  at  100  mph. 
For  the  first  three  figures,  flatwise  root  flexibility  has  been  varied, 
while  edgewise  restraint  is  held  to  a  value  representative  of  the  actual 
CL- 475  design. 

The  first  plot,  Figure  (8),  gives  the  distribution  of  flatwise 
vibratory  bending  moments  along  the  blade  for  an  articulated  version  of 
the  CL- 475.  The  curve  is  typical  of  articulated  blades  in  that  the  maxi¬ 
mum  vibratory  moment  occurs  at  about  two-thirds  blade  radius.  Maxi¬ 
mum  value  is  ±2100  inch  lbs.  The  second  plot,  Figure  (9)  shows 
corresponding  flatwise  vibratory  moments  for  a  theoretically  rigid  root 
restraint.  In  this  case  the  outboard  blade  moment  is  seen  to  increase  to 
±2700  inch  lbs. ,  while  at  the  root,  moments  reach  a  value  in  excess  of 
±8000  inch  lbs.  Finally,  shown  in  Figure  (10)  are  calculated  moments 
for  a  blade  with  selected  flatwise  root  restraint.  For  this  condition, 
outboard  blade  moments  are  reduced  to  ±650  inch  lbs.  while  root  vibra¬ 
tory  moments  increase  to  about  double  the  maximum  level  for  the 
articulated  blade.  Since  more  blade  area  is  required  at  the  root  to 
carry  centrifugal  loads,  a  section  modulus  can  be  selected  for  die  semi¬ 
rigid  blade  that  will  yield  vibratory  str  ss  levels  comparable  to  those  of 
an  articulated  design. 

Sensitivity  of  edgewise  stresses  to  proper  tuning  is  illustrated  by 
Figure  (11).  Reference  (7)  recommends  a  frequency  ratio  of  1. 4  which 
is  apparent  from  the  plot  where  vibratory  stresses  are  sewn  to  become 
excessively  large  as  the  first  edgewise  mode  of  the  blade  approaches  one 
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per  rev.  Results  given  are  for  a  neutral  c.  g.  fuselage  trim  condition. 
Center -of-gravity  shifts  can  increase  edgewise  stresses  due  to  increased 
first  harmonic  blade  bending  and  associated  Coriolis  forces. 

b.  Twist  and  Planform  Variations 


Shown  in  Figures  (12)  through  (17)  are  results  of  varying  blade 
twist  and  planform  for  the  12, 000  and  33, 000  lb.  helicopters  at  150  knots. 
The  effect  of  blade  twist  on  performance  is  given  in  Figures  (16)  and  (17). 
Only  one  power  curve  is  shown  for  both  rigid  and  articulated  helicopters 
of  the  same  gross  weight.  Calculated  power  was  the  same  for  each  type 
helicopter,  since  the  rigid  rotor  was  preconed  for  this  airspeed,  giving 
both  articulated  and  rigid  systems  the  same  coning  angle.  Power  shown 
is  that  required  by  the  rotor  and  does  not  take  into  account  tail  rotor 
power  and  accessory  drives.  Coefficients  of  drag  for  the  blade  were 
incremented  by  A CD  =.  002  to  account  for  blade  roughness.  The  non¬ 
linear  planform  blade  gave  the  b^st  performance  results  for  both  gross 
weight  machines,  but  was  the  first  to  stall  out  at  low  twists  or  if  forward 
speed  was  increased.  The  marked  increase  in  power  for  the  non-linear 
blades  from  -8  to  -4  degrees  in  Figure  (17)  is  attributed  to  a  near  stall 
condition  and  associated  high  drags.  From  Figures  (16)  and  (17)  it  can 
be  seen  that  change  in  blade  twist  can  result  in  a  maximum  power  gain 
of  about  500  hp  for  the  33, 000- ib.  helicopter,  about  200  hp  for  the 
12,000-lb.  aircraft. 

Presented  in  Figures  (12)  through  (15)  are  maximum  flatwise 
vibratory  bending  moments  versus  twist  for  the  four  blade  planforms. 
Results  are  given  for  both  rigid  and  articulated  rotor  systems.  It  is 
interesting  to  note  that  the  planform  order  from  lowest  to  highest  bending 
momenta  is  fairly  consistent  between  both  gross  weights  and  rigid  and 
articulated  systems.  For  both  gross  weight  machines  highest  blade 
bending  moments  were  for  the  relatively  stiff  3:1  planform  blade;  lowest 
were  for  the  unloaded-tip,  non-linear  blade.  The  relatively  flexible  1:2 
blade  while  showing  low  moments  would  present  root  design  problems  as 
a  rigid  rotor. 

c.  Variation  of  Blade  Stiffness 


A  requirement  of  the  aeroelastic  investigation  was  to  explore 
effects  of  blade  stiffness  on  vibratory  bending  moments.  Figures  (21), 
(22),  (23),  and  (24)  give  tesults  of  this  study  for  the  150-knot  33,000-lb. 
helicopter.  In  Figures  (21)  and  (22),  flatwise  stiffness  of  the  basic  blade 
has  been  both  doubled  and  reduced  to  one-half.  Plots  for  both  rigid  and 
articulated  blades  are  essentially  what  would  be  expected.  For  the 
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proportionately  with  increase  in  stiffness.  The  same  is  observed  for 
the  rigid  blade  in  the  root  region.  For  this  case  it  is  interesting  to  note 
that  the  point  of  minimum  bending  moment  moves  inboard  as  the  blade 
becomes  more  flexible. 
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Effects  of  linear  variation  of  stiffness  about  the  basic  blade  design 
are  given  by  Figures  (23)  and  (24).  Here,  for  both  articulated  and  rigid 
blades,  flatwise  stiffness  has  been  multiplied  by  factors  ranging  linearly 
from  2  at  the  root  to  1/2  at  the  tip  or  vice  versa.  Lowest  vibratory 
moments  are  achieved  for  the  articulated  blade  of  standard  design  as 
given  by  Figure  <,23).  Here,  stiffened  outboard  blade  section  and  softened 
cuff  appear  most  detrimental.  For  rigid  blades  (Figure  (24)),  where  the 
critical  region  is  at  the  root,  the  reverse  is  true.  Low  root  stiffness 
contributes  to  lower  moments,  and  the  stiffer  outboard  section  does  not 
have  a  significant  detrimental  effect. 

d.  Blade  Mass  Distribution 

Given  in  Figures  (18)  and  (19)  are  effects  on  vibratory  stress  for 
the  33,000-lb.  helicopter  by  separate!/  introducing  a  concentrated  36-ib. 
weight  at  five  radial  locations  along  the  blade.  Radial  stress  envelopes 
for  the  five  positions  for  an  articulated  blade  are  presented  in  Figure  (18). 
The  standard  unweighted  blade  (not  shewn)  had  a  maximum  vibratory 
stress  of  4600  p.  s.  i.  Results  shown  in  Figure  (18)  indicate:  (1)  a  signif¬ 
icant  reduction  in  vibratory  stress  by  introducing  a  weight  at  the  blade 
tip;  (2)  an  increase  in  vibratory  stress  for  weights  at  the  .  75R  and  .  50R 
positions;  and  (3)  little  effect  on  vibratory  stress  for  concentrated  weight 
added  to  blade  inboard  regions.  Effect  on  flatwise  and  edgewise  maximum 
stress  by  varying  the  amount  of  weight  added  to  the  blade  tip  is  giver,  by 
curves  shown  to  the  right  of  Figure  (20).  Here,  vibratory  stress  reduc¬ 
tion  is  greatest  with  the  first  10  lbs.  added,  then  levels  off  in  the  30-40 
lb.  region. 

Effects  of  added  concentrat-d  weights  for  rigid  blade  stresses  are 
plotted  in  Figure  (19).  In  this  case  the  root  region  is  most  critical. 
Vibratory  stress  for  the  standard  blade  (not  shown)  was  8750  p.  s.  i. 

Here,  the  tip  weight  again  produced  a  significant  reduction,  but  weights 
at  other  locations  gave  either  a  slight  decrease  or  an  increase.  The  left- 
hand  plot  of  Figure  (20)  shows  results  of  varying  tip  weight  on  rigid  blade 
stresses.  Flatwise  stress  reduction  follows  the  same  trend  as  the 
articulated  system.  The  sharp  continuous  reduction  shown  for  edgewise 
stresses  is  attributed  to  detuning  the  first  edgewise  mode.  The  rigid 
blade  in  this  case  (aspect  ratio  >*  18)  was  designed  for  ~  .  65-TL , 

Added  weight  brings  this  mode  further  below  one  per  rev.  Were  a  stiff 
inplane  design  used,  increased  added  tip  weight  would  have  the  opposite 
effect.  Edgewise  stresses  would  increase  for  the  first  edgewise  mode 
would  be  reduced  from  1.  4  XL  in  the  direction  of  one  per  rev. 

e.  Compound  Helicopters 

Interest  centers  today  in  extending  the  speed  frontier  of  helicopters 
to  the  200-300  knot  region.  The  compound  helicopter  is  a  possible  con¬ 
figuration  to  meet  this  requirement.  Here,  the  rotor  is  substantially 
unlo?deH  high-speed  flight  by  external  wings,  and  jet  or  prop-jet 
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engines  provide  the  necessary  propulsive  force.  At  the  time  of  writing 
this  paper  seme  preliminary  work  had  been  done  under  the  contract  on 
this  type  of  design.  Results  are  presented  to  shpw  early  trends,  while 
it  is  recognized  that  many  questions  still  remain  to  be  answered. 

Tne  basio  aircraft  is  defined  in  Table  I.  For  the  high-speed  flight 
condition  the  rotor  was  required  to  carry  1500- lbs.  of  lift.  Rotor 
rotational  speed  was  reduced  to  168  rpm,  which  gave  an  advancing  blade 
Kach  number  of  0.  VO  at  300  knots.  For  first  studies,  both  rigid  and 
articulated  biades  were  taken  with  -4  degrees  of  twist.  Were  the  design 
to  be  further  optimized  it  is  recognized  that  additional  stress  relief 
could  be  achieved  for  both  rotor  systems  by  a  backoff  in  blade  twist. 

An  important  variable  in  compound  design  is  rotor  angle-of-attack. 
This  can  range  from  having  the  rotor  tilt  well  forward  to  provide 
propulsive  force  as  for  the  helicopter,  to  an  aft-tilted  windmilling  rotor 
which  extracts  energy  from  the  air  Rfream.  For  the  results  shown  here 
the  rotor  has  been  placed  at  zero  angle-of-attack.  This  gives  a  very  low 
coning  angle,  blades  at  relatively  flat  pitch,  and  requires  some  power 
from  the  engines  to  drive  the  rotor. 

Since  coning  angles  were  «imi>.r,  power  and  drag  plots  are  the 
same  for  both  rigid  and  articulated  systems.  These  are  shown  in 
Figures  (23)  and  (24).  In  Figure  (23),  rotor  power  required  is  observed 
to  he  only  412  hp  at  300  knots.  The  rotor  drag  or  H- force  ploftqd  In 
Figure  (24)  is  seen  to  be  approaching  1QQ0  lbs.  at  this  airspeed.  While 
this  value  may  be  considered  hjfgh,  it  is  reduced  over  the  drag  of  g  pure 
helicopter  configuration  through  holding  rotor  angle-of-attack  flat  frnd 
blade  pitch  low,  : 

Some  interesting  trends  are  noted  in  the  plot  showing  flatwise 
vibratory  bending  moments  versus  airspeed  (Figure  (2$)).  Maximum 
vibratory  momenta  are  seen  to  increase  more  sharply  with  airspeed  for 
the  rigid  than  for  the  articulated  case.  This  trend  is  partially  explained 
by  the  diagrams  given  on  the  plot  showing  spanwise  disjbrjbiiticin  of 
bendi.g  moments  for  both  rotor  systems  at  200  and  3CjQ  knots,  jiere,  it 
c  an  be  seen  that  outboard  moment  envelopes  are  similar  for  both  systems 
and  build  up  at  the  same  rate  with  airspeed.  However^'  ropt  moments, 
which  govern  rigid  blade  design,  are  seen  to  digress  r^iore  rapidly. 


6.  CONCLUSIONS 

. . 

gar»y  results  have  been  presented  from  an  analytical  program  to  explore 
ffev A  parametric  variations  on  blade  stress  and  performance  of. rigid  and 
.  m*.  uiared  rotor  systems.  The  analysis  considers  the  high-speed  condition 
in  to  account  the  fully  coupled  flaiwise-edgewise-torsioncl  response 
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of  the  rotor  blades.  For  the  studies  constant  inflow  was  taken.  Conclusions 
which  may  be  drawn  from  the  results  include: 


a.  Calculated  performance  and  one-half  peak-to-peak  stresses  are 
shown  to  give  good  agreement  with  flight  measured  values.  Good 
stress  correlation  is  shown  for  both  rigid  and  articulated  rotor 
helicopters  at  higher  airspeeds.  Flight  measured  data  is  taken 
for  the  Sikorsky  S-58  and  Lockheed  CL- 475  helicopters. 

b.  Rigid  rotor  flatwise  and  edgewise  vibratory  stresses  are  observed 
to  be  highly  sensitive  to  design  of  the  blade  root  region.  For  this 
rotor  system  properly  designed  flexibility  gives  a  root  restraint, 
which  is  between  a  theoretical  rigid  and  hinged  condition. 

c.  For  the  articulated  blade,  maximum  flatwise  vibratory  stress 
generally  occurs  at  about  two-thirds  blade  radius.  Location  will 
vary  with  planform,  twist,  and  whether  the  aircraft  is  a  helicopter 
or  compound.  For  the  rigid  rotor,  maximum  flatwise  vibratory 
stress  occurred  at  the  blade  root  for  all  cases  studied. 

d.  Results  of  calculations  for  planform- twist  variation  indicate  a  rapid 
rise  in  vibratory  bending  moments  with  increase  in  blade  twist  for 
both  rigid  and  articulated  rotor  systems.  Highest  vibratory  bend¬ 
ing  moments  were  noted  for  the  3:1  tapered  blade. 

e.  Studies  on  variation  of  blade  stiffness  indicate  that  outboard  blade 
stiffening  is  detrimental  for  an  articulated  blade,  whereas  inboard 
blade  stiffening  shows  the  largest  increase  in  vibratory  moments 
for  a  rigid  blade. 

f.  Reduction  ir.  vibratory  bending  moments  was  observed  for  both 
rigid  and  articulated  blades  by  introducing  a  concentrated  weight 
at  the  blade  tip.  For  the  33,000-lb.  helicopter  considered, 
greatest  reduction  was  achieved  for  both  systems  with  the  first 
10  lbs.  added,  lesser  reduction  thereafter. 

g.  Preliminary  analysis  for  compound  helicopters  indicates  a  more 
marked  increase  in  flatwise  vibratory  bending  moments  with  air¬ 
speed  for  the  rigid  than  for  the  articulated  rotor  system. 
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Aerodynamic 
Angles  of  Attack 


trim  Ik 


Helicopter  Basic  Blade  Data 
Gro*§  Weight  33, 000  lbs. ,  150-Knot  Deaign 


Articulated  Blade 


Radii 

Weight  Per 

Blade 

Blade 

MidUne 

iM 

Segment  (Lb. ) 

Twist  (Deg. ) 

Chord  (In.) 

Area  (Sq.  In. ) 

432,0. 

6.97 

-1.79 

23.65 

9.55 

410;  4 

19.60 

-1.52 

23.65 

19.20 

3SM 

14.51 

-1.22 

23.65 

19.40 

367.2 

15. 77 

-0,91 

23.65 

19.73 

345.6 

15.74 

-0.74 

23.65 

19.86. 

324.0 

17. 14 

-0.30 

23.65 

20,12 

302.4 

10.29 

0.00 

23.65 

20. 15 

280.8 

17.63 

,30 

23.65 

20.20 

259.2 

17.02 

.60 

23.65 

20.23 

237.6 

18.04 

.91 

23.  £5 

20.28 

216.0 

17.71 

1.2 

23.65 

20.30 

194. 4 

16. 61 

1.51 

23.65 

20;  36 

172.8 

16. 42 

1.81 

23.65 

20.40 

151.2 

16.42 

2.13 

23.65 

20.43 

i29. 6 

17.33 

2.39 

23.65 

20.57 

108.0 

18.46 

2.59 

9;  71 

21.00 

£6.4 

38.11 

2.59 

9.71 

24.60 

64.8 

124.62 

2.59 

9.71 

100.00 

43.2 

77. 12  , 

2.59 

9.71, 

100.00 

24.0 

84.59 

2.59 

9.71 

100.00 
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1 

8 

I 

8 

0 

B 

0 

0 

0 

Q 

0 


AHBMDDC I  (count) 
Articulated  Blade 


‘  s‘  ■> , 


Spar  Wall 
Thickness  (In. ) 

Flatwise 
Ixx  (In^) 

.  Edgewise 

iyy  (toft 

Torsion 

M 

Section  * 
Zxx  (In?) 

modulus 

mMl 

,081 

1.87 

19.9 

5.90 

1*38 

3.8 

!  162 

3.75 

39.8 

12.70 

2.92 

7*8 

.183 

4.20 

42.9 

15.00 

3.  $4 

8.5 

,204 

4.75 

53.0 

14.90 

3.62 

9.8 

.212 

5.05 

55.3 

16.55 

3.80 

10.2 

.216 

5.42 

64.8 

18*00 

4 .  15 

11.3 

.$21  - 

5.58 

65.3 

19.05 

4.20 

11.5 

.224 

5.68 

66.0 

19.30 

4.28 

11.6 

.221 

5.80 

66.8 

19.70 

4.35 

11.7 

.232 

5.95 

67.3 

20.03 

4.40 

11.9 

.237 

6.05 

68.0 

20.30 

4.48 

12.0 

.241 

6.18 

68.8 

20.75 

4.55 

124 

.246 

6.25 

69.3 

21.05 

4.60 

12.3 

.249 

6.40 

70,1 

21.40 

4.7Q 

12.4 

.258 

6.45 

72.0 

21.85 

4.9Q 

12.8 

.290 

8. 10 

80.8 

23.40 

5.80 

14.3 

.375 

24. 10 

124.0 

48.80 

13*20 

22.2 

3.0 

100.00 

100.0 

100.00 

25.00 

25.0 

3.0 

100.00 

100.0 

150.00 

25.0 

25.0 

3.0 

100.00 

100.0 

150.00 

is.o 

25.0 

0 

B 

0 
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APIPBNDIX  Kcontd;) 


R*g*d  Blade 


Radii 

Weight  Per 

Blade . 

Blade 

Midline 

(in.) 

Segment  (Lb. ) 

Twist  (Deg, ) 

Chord  (In. ) 

Are*  (Sq.  In. ) 

432.0 

6.97 

0. 

23.65 

9.55 

410.4 

19.60 

0. 

23.65 

19.20 

388.  $ 

14.51 

0. 

23.65 

19.40 

367.2 

15.77 

0. 

23.65 

19.73 

345.6 

15.74 

0. 

23.65 

19.86 

324.0 

17. 14 

0. 

23.65 

20. 12 

302.4 

10.29 

0. 

23.65 

20. 15 

280.8 

17.63 

0. 

23.65 

20.20 

259.2 

17;  02 

0. 

23.65 

20.23 

237.6 

18.04 

0. 

23.65 

20.28 

216.0 

17.71 

0. 

23.65 

20.30 

194.4 

16.61 

0. 

23.65 

20.36 

172. 8 

16.42 

0. 

23.65 

20;  40 

151.2 

16.42 

0. 

23.65 

20.43 

129.6 

17.33 

0. 

23.65 

20.57 

108.0 

18.46 

0. 

9.71 

21.00 

•M 

38.11 

0. 

9.71 

24.60 

64.8 

124.62 

0. 

9.71 

100.00 

43.2 

77. 12 

0. 

9.71 

100.00 

24.0 

84.59 

0. 

9.71 

100.00 

1$ 


I 

I 


"  */ 
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*  *» j  Vt.,  »”  MgS 
^a'  iv.  £»  .,v£ 

"r  1  .  *is; 


\  V  ‘ 

*?.;  ivy 
“  v  5$| 


I 

I 

I 

I 

I 

I 

I 

K 

K 


Rigid  Blade 


Spar  Wall 

Flatwise 

Edgewise 

Torsion 

Section 

Modulus 

Thickness  (In. ) 

Ixx  (In*) 

JffJjaa 

1  «"■*) 

Sex  (In?) 

gyy<^» 

.081 

1.87 

19.9 

5.9 

1.38 

3.8 

.162 

3.75 

39.8 

12.7 

2.92 

7.8 

.183 

4.20 

42.9 

15.0 

3.24 

8.5 

.204 

4. 75 

53.0 

14.9 

3.62 

.9.8  . 

.212 

5.05 

55.3 

16.6 

3.80 

10.2 

.216 

5. 42 

64.8 

18.0 

4. 15 

11.3 

.221 

5.58 

65.3 

19.1 

4.20 

11.5 

.224 

5.68 

66.0 

19.3 

4.28 

11.6 

.228 

5.80 

66.8 

19.7 

4.35 

11.7 

.232 

5.95 

67.3 

20.0 

4.40 

11..  9 

.237 

6.05 

68.0 

20.3 

4. 48 

12.0 

.241 

6.18 

68.8 

20.8 

4.55 

ill 

.245 

6.25 

69.3 

21.1 

4.60 

12.3 

.249 

6.40 

70.1 

21.4 

4.70 

12.3 

.258 

6.45 

72.0 

21.9 

4.90 

118 

.290 

8. 10 

80.8 

23.4 

5.  $Q 

14.3 

.375 

10.00 

124.0 

48.8 

6.00 

212 

3.0 

10.00 

43.0 

100.0 

7.00 

10.0 

3.0 

1Q.00 

43.0 

150.0 

7.00 

10.0 

3.0 

10.00 

43.0 

150.0 

7.00 

10.0 

I 


I 


I 

8 

1 

I 

I 

I 
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.L 
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Q-IOCKHEKD 


/ 


MADE  STATION  (INCHES) 


mm  5 
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I 

E 

I 

i 

K 

6 
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vs 
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I 

I 
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IN-FLIGHT  MEASUREMENT  OF  STEADY  AND 
OSCILLATORY  ROTOR  SHAFT  LOADS 
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Vortol  Division,  Booing  Aircraft 


IN-FLIGHT  MEASUREMENT  OF  STEADY  AND  OSCILLATORY 
ROTOR  SHAFT  LOADS 


Introduction 

A  set  of  comprehensive  tost  data  on  alternating  shaft  loads  Is  an  obvious 
necessity  In  the  development  of  improved  methods  for  calculating  such  loads. 

A  program  of  in-flight  shaft  load  measurement  was  conducted  as  part  of  an 
overall  vibration  research  effort  sponsored  by  the  Amy  and  Air  Force  under 
Contract  AF  33(616) -5240  and  published  as  ASD-TDR-62-284.  This  paper  describes 
the  load  measurement  system  and  the  results  obtained  from  the  flight  tests. 

Load  Measurement  System 

This  development  of  a  shaft  lead  measurement  system  represented  a  culmination 
of  several  previous  attempts  with  very  limited  success.  However,  from  the 
results  of  the  previous  pioneering  attempts  it  was  evident  that  the  beet 
approach  lay  in  further  development  of  the  shaft  strain  meaaurament  system 
which  had  been  used  previously.  Additional  work  indicated  that  through  the 
use  of  (1)  strain  gages  with  high  sensitivity,  (2)  high  voltage  in  the  strain 
gaga  bridges,  (3)  an  improved  bridge  network,  end  (4)  separate  load  measuring 
systems  for  the  eteady  end  alternating  components  of  both  lift  end  torque, 
the  desired  sensitivity  and  accuracies  could  be  achieved. 

An  H-21  helicopter  wee  equipped  with  strain  gages,  as  illustrated  in  Figure  1, 
to  measure  the  following  loads. 

Foxwerd  end  eft  transmission  end  rotor  shafts: 

1.  Fg,  shear  in  the  rotor  plane  in  the  direction  of  the  master  spline. 

2.  Fy,  shear  in  the  rotor  plane  in  the  direction  normal  to  the  master  spline. 

3.  Mx,  moment  in  the  rotor  plane,  about  the  rotating  R  axle. 

4.  My,  moment  in  tfca  rotor  plane,  about  the  rotating  Y  axis. 

5.  ?z>  lift  force  in  the  shaft  vertical  direction. 

6.  Mg,  torque  about  the  shaft  axis. 

A  developed  view  of  the  rotor  shaft  strain  gage  instrumentation  is  shown  in 
Figure  2. 

Ifil&JBftyfemiga 

Load  calibration  of  the  instrument ed  forward  and  aft  rotor  shaft  was  per¬ 
formed  in  a  specially  constructed  fixture,  and  known  loads  were  applied 
individually  et  the  hub  in  ell  directions.  Incremental  loads  were  applied 


MASTER  SPUME 
CPosmoM  or  red  blade) 


to  tha  shaft  and  the  output  of  all  gagoe  to  such  load  was  r«cocd«d>  regardless 
of  tha  direction  of  tha  loada^ln  order  to  obtain  Interaction  aa  vail  as  primary 
load  calibration  factors. 

JBESm  -  Initially,  tha  torqua  calibration  uaed  a  load  balance  systan  as  sbovn 
in  tba  lower  half  of  Figure  3,  but  this  ayatea  was  replaced  by  a  closed  loop 
ayatea  to  improve  tha  accuracy.  Final  torque  calibration  was  perforated  uei.ig 
the  cast  setup  shown  in  tha  upper  half  of  Figure  3.  A  schematic  of  this  loading 
system  shoeing  tha  closad  loop  cable  arrangement  is  given  below. 


FIXTURE 


The  torque  calibration  was  performed  usirg  the  tumbuckle  to  apply  tha  lncre> 
•dotal  cable  tension  in  eonbinatlon  with  a  load  cell  for  measuring  the  applied 
load. 


MU  -  Uft  calibration  loads  ware  applied  with  tha  hydraulic  jack  shown  la 
figure  4.  Accurate  altgnnnt  of  the  load  through  tha  canter  of  tha  shaft 
was  daeahiianed  using  cn' optical  myiimm  consisting  of  s  sighting  scops  and 
target.  X*  tha  calibration,  tha  incremental  loads  applied  to  tha  shaft  wars 
in  addition  to  the  swan  Uft  load  sinulatlng  the  aircraft  weight. 


!|JV 


TORQUE  LOADING 

CLOSED  LOOP  CABLE  LOAD  SYSTEM 


TORQUE  LOADING 

EXTERNAL  WEIGHT  COUPLE  LOAD  SYSTEM 
FIGURE  3,  ROTOR  SHAFT  CALIBRATION 
TORQUE  LOADING 


LIFT  LOADING 
LIFT  LOAD  AT  ROTOR  SHAFT 
FIGURE  4.  ROTOR  SHAFT  CALIBRATION 
LIFT  LOADING 


-3fr**r  *  Application  of  shear  loads  for  calibration  of  the  strain  gage  bridges 
vas  performed  as  shown  in  Figure  5.  Rotor  shaft  shear  loads  were  applied 
using  a  hydraulic  jack  Mounted  at  its  base  on  a  ball* joint  to  insure  loading 
through  the  center  line  of  the  shaft.  In  order  to  obtain  calibration  curves 
for  both  of  the  shear  bridges  used  in  the  test,  the  shaft  was  loaded  on  the  0°- 
180°  axis  in  addition  to  the  90°-270°  axis. 

Honest  -  Calibration  loads  for  bending  acewnt  were  imposed  on  the  rotor  shaft 
using  the  test  arrangement  shown  below, 


Using  the  load  calls  located  at  each  and  of  tha  meant  arm  permitted  one  side 
to  be  loaded  in  tension  while  the  other  was  In  compression,  thus  producing  a 
pure  moment  at  tha  centerline  of  tha  shaft.  Calibrations  ware  performed  for 
the  Dm  banding  bridges  by  applying  equal  tension  and  compression  loads  along 
the  ails;  the  Mgg,  bending  bridges  were  calibrated  in  the  sms  msnuer  after 
tha  shaft  had  bean  rotated  90°. 


SHEAR  LOADING 

HYDRAULIC  JACK  LOAD  SYSTEM  AND  LOAD  CELL 


SHEAR  LOADING 
INCLINOMETER  CHECK 
FIGURE  5.  ROTOR  SHAFT  CALIBRATION 
SHEAR  LOADING 


■vwuxmw  -act— ru> <i  mnei  ammmmmiamm  —  r  •  —  -rimfi  'if  nwuWtit  ftiitGMffiSWBr' 


Load  Interaction  -  A  typical  calibration  chart  is  illustrated  in  Figure  6  for 
shear  load.  The  calibration  constant  for  shear  load  (9Q°-270°)  is  obtained  by 
dividing  the  load  variation  (1200#  -  0)  by  the  trace  deflection  (364-210)  counts. 

K  -  x  131  eta  -  1040  #/in. 

154  cts  in,  trace 

This  chart  indicates  that,  as  intended,  the  predominant  gage  response  occurs 
between  the  applied  shear  load  and  the  strain  gage  bridge  representing  shear. 

This  is  the  pair  of  bending  gages  which  are  wired  together  to  measure  dif¬ 
ferential  moment,  and  designated  as  a  shear  bridge.  It  is  also  apparent  that 
other  gages  respond  which  are  not  in  the  direction  of  the  calibration  load. 

Such  aide  effect  may  be  due  to  an  unfavorable  location  of  a  particular  gage 
with  respect  to  a  load  direction  it  is  not  Intended  to  measure,  or  due  to 
Poisson's  ratio  effect.  These  side  effects  are  accounted  for  by  an  "Interaction 
Matrix"  which  also  contains  the  main  load  calibration.  This  matrix  is  the  com¬ 
plete  relationship  between  the  measurement  and  the  loads  on  the  shaft  in  all 
directions  at  the  rotor  plane.  This  result  is  an  Interaction  matrix  of  the  form. 


Fx 

iax  4ay  1a*  iac*  *ajf  ^a* 

FX 

Fy 

*bx  *by  *bx  ibmt  ^bp 

Fy 

h 

TSt 

*xx  *cy  *co< 

• 

FX 

Hex 

^dx  ldy  id* 

Ho, 

M* 

*ex  *ay 

M* 

lfx 

M* 

Fy,  Fy,  Fz  ■  actual  loads  in  the  X,  Y,  Z  direction 
My,  M2  *  actual  moments  about  the  X,Y,  Z  axes 
*ax»  iay>  *as  "  interaction  of  gage  "a"  to  unit  X,  Y,  Z  forces 
*aw»  iag»  *ag  m  interaction  of  gage  "a"  to  unit moments 
Fxt  Fy>  h  *  apparent  forces  from  the  strain  measurements 
Me, ,  Mjj,  Mg  -  apparent  moments  from  the  strain  measurements 


The  numerical  interaction  matrices  from  the  load  calibration  are  shown  below. 


ntGQBB  6 


Steady  Loads 


1  +0.0633  +0.0021  +0.0001  -0.0089  -0.0016 

-0.1171  1  +0.0010  +0.0067  -0.0030  +0.0006 


-0.0379  +0.0416 


+0.0015  +0,0091  -0.0151 


0  -9.517  -0,0124  1  -0.0104  +0.0005 

+9.415  -0.0256  +0.0881  +0.0156  1  +0.0074 


L"*J 

Measured 


-0.0658  -0.0144  -0.0068  -0.0205  -0.0331  1 


+0.0633  +0.0021  +0.0001  -0.0089  -0.0016 


-0.1171  1  +0.0010  +0.006/  -0.0030  +0.0006 


+0.0583  +0.1008 


+0.0052  -0.0088  +0.0029 


-9.517  -0.0124  1  -0.0104  +0.0005 


+9.415  -0.0256  +0.0881  +0.0156  1  +0.0074 


L*J 

Measured 


-0.1567  +0.1142  -0.0343  -0.0289  -0.0467 
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Actual 


It  la  acted  that  the  interaction  Matrices  for  the  steady  and  alternating 
loads  are  coapoaad  of  the  seas  eleaents  in  the  rove  and  coluans  Involving 
interaction  between  tba  she«r  and  bending  ament.  However,  the  e levants 
Mandated  with  lift  and  torque  loads  usually  exhibit  unlike  eleaents, 
fjtikttly  resulting  from  the  nonlinear  load  characteristics  existing  between 
tie  steady  and  alternating  load  ranges  used  in  the  calibration,  and  partly 
because  independent  strain  gaga  bridges  were  used  for  the  steady  and 
alternating  lift  and  torque  loads. 

IsttrumttatlMi  Evaluation 

■  - "  . . . — 

1  -  following  the  load  calibration,  the  eft  treneaisslon  was 
v ,  voter  stand  with  a  standard  H-21  hub,  wood  rotor  blade,  and 
upper  eoatrol  ays tea.  the  purpose  of  this  test  was  to  evaluate  the 
Uty  end  accuracy  of  the  aystea  under  operating  conditions. 


With  the  full  strain  recording  system  operative,  the  rotor  was  whirled  through 
a  full  range  of  rotor  speeds  from  240  to  278  rpo,  varying  collective  pitch 
setting  fxoa  flat  pitch  to  12° ,  with  three  cyclic  pitch  settings.  The  data 
obtained  from  this  test  showed  that  the  strain  gage  system  was  operating  satis¬ 
factorily. 

Accuracy  gstlaate  -  The  results  of  the  pr'^ery  load  calibration  are  tabulated 
below,  along  with  an  estiaate  of  the  load  trough-record  accuracy  based  on  an 
assuaed  oscillograph  trace  reading  error  of  +  0.02  inches. 

Estimated 


Item 

?X 

Trace 

Calibration 

Factor 

Reading 

Accuracy 

9  +  0.02" 

Estimated 

Overall 

Accuracy 

Typical 

Load 

Value 

X 

Accuracy 

1,000  lb/in. 

± 

20  lb 

± 

36  lb 

1,000  lb 

4 

fy 

1,115  xo/in. 

+ 

22  lb 

± 

36  lb 

1,000  lb 

k 

*pz 

6,050  lb/in. 

± 

120  lb 

± 

360  lb 

7,000  lb 

5 

Steady 

pz 

Alt. 

721  lb/in. 

± 

20  lb 

+ 

30  lb 

500  lb 

6 

Mx 

22,050 

in. 

± 

440  in. lb 

+ 

500  in. lb 

5,000  in. lb 

10 

19,290  iHtlfe 
in. 

± 

400  in. lb 

+ 

450  in. lb 

5,000  in. lb 

9 

“z 

Steady 

92,500  ABult 
in. 

+ 

1,860  in. lb 

± 

2,400  in. lb 

100,000  in. lb 

3 

*2 

14,710  iBoiil 
in. 

+ 

300  in. lb 

350  in  lb 

5,000  in. lb 

7 
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As  noted  above,  strain  gage  bridges  on  the  forward  and  aft  rotor  shafts  asasure 
axial  shaft  loads  and  torque,  and  shear  and  bending  moments  in  two  directions.  A 
square  matrix  of  ituaerical  coefficients,  froa  the  load  calibration,  is  then  applied 
to  the  data  to  account  for  interaction,  i.e.,  apparent  response  of  a  load  aaasuring 
device  to  loads  it  should  not  theoretically  respond  to.  The  corrected  loede  in  the 
rotatina  shaft  are  next  transferred  into  esutw*!***  losdt  in  the  fixed  eecrdisatt 
•paten  of  the  fuselage . 


Figure  7  presence  e  flew  diagram  Illustrating  the  steps  used  In  converting  the 
■eaaured  strain  dace  to  the  correct  Fourier  coefficients  for  the  steady,  and 
first  four  harmonies  of  rotating  shaft  loads .  Measured  flight  data  whs  har¬ 
monically  analysed  using  a  24  ordinate  analysis  ptogrsaiacd  on  a  digital  computer 
to  obtain  the  steady  end  harmonic  losd  component*.  Following  the  harmonic 
analysis ,  each  steady  tar*  of  tha  aarlas  was  adjuatsd  to  corract  for  the 
difference  between  the  base  line  ueed  In  the  harmonic  analysis  and  the  aero  load 
position  from  the. flight  calibrations. 

Continuing  to  tbs  right  on  ths  flow  diagram,  the  hernonlc  components  srs  con¬ 
verted  to  shaft  loads  by  applying  ths  aaplitude  and*  in  too*  cases*  phast 
calibration  constants  obtalnsd  from  tha  load  calibration  teat.  Following 
these  corrections*  tha  matured  data  la  ragrouped  into  steady  end  alternating 
tame  in  Item  1-4,  alternating  term  only  In  Item  5  end  6,  and  eteedy 
coefficients  In  I'.eas  7  end  8.  The  nest  step  in  the  processing  consists  of  an 
intaracelon  corrtction  using  eha  lntaractlon  Matrices  obtained  during  the  load 
calibration.  After  the  inversion  of  tha  staady  and  alternating  Interaction 
Matrices,  tha  corrtctad  coefficients  for  steady  and  alternating  loads  can  ba 
obtalnsd  by  adjustment  of  tha  Fourlar  coafflcients  using, 
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ifa 
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4«r 

Actual  Inversion  of  lntaractlon  Matrix 

Loads 


*Y 

Mm 

* 

Measured 

Loads 


The  rotor  system  la  rsprasantad  by  tha  schematic  diagram  of  Flgura  •  which 
•how  the  rotating  end  fixed  axes  system.  Fixed  system  loads  can  he  wrifcte** 
as  trigonometric  functions  of  tha  relative  angle  At  between  the  rotting 
axis  aystam  and  fixed  exit  system  as  shown  in  Figure  9.  forward  rotor  eteedy, 
first  harmonic  and  third  harmonic  fixed  system  coefficients  in  term  of 
rotating  system  coefficients  are  shown  in  Figure  10.  Similar  expressions  fur 
tha  fixed  system  shaft  loads  of  the  eft  rotor  derived  using  e  negative  direetlen 
of  rotation  ere  included  in  the  date  processing  computer  program* 


Shaar  In  tha  rotor  plana  in  tha  direction  of  the  neater  spline 
Fyg,  Shaar  in  tha  rotor  plane  In  a  direction  normal  to  tha  master  spline 
*ZR*  hiit  *orc*  in  abaft  vertical  direction 
Mqi  Moment  in  tha  rotor  plana,  about  the  rotating  x  axis 

Myg,  Moment  in  tha  rotor  plana,  about  the  rotating  y  axis 

Mjy,  Torque  about  tha  shaft  axis 

Pgy,  Shear  in  tha  rotor  plana  In  the  aft  direction 
Fyp,  Shear  in  the  rotor  plana  in  tha  lateral  direction 
Vgf,  Lift  force  in  the  shaft  vertical  direction 
Hgy,  Moment  in  the  rotor  plana,  about  tha  fixed  x  axis 

Hyp,  Moment  in  the  rotor  plane,  about  the  fixed  y  axis 

Hgpir  Torque  about  tha  shaft  axis 

FIGURE  *.  MSAStHf)  F1XH>  STS  TIM  LMDS 


i 

JdMttA 

*01  *  1/iAlX  “  1/2»1T 

*01  “  */2BlX  +  l/2AlT 

First  Subscript  Denote*  Rotor 
Harnoaic 

Second  Subscript  Denotes  Load 
Direction  at  Shown  Belov, 

*552  *  *OZ 

Xs  Force  Along  X  Axis 

*0*  "  1^2A1CX  '  l/2*tjg 

Y,  Force  Along  Y  Axis 

*Op  "  1/f2Bia  +  \Hk\p 

Z,  Force  Along  Z  Axle 

* 0 1  “ 

<*,  Moment  Akiut  X  Axis 

(Fixed)  (Rotating) 

Moment  About  Y  Axia 

r»  Moment  About  Z  Axis 

Hat  Hwanie 

SaLJlflML 

Sin  Term 

A1X  *  *0X  *  */2A2X  ’  *^2B2Y 

i  * 

B1X 

-  -Aoy  +  1/2B2X  +  1/2A2Y 

A1Y  “  Hy  +  1/2a2Y  +  1^2B2X 

*1Y 

■  A q,  -  1/2A2x  +  i/2B2Y 

A1Z  "  A1Z 

biz 

"  B1Z 

Alc*  "  *001  +  i/2A2o<  *  1/2b2 fi 

Blct 

-  -*qu  +  l/2B2o (  +  UlKifi 

*lp  *  A00t  +  1/2B20C  +  l/2A2Jf 

Bi£ 

-  *0  J3  -  1/2A20C  + 

Al*  "  AU 

Bir 

•  Bir 

(Fixed)  (Rotating) 

(Fixed) 

(Rotating) 

IMid  ftttlMdUL 

fiat  Him 

Sin  Terns 

Ajjj  -  l/2Ajx  +  1/282y  +  1/2*4,  '  1/2*4Y 

*3*  •  1/2*2X  -  1/2A2Y  ♦  1/2*4, 

a3Y  ■  *1/2*2,  +  1/2A2y  *  1/2*4,  + 

1/2*4* 

+  I/2A4 y 

*3Y  "  X  *  l^2B2Y  “  l/2*AX 

*3Z  "  a3Z 

+  1/2R4Y 

a**  -  i/2A**  ♦  i/2%»  +  l/2a*w 

.17211. «. 

• 

Hz  “  Hz 

*3^  T  +  V**2fi  ♦  l/2»4Ct 

+  Ut**fi 

H*  *  1/2***.  +*/*Hot 

Hr  "  Hr 

ha  *  «'»!»  -*/»*** 

(Fixed)  (Rotating) 

Hr  "  Hr 

(Fixed)  (Rotating) 

wmu  10.  rtm  mm  cmnmmn 


Rotor  Shaft  Load  Re/n'li » 


Using  th«  measurement  technique  described  above,  H-21  rotor  loads  were  obtained 
et  a  normal  gross  weight  of  13,500  lb  for  (1)  an  airspeed  sweep  at  normal  258 
rpm,  (2)  a  rotor  sweep  at  40  knots  airspeed  and  (3)  a  rotor  speed  sweep  at  90 
knots  airspeed.  Load  plots  are  shown  herein  for  the  steady  and  third  harmonic 
fixed  system  results,  and  second  fourth  harmonic  rotating  results.  Each 
*et  of  force  data  contains  forward  and  aft  rotor  loads  as  sine  and  cosine  com- 
‘  ments  and  their  resultants.  The  sine  and  cosine  time  reference  is  when  the 
rjd  blade  and  the  rotor  hub  master  spline  are  in  trail  position. 

■Steady  Loads.  Instrumentation  Check  -  Ho  rational  method  was  known  for  positively 
defining  the  accuracy  of  the  measured  loads.  Apart  from  an  examination  of  the 
reasonableness  of  the  flight  data  compared  with  calculated  shaft  loads,  the 
best  check  of  the  flight  data  was  a  steady  load  balance.  This  provides  a 
gross  check,  but  does  not,  of  course,  provide  a  check  down  into  the  reading 
accuracy  of  the  data.  Even  this  is  appropriate  because,  in  addition  to  the 
measured  loads  and  the  aircraft  weight,  other  unknown  loads  act  on  the  heli¬ 
copter,  namely,  rotor  downwash  and  aerodynamic  drag.  Further,  the  actual  eg 
location  at  any  Instant  during  a  flight  when  data  is  recorded  is  not  known 
exactly  but  is  estimated  from  the  takeoff  eg  location. 

To  illustrate,  a  steady  load  balance  is  performed  below  using  measured  data 
from  a  hover  condition  In  which  the  unknown  factors  are  minimised. 


FLT.  C96  X-121  RPM  -  258  HOVER 


RUM  1 

G.W.  13,500  lb 


fitsL 

FX  *  +476# 

Fy  **  +5# 
p z  -  +7691# 

Mx  «  +2268  in-# 
My  -  +8821  in-# 
MZ  -  -124001  in-# 


For  Che  gross  weight  end  eg  condition  shown  end  assuming  vertical  trim  and  no 
wind  the  rotor  shafts  are  1.5°  forward  of  the  vertical  in  the  pitch  direction. 
For  force  balance  in  either  the  lateral  or  longitudinal  directions,  the 
appropriate  forces  from  each  rotor  are  combined  to  obtain  the  total  force  acting 
on  the  helicopter  in  each  direction  which  should,  of  course,  be  xero.  In  the 
vertical  direction,  the  total  should  equal  the  helicopter  given  weight  plus 
any  download  due  to  rotor  downwash  impinging  on  the  fuselage  and  pitch  link 
compressive  loads  but  minus  the  rotor  weights. 

Summing  the  loads  in  the  lateral  direction, 

5IFY  "  PY  Fwd  +  FY  Aft 

]>“  Fy  -  5#  +  28#  -  33# 

In  the  longitudinal  direction,  the  components  of  lift  and  shear  are  combined 
as  illustrated  below: 


A&l 

Fx  -  -174# 

Fy  *=  +2C# 
pz  -  +5569# 

Mx  «  +1103  in-# 
My  *  -1451  in-# 

“  +122866  in-# 


X'Lonj  ‘  <FW  +  ftA>  +  »**>  c®*  15° 

-  (7691  +  5569). 0262  -(476  -  176). 9997 
•  +44# 


Similarly,  for  the  vertical  force  at  each  rotor 

FVert  *  FZP  cot  1.5°  +  FXF  «in  1.5° 

FWd  7691 (.9997)  +  476 (.0262)  -  7703# 

FVert  -  fza  coa  1.5°  +  Fxa  sin  1.5° 

Aft  5569( .9997)  -  174(.0262)  -  5574# 

Further,  coe&ining  the  vertical  loads  at  each  rotor, 

XFVert  "  7703  +  5574  "  13 * 277 

Before  conaidering  the  comparison  of  measured  forcea  and  known  forcea  acting 
on  the  helicopter,  the  unknown  lift  load  at  each  rotor  la  eatimated  from  the 
eg  location. 


|Lp 

L - 254" 


♦ 


13,500# 


La 


Forward  Rotor, 
508  5 


Aft  Rotor, 

&&--■  21 


UA 


508 


13,500 

13,500# 


7350# 

6150# 


From  this  distribution,  the  ateady  lifts  are  defined  by  adding  the  estimated 
downwash  and  pitch  link  compressive  load,  and  then  subtracting  the  rotor  weight. 


Lp  -  7350  +  13,500  (.02) 
Eat.  Downwash 

La  -  6200  +  13,500  (.02) 


+  3(70)  -  667  *  7163# 

Fitch  Link  Rotor 

+  3(70) 


667  -  6003# 


The  table  below  presents  a  comparison  of  the  measured  forces  and  the  estimated 
steady  loads: 


Forward  Rotor  Aft  Rotor  Forward  &  Aft  Rotor 

Load 

Direction  Measured  Est.  Error  Measured  Eat.  Error  Measured  Est.  Error 
Longitudinal  -  -  440  +44 

Lateral  ...  33  0  +33 

Vertical  7703  7163  -540  5574  6003  +429  13169  132-»7  +108 


Lateral  and  longitudinal  loads  are  close  to  the  estimated  accuracy  shown 
previously.  Further,  the  difference  between  the  estimated  total  vertical 
lead  and  the  measured  loads  are  within  estimated  accuracy  of  360  lb. 
Individually,  the  forward  and  aft  rotors  are  slightly  above  the  specified 
accuracy,  but  this  variation  reflects  some  inaccuracy  resulting  from 
*s':laation  of  the  eg  location  and  downvash  distribution. 

From  the  measured  torque  the  shaft  horsepower  is  calculated  as: 

Fwd  Rotor  508  H.P. 

Aft  Rotor  503  H.P. 

Total  Shaft  Power  1011  H.P. 

This  measured  power  shows  reasonable  agreement  with  a  calculated  power 
requirement  of  1100  H.P.  for  the  same  flight  condition. 

From  the  results  of  this  steady  load  check,  it  was  concluded  that  the 
accuracy  of  the  steady  load  measurement  was  adequate  for  the  flight  program. 

Fixed  System  Steady  Loads  -  Figure  11  presents  the  fixed  system  hub  forces 
and  moment  as  an  airspeed  sweep  at  258  rpm  and  rpm  sweeps  for  40  and  90 
knots. 

For  reference,  the  H-21  forward  and  aft  rotor  shafts  are  parallel,  but 
the  rotor  is  rigged  to  produce  a  dihedral  effect  which  promotes  better 
speed  stability.  A  side  effect  of  this  feature  is  that  in  order  to  maintain 
a  hover  attitude,  the  combination  of  longitudinal  cyclic  and  differential 
collective  pitch  control  must  provi.’ ».  a  greater  aft  longitudinal  force 
at  the  forward  rotor  than  the  aft  rotor.  This  is  illustrated  in  Figure  11a 
showing  a  500  lb  aft  force  at  the  forward  rotor,  and  a  200  lb  forward  force 
at  the  aft  rotor  In  hover.  As  speed  increases,  the  forward  force  at  the 
aft  rotor  decreases,  and  above  80  knots  reverses  direction.  The  aft  longi¬ 
tudinal  force  at  the  forward  rotor  Increases  with  airspeed. 


•  M 


Rotor  shaft  torque,  Mg  in  Figures  lib  and  lie  averages  about  160,000  inch- lb 
for  the  aft  rotor  and  about  60,000  inch- lb  for  the  forward  rotor.  Th*  rotor 
directions  are  as  shown  in  Figure  12.  The  torque  differential  of  100,000 
in- lb  duo  to  the  larger  aft  rotor  torque  acts  counter  clockwise  on  the  aircraft, 
and  Mist  be  balanced  by  approximately  equal  and  opposite  lateral  loads  at 
both  rotors.  The  necessary  load  magnitude  is  approximately  200  lb.  In  Fig¬ 
ures  lib  and  lie  it  is  seen  that  these  loads  do  occur  in  the  measured  Fy, 
averaging  about  200  lb  to  the  right  at  the  forward  rotor  and  to  the  left  at 
the  aft  rotor.  Roll  moment  Mg  is  also  related  to  this  effect  since  the 
forward  rotor  roll  moment  is  roll  right,  and  the  aft  rotor  roll  moment 
is  roll  left.  In  hover,  the  rotor  torques,  become  nearly  equal. 

■IhUd  Htnwmlc  Fixed  System  loads  -  With  a  three  bladad  helicopter,  third 
harmonic  vibration  levels  arc  usually  tha  dominant  portion  of  the  vibration 
environment.  Figure  13  presents  the  measured  third  harmonic  fixed  system 
loads.  The  second  and  fourth  harmonic  rotating  shaft  loads  from  which 
the  third  harmonic  in-plane  loads  result  are  presented  in  Figures  15'  and 


Figure  13a, 
average , 

Rotor 

Forward 

Aft 


1  airspeed  sweep, 

shows  that  shaft 

force  resultants 

Fg  lb 

Fy  lb 

Fz  lb 

Lateral 

Vertical 

800 

250 

200 

1300 

300 

400 

The  longitudinal  loads  are  easily  the  larger,  and  the  aft  rotor  loads  are 
larger  than  those  at  the  forward  rotor.  This  is  particularly  marked  at 
60  knots  where  the  aft  load' is  almost  three  times  the  forward  load.  Mo 
clear  general  trend  with  airspeed  is  evident.  Tha  vertical  force  peaks 
at  60  knots  on  the  aft  rotor,  however, the  forward  rotor  vettlcal  has  a 
minimal  in  .this,  speed  region  and  shows  a  rapid,  fairly  staep  rise  to 
103  knots . 


Figures  13b  and  13c,  rotor  speed  sweeps  at  40  knots  and  90  knots,  indicate 
an  extreamly  rapid  load  buildup  with  rpm  for  longitudinal  load  at  both 
rotors  and  for  vertical  load  at  the  aft  rotor.  The  sine  component  is 
largely  responsible  for  this  trend  at  40  knots,  while  both  the  sine  and 
cosine  contribute  at  90  knots  aft  rotor  longitudinal  load  was  nearly  twice 
as  large  as  the  corresponding  forward  rotor  load. 

Third  harmonic  vertical  loads  arise  from  direct  addition  of  third  harmonic 
blade  root  vertical  shears.  Third  harmonic  ln-plane  loads  result  from 
second  and  fourth  harmonic  bfada  root  horisontai  shears  as  shown  in 
Figure  14. 


ry 
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Floras  14.  THIRD  HARMONIC  W-PLAHE  LOADS 


Stsgad  isUauaJi  Jteaaai£.  EataUai.  Ssuai  teak  -  n$wc*»  is  and  u  pr«*nt 

tha  rotating  second  and  fourth  hansonlc  in-plans  loads.  Tht  second  harmonic 
rasultant  loads  show  a  tendency  toward  the  high  hover,  low  cruise  characteristics 
In  the  forward  rotor,  but  not  in  the  aft  rotor.  The  fourth  harmonic  airspeed 
sweep  produces  relatively  constant  resultants  within  nose  cases,  a  tendency 
to  increase  sharply  at  high  forward  speed. 

Numerical  combination  of  these  loads  Into  the  fixed  system  third  harmonic 
Is  illustrated  in  Figure  14  for  the  forward  rotor  at  20  knots.  The  second 
and  fourth  harmonic  rotating  loads  are  additive  in  forming  the  longitudinal 
fixed  system  third  harmonic  load,  subtractive  in  forming  the  lateral  fixed 
system  harmonic  load,  hence  longitudinal  loads  are  in  almost  all  instances 
larger  than  lateral  loads. 

It  is  of  interest  to  note  certain  similarities  among  the  second  harmonic 
airspeed  sweep  curves  in  Figure  15s.  The  Fx  sine  curve  of  the  foxwerd 
rotor  is  the  negative  of  the  Fy  cosine  curve;  the  Fy  eine  curve  is  the 
tame  as  the  Fx  cosine  curve  of  the  forward  rotor.  Tha  aft  rotor  shows 
similar  characteristics,  however  the  pattern  differs  in  sign  relation  because 
of  the  opposite  direction  of  shaft  rotation.  The  fourth  hsrmonlc  airspeed 
sweep  illustrates  its  own  pattern,  Figure  16a.  The  Fx  sine  curve  is  the 
negative  of  the  Fy  cosine  curve.  The  tarns  similarity  patterns  between 
longitudinal  and  lateral  second  and  fourth  harmonic  loads  are  also  evident 
in  the  rpm  sweeps. 


FIGURE  16a.  ROTATING  SYSTEM  FOURTH  HARMONIC  SHAFT  LOADS 
AIRSPEED  SHEEP  AT  258  RPM 


FIGURE  16c.  ROTATING  SYSTEM  FOURTH  HARMONIC  SHAFT  LOADS 
RIM  SHEEP  AX  90  KNOTS 


Conclusions 


An  extensive  rotor  hub  load  measurement  program  has  been  conducted t  giving 
a  complete  picture  of  all  the  steady  and  oscillatory  loads  acting  in  the 
rotor  shafts  of  a  tandem  helicopter.  Shaft  strain  gage  instrumentation 
was  developed,  calibrated  and  flown  which  had  average  accuracies  of  6% 
for  steady  load  and  4%  for  oscillatory  load. 

With  a  three  bladed  helicopter,  third  harmonic  vibration  levels  are  usually 
the  dominant  portion  of  the  vibration  environment,  so  that  the  third 
harmonic  'oads  measured  here  are  particularly  important.  Longitudinal 
shaft  loaas  were  the  largest,  averaging  800  lb  at  the  forward  rotor  and 
1300  lb  at  the  aft;  lateral  loads  averaged  250  and  500  lb  forward  and  aft 
respectively;  vertical  loads  averaged  200  and  400  lb  forward  and  aft 
respectively.  Roth  roll  and  pitch  moments  averaged  2,500  in.  lb  at  the 
forward  rotor  and  5,000  in.  lb  at  the  aft  rotor. 

In-plane  loads  exhibited  a  trend  with  airspeed,  low  in  hover,  high  at  20 
knots  transition,  low  at  80  knots  cruise,  and  then  a  rise  toward  the  top 
speed  at  105  knots.  Vertical  loads  differed;  the  forward  rotor  was  high 
at  hover  and  105  knots,  low  at  60  knots;  the  aft  rotor  was  low  at  hover 
and  105  knots,  and  peaks  at  60  knots.  Rotor  speed  sweeps  at  40  and  90 
knots  indicated  an  extremely  rapid  load  buildup  with  rpm  for  longitudinal 
load  at  both  rotors  and  vertical  load  at  the  aft  rotor.  Lateral  loads 
and  the  moments  were  not  strongly  affected  by  rpm. 

Third  harmonic  vertical  loads  arise  from  direct  addition  of  the  vertical 
components  of  third  harmonic  blade  root  normal  flap  shears.  Third  harmonic 
in-plane  loads  come  from  second  and  fourth  harmonic  blade  root  lag  shears 
and  from  horisontal  components  of  second  and  fourth  harmonic  blade  root 
normal  flap  shears.  The  measured  second  and  fourth  harmonic  resultant 
rotating  loads  were  the  same  in  the  longitudinal  and  lateral  directions; 
in  the  rotating  system  these  directions  refer  to  the  master  hub  spline 
aslmuth  and  its  perpendicular.  The  equality  of  these  rotating  longitudinal 
and  lateral  loads  is  a  proof  that  the  causative  blade  root  shears  are 
harmonically  repetitive  at  each  of  the  three  blades.  Second  harmonic 
rotating  loads  averaged  500  lb  at  the  forward  rotor;  900  lb  at  the  eft 
rotor.  They  displayed  a  tendency  toward  the  high  hover,  low  cruise, 
force  levels  which  are  characteristic  of  the  fixed  third  harmonic  airspeed 
sweep,  and  showed  a  definite  Increase  with  rotor  speed.  Fourth  harmonic 
rotating  loads  averaged  250  lb  at  the  forward  rotor  and  500  lb  at  the 
aft  rotor.  There  were  moderate  changes  with  airspeed,  but  large  increases 
with  rotor  speed.  In  converting  to  the  fixed  third  harmonic  system, 
second  and  fourth  add  to  produce  large  longitudinal  loads  and  subtract 
to  produce  esmll  lateral  loads. 
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INTRODUCTION 


During  the  development  of  the  HH-43B  Huskle  helicopter  rotor  system, 
which  uses  a  servof lap-controlled,  torsionally  flexible  rotor  blade  in 
a  teetering  hub  configuration,  the  opportunity  arose  to  obtain  a 
United  amount  of  flight  test  data  in  areas  related  to  rotor  unsteady 
aerodynamics  and  dynamic  loads.  The  data  obtained  allowed  some  insight 
into  the  significance  of  unsteady  aerodynamics  on  forward  flight  rotor 
stability,  of  the  influence  of  rotor  stability  margins  on  rotor  loading 
in  normal  operating  regimes,  and  of  possible  ground  test  methods  for 
substantiating  rotor  flutter  stability  margins.  This  information  is 
presented  here  in  summary  form. 


BACKGROUND 

During  the  conduct  of  the  HH-43B  structural  demonstration  program  in 
1959,  flutter  was  encountered  during  the  buildup  to  high-speed  dive 
test  conditions.  Two  distinct  flutter  modes  were  found,  one  involving 
servoflap  bending  and  rotational  deflection  (against  control  system 
stiffness)  and  the  other  Involving  blade  bending  and  servoflap  rota¬ 
tional  deflections.  These  flutter  modes  were  found  to  be  critically 
dependent  on  rotor  lift  coefficient  and  dynamic  pressure,  occurring 
only  in  high  speed  partial  power  descents. 

Subsequent  whirl  testing  at  speeds  well  beyond  normal  overspeed  whirl 
test  requirements  showed  that  it  was  possible  to  duplicate  the  in¬ 
flight  occurrences  during  static  whirling.  This  enabled  standard 
flutter  elimination  techniques  (mass  balance,  stiffness,  and  hinge 
moment  changes)  to  be  used  to  raise  the  flutter  speed  well  above  the 
operating  range,  and  to  substantiate  this  change  prior  to  flight. 
Subsequent  limit  speed  flights  in  the  critical  power  ranges  confirmed 
the  satisfactory  nature  of  the  changes  made,  showing  no  re-occurrence 
of  flutter  in  the  entire  flight  envelope. 


DISCUSSION 


Figure  i  shove  a  typical  variation  of  vibratory  servoflap  load  with 
collective  control  position  at  constant  airspeed.  The  data  vas  obtained 
by  establishing  an  airspeed  in  level  flight  and  then  slowly  lowering 
collective  while  increasing  rate  of  descent  so  as  to  maintain  airspeed. 
Collective  variation  vas  sufficiently  slow  as  to  sake  each  data  point 
shown  a  quasi-steady  state  point.  When  the  collective  control  position 
reached  48  percent  of  the  level  flight  requirement,  a  very  sudden  in¬ 
crease  in  flap  vibratory  load  vas  noted.  Detailed  exaainatlon  of  the 
oscillograph  data  shows  this  increase  in  load  to  have  resulted  from  the 
sudden  appearance  of  a  new  non- harmonic  frequency  in  the  loading.  This 
frequency  appeared  each  time  the  blade  reached  the  maximum  dynamic 
pressure  (advancing  blade)  region,  then  damped  out  rapidly  as  the  blade 
total  velocity  vas  reduced.  Recovery  vas  accomplished  by  a  further  re¬ 
duction  in  power,  which  again  returned  flap  loading  to  normal  levels. 

Whirl  testing  shoved  this  particular  flutter  mode  to  be  a  coupling  of 
rotor  blade  third  symmetrical  bending  deflection  with  servoflap  rota¬ 
tional  deflection  (against  control  system  stiffness).  Flutter  only 
occurred  vhen  rotor  leading  vas  very  close  to  zero  lift  (approximately 
+.10%  of  one  "g"  thrust).  Flutter  did  not  occur  at  integer  multiples  of 
rotor  speed  nor  at  any  fixed  multiple  of  rotor  speed. 

Analysis  of  blade  loading  in  flight  also  shoved  that  the  lift  coeffi¬ 
cient  on  the  advancing  blade  at  the  flutter  point  vas  very  close  to 
zero.  Flutter  frequency  vas  also  equal  to  third  symmetrical  blade  bend¬ 
ing  frequency.  Apparently,  then,  a  wake  coupling  existing  in  hover 
operation  and  shoving  a  strong  variation  with  lift  coefficient  persisted 
with  similar  lift  coefficient  effects  out  to  advance  ratios  of  at  least 
0.29. 

Since  standard  flutter  analysis  methods  shoved  extremely  poor  absolute 
correlation  with  the  test  results,  the  trends  indicated  to  be  signifi¬ 
cant  by  analysis  were  explored  on  the  whirl  test  facility.  Use  of 
servoflap  mass  balance  vas  shown  to  provide  the  most  expeditious  means 
of  eliminating  the  flutter  mode.  This  approach  vas  substantiated  on 
the  whirl  rig  and  then  in  flight,  yielding  the  data  shown  in  Figure  2. 
This  plot  shows  the  variation  of  maximum  vibratory  flap  load  with  air¬ 
speed  before  and  after  the  incorporation  of  the  servoflap  mass  balance. 
Belov  the  flutter  speed  point  no  influence  of  the  flutter  mode  on 
vibratory  flap  loads  can  be  detected,  both  the  the  original  and  the 
revised  configuration  shoving  identical  loads.  Similar  trends  wero 
noted  in  blade  torsional  loadings. 

It  might  be  of  Interest  to  note  here  that  flutter  loadings  vere  consid¬ 
erably  higher  during  static  whirl  testing  than  in  flight.  This  appeared 
to  be  due  to  the  intermittent  entry  and  departure  from  the  flutter  zone 
in  flight  as  compared  to  the  continuous  penetration  of  the  boundary 
during  static  whirl.  The  intermittent  operation  appearied  to  substan¬ 
tially  limit  flutter  load  buildup. 


The  Incorporation  of  servoflap  nans  balance  shifted  the  blade  chordwise 
center  of  gravity  aft.  An  predicted  in  Reference  Xj  thin  change  In 
chordwise  c.g.  resulted  in  an  Increase  in  high  speed  rotor  vibratory 
shear  and  bending  moment.  it  vas  therefore  necessary  to  add  additional 
blade  mass  balance  to  compensate  for  the  change  produced  by  the  flap 
counterweight .  Figure  3  shows  the  effects  of  this  change  on  blade  root 
flatwise  vibratory  loading,  and  compares  the  results  with  the  load 
variation  predicted  in  Reference  1.  Somewhat  larger  flight  test  varia¬ 
tions  are  seen  than  were  calculated.  It  should  be  pointed  out  that 
this  change  in  blade  root  loading  increased  rapidly  with  airspeed, 
indicating  the  Increasing  Importance  of  blade  torsional  deflections  on 
flatwise  loading  as  airspeeds  are  Increased.  Using  cg-ac  offset  as  a 
means  of  vibratory  load  reduction  at  high  speed,  as  suggested  by  Reference 
1,  appears  to  be  a  powerful  tool;  a  1/2%  change  in  cg-ac  offset  produced 
an  18%  reduction  in  flatwise  vibratory  loading. 

It  has  already  been  pointed  out  that  static  whirl  testing  was  successful 
in  reproducing  the  characteristic  frequency  and  lift  coefficient  sensiti¬ 
vity  of  the  in-flight  flutter  modes.  It  is  interesting  to  also  note, 
as  shown  in  Figure  4,  that  whirl  test  flutter  speeds  and  flight  test 
flutter  speed  correlated  well  on  a  servoflap  dynamic  pressure  basis. 

While  this  result  appears  to  be  somewhat  different  than  results  obtained 
on  model  rotors  in  the  wind  tunnel,  and  hence  may  only  be  applicable  to 
rotor  systems  with  localized  control  devices,  it  has  provided  a  convenient 
means  of  providing  successful  ground  substantiation  of  flutter  integrity 
in  subsequent  servoflap  rotor  configurations.  Flight  test  investigation 
on  these  rotor  systems  has  continued  to  prove  that  freedom  from  flutter 
on  the  whirl  rig  means  freedom  from  flutter  in  flight  up  to  similar 
servoflap  dynamic  pressure.  Combinations  of  overspeed  and  propeller 
slipstream  velocity  were  used  in  achieving  desired  test  "q".  Additional 
testing  to  define  the  usefulness  of  this  method  for  general  flutter 
substantiation  would  be  desirable. 


SUMMARY 

In  summary  the  data  shown  here  indicates  that: 

1.  Reduced  damping  at  zero  lift  coefficient  (wake  effects)  apparently 
persists  at  substantial  advance  ratios  (up  to  0.3) 

2.  Torsional  coupling  and  relation  between  operating  chordwise  c.g. 
and  flutter  boundary  can  have  a  significant  Influence  on  blade 
vibratory  loading. 

3.  Constant  dynamic  pressure  testing  offers  some  promise  of  providing 
a  means  for  static  flutter  margin  substantiation. 
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1.  Helicopter  Blade  Vibration  and  Flutter,  R.  K.  Miller  and  C.  W.  Bills, 
Journal  of  the  American  Helicopter  Society,  July  1936 
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If  mt  have  to  evaluate  the  aerodynamic  loads  ever  a  helicopter  rotor 
„  <3,  w*  must  take  in  account  the  effect*  produced  by  wake  vortices  issuing 
jm  '.tv a  blades, 

Our  purpose  is  to  define  these  wake  vortices  in  permanent  state  by  using 
.its  at  tab  to  performed  in  the  wind  tunnnl.  This  will  be  done  in  three 

.  .1  i 

'-at  we  shall  show  the  location  of  helix  vortices, 

} 

’na->  .<&  shall  apply  these  results  xo  e  sample  calculation  in  havering  'light, 
i-'ally  we  shall  give  some  ,-ast xts  for  the  transient  dynamic  case. 

T-.a  most  extensive  part  of  the  results  presented  in  this  paper  comes  from 
f  ’f *perimentsl  and  theoretical  study  of  local  induced  velocities  over  a  rotor 
«k  •■xc  analytical  evaluation  of  tha  primary  loada  acting  on  helicopter  rotor 
•  ;ssB,  This  work  was  carried  ouc  by  “GIRAVIONS  DOhAND '-to."  under  the  sponsor- 
u-  f  trie  EUROPEAN  RESEARCH  OFFICE  US.  DEPARTMENT  OF  THE  ARMY. 

i  tane  hare  tha  occasion  to  thank  the  U.S.A.  TRECOM  who  made  it  possible 
.  ”...xa  work  to  be  accomplished  and  to  be  presented  nere  to  day. 

Ail  the  wind  tunnel  teats  worn  performed  wind  tunnel  of  the  Institut 
•'f-  an  10 uas  de3  Fuidas  at  Marseilla. 

■option  _qf  free  vortices 

-  Slide  1  shows  the  case  where  the  circulation  1"*  is  constant  along 
tna  span  but  varies  with  time.  The  straight  bides  of  the  classical  hor^e 
shoe  vortices  are  now  helix  like  curved  lines.  Tha  radial  vortices 
appear  as  tha  rate  of  circulation  P  changes  with  time. 

Slide  ?  shows  the  theoretical  case  where  the  circulation  P  varies 
along  thr  «,  an.  The  rapid  rate  of  change  of  P  et  the  tip  of  the 
tlada,  between  points  2  and  3,  produces  a  well  located  vortices  which 
-an  easily  be  assn  in  wind  tunnel.  The  slow  rate  of  change  of  circu- 
iation  P  at  tha  root  of  the  blade  between  1  and  2  givse  a  more 
diffusnr  vortice. 


S-  ide  j  effyrs  a  thaoretical  ej'planation 


Slide  4  shows  the  tip  vortiee.  The  helix  is  strongly  deformed  at  the 
rear  section  under  the  disc* 

Slides  5  and  6  show  the  root  helix  vortiee  which  is  moving  in  a  tur¬ 
bulent  region. 

Slides  7  and  8  show  the  tip  helix  vortices.  1  have  drawn  the  contour 
line  of  the  root  vortiee.  The  interference  between  the  two  vortices 
is  the  cause  of  the  strong  deformation  of  the  helix  at  the  rear 
sector  of  the  disc. 

For  this  reason  we  can  not  aasune  a  regular  helix  form  and  we  are 
obliged  to  consider  the  real  helix. 


1 .2  -  Tip  vortices 

Slide  9*  For  an  instantaneous  azimuthal  position  1  of  the 

blade  we  have  measured  the  vertical  distance  »AtA%  between 

two  successive  vortices*  In  this  way  we  know  the  correct  distance 
of  the  vortics  in  the  neibsurhood  of  the  instantaneous  blade  position* 
This  measursment  was  made  on  corresponding  pairs  of  photos* 


Ws  have  achieved  this  in  two  different  ways. 


In  the  first  method  smoke  was  emitted  from  the  tip  of  the  blade. 
We  have  already  seen  a  pair  of  these  photos* 


In  the  second  method  (slide  10  and  11)  the  smoke  emission  was  sxternsl 
and  was  directed  so  that  the  tip  of  the  blade  cut  the  smoke  at  different 
given  azimuthal  angles. 


-  For  every  "flight"  we  repeated  the  experiment  for  24  different  azi¬ 
muthal  values. 


-  We  may  write  that  the  vertical  distance  is  product  of 

inflow  velocity  time  the  time  interval  which  separate  two 
successive  vortices* 


(1,1) 


b 

T,  M 

12 

b=.  7 


velocity  is  the  inflow  velocity  induced  by  the  wakes. 

This  velocity  is  not  the  induced  velocity  to  be  introduced  in  the 
thrust  formula  but  it  helps  us  to  locate  correctly  the  wake  vortiee 
in  the  naighJbourhood  of  the  instantaneous  blade  position. 


We  have  compared  this  inflow  Velocity  to  the  mean  induces  FROUDC 
velocity  Yi*  given  in  the  text  book  i 

(1.2) 
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Slide  12  -  The  typical  variation  of  non  dimanaionn&l  inflow  velocity 
with  azimuth  is  shown  in  elide  12. 

Slide  13  ahowa  the  same  variations  around  the  disc. 

1.3  -  We  have  studied  the  different  influence*  of  t 

-  collective  pitch 

-  tip  speed  ratio 

-  angle  of  attack  of  the  rotor 

(  oC  is  considered  negative  for  nose  up) 

-  cyclic  pitch  (pitch  and  roll) 

-  solidity 

-  Lock* a  constant* 

Different  curves  are  given  in  the  tsxt*  The  moat  important  influences 
come  from  the  tip  spaed  ratio,  angle  of  attack  of  the  rotor  and  cyclic 
pitch.  The  collective  pitch  and  solidity  influence  mainly  the  mean 
value  . 

1.4  -  We  era  now  able  to  give  coordinates  of  an  slamsnt  of  a  helix  vortice. 

Slide  14.  Let  us  consider  at  time  t  a  vortice  element  wich  had  left 
the  blade  ^  seconds  ego  thet  is  to  say  at  t  - 

-  Far  an  element  of  tip  vortice  we  propose  * 

'  02r  =  -Vc^y.ft .-Ta)  ■+ 
ii,3)  y  r  =  Rs<tia..(t-z,) 

zT=  ~Rf>  It-r.) 

t,  =  jtt(t  -zrj 

This  aene  that  after  leaving  the  blade  the  vortice  element  has  been 
taken  by  tits  resultant  wind  due  to  flow  velocity  plus 
inflow  velocity  is  shown  in  tha  curves  givsn  in  the  text. 

-  For  an  slamsnt  of  helix  vortics  we  propose  the  following  formulas t 

r  oCq  =  £^)  4  r0  cosJ2 .  ( t-T-, } 

,  u  r0  Sinfl.lt-  V.) 

<M>  ,  . 

**  =  -r0ft  Lt-zJ  (t-rj 

After  a  period,  T#  is  situated  approximately  at  the  center  of  gravity 
of  the  poaitiva  derivatives  -ML  erf  the  circulation. 

Dr 
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US  »  Slid*  15  -  The  stxsnghs  P  of  tho  tip  vertices  »re  equal 
and  opposite  < 

(US)  «[-|.  Ct.  C  .  (Jlr  ♦  V*m  ^jJ  ^ 

Tha  value  can  be  dovolopped  aa  a  FOURIER  aariaa  and  the  firet 

tax*  can  be  estimated  by  tha  formula  of  MEIJER  DREES. 

1.6  -  Wa  may  now  calculate  the  location  of  waka  vorticaa. 

-mt 

2.1  -  How  may  we  apply  the  previoua  raaulta  to  calculate  aerodynamic  loada 

along  the  blada. 

Slide  16  -  We  shall  treat  tho  wsks  of  the  blade  in  two  zones.  Tha 
first  zone  is  in  the  plane  of  the  disc}  it  is  defined  by  two  azimu¬ 
thal  positions  t 

instantaneous  ^  *  SL  t 

and  <NJ2/fc-X 

T  _  2K 
12 

This  zone  can  be  treated  by  methods  similar  to  those  used  for 
fixed  wings, 

-  The  second  zone  of  the  wake  comprises  two  helix  vortices  and  possibly 
radial  vortices.  We  estimate  that  if  the  disc  loading  is  great 
enough  we  can  teke  into  account  only  helix  vortices  with  a  constant 
T\  circulation  strength. 

2.2  -  As  a  numerical  example  we  took  the  two  blade  rotor  of  Mr  John 

R.RABBOTT  Jr.  (NACA  T.N.  3686  July  1956,  and  we  shall  compare  the 
theoretical  results  to  the  axpsrimentol  results  found  by  Mr  RABBOTT. 

Slide  18.  The  local  circulation  P  (  r )  at  a  station  r  along  the 
A  blade  is  given  by  the  lift  line  equation  t 


(2,2) 


The  induced  velocities  are  produced  by  the  wake  of  blade  A  and 
blade  2.  T^e  velocity  produced  by  the  weke  of  blade  A  is  divided 
into  two  zones  (1  and  2)  the  new  equation  can  ba  written  ae  follows  t 
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The  unknown  quantities  in  t 

PA  (r)  8  function  of  T 


I 

I 

I 

I 

I 

I 


(2,4) 
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maxim*  value  of  circulation 


inflow  velocity 

The  inflow  velocity  value  is  calculated  by  the  elementary  PROUDE  strip 
f annul e  et  r  «  0,5. 


-  To  resolve  the  equation  we  used  an  iterative  method.  We  estimate  the 
value  of  r*m  and  we  calculate  the  P*  ( t~)  function.  Then  we  draw 
the  curve  P*  (  r)  to  obtain  a  new  maxima  value;  the  process  is 
repeated  and  in  the  case  presently  considered  we  found  we  required 
four  iterations  to  obtain  a  good  solution  (Slide  19). 

Tha  comparison  between  the  r"^opoead  method  (Slides  20  and  slide  21)  and 
the  experimental  results  is  valuabia.  However  it  will  bs  nacssaary  to 
examine  other  exaaplf>c  to  verify  tha  theory. 
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3  -  are  k,suut§  rep  the  vmm.  stats  amtiwwcs 

3.1  -  These  a  few  underlying  results  of  transient  stats  aerodynamic  will 

help  us  to  propose  a  treatsmeqt  which  take  into  account  the  effect 
of  radial  vortices  in  permanent  state. 

3.2  .  If  we  change  the  collective  pitch  from  0to  to  (Slides  22,  23, 

)  the  response  in  thrust  of  the  rotor  can  present  very  different 
aspects  because  there  are  many  phenomena  which  may  interfere. 

-  Primary  factors  seems  to  be  4 

a)  the  rata  of  change  of  collective  pitch 

b)  flapping  movement 

c)  disc  loading. 

We  shall  now  consider  a  few  results  experimentally  observed. 

3.3  «  The  two  blade  rotor  had  heavy  blades  (the  LOCK  constant  was  approxima- 

tively  0.64).  This  enabled  the  Mika  aerodynamics  to  be  observed.  The 
flapping  movement  is  small  but  poorly  damped  (Slides  24  end  25). 

The  circulation  P  around  the  blade  is  1 
(3,1)  P  -  K  (  9c  -  -  f.  ij  j 

&£  is  the  induced  angle  by  the  first  rone  of  the  weka. 
ijk  is  the  induced  angl*  by  the  second  rones  of  the  wska. 
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K*  the  circulation  P  changes  with  time  tha  tip  helix  vaartica  is  fad 
with  new  vertices  .  Thaaa  vorticaa  can  ba  saan  in  tha  wind  tunnel 
by  aaanc  of  aaoka  trocars,  Tha.ir  number  inexaasaa  with  tha  rata 
of  change  of  collsctive  pitch  (Photos). 

Thaaa  vorticaa  produce  c  double  effect  s 

t  -  They  pass  under  tha  blade  and  aay  induce  tha  LOtWY'S  (Van  VOOREN) 
affect. 


2  -  They  deform  the  wake  and  change  its  intensity. 


We  have  effected  i..  hovering  "flight"  a  alow  rata  of  change  of 
collective  pitch,  the  distance  between  two  auocneeiva  vorticaa  was 
greater  than  1  —1,5  chard  lengthe.  Within  bo under iee  of  this  expe¬ 
riment  remarks. 


1  -  If  we  change  collective  pitch  from  Gt#  to  the  final  perma¬ 
nent  state  is  established  after  the  vortice  has  travailed 

a  vertical  distance  equal  to  approximative^  ana  blade  chord. 


2  -  The  thrust  variation  with  time  is  not  exactly  tha  same  if  we 

diminish  or  increase  the  pitch. 

3  -  If  the  change  of  pitch  is  sinusoidal  having  a  frequency  from 

Q.5  to  3.5  cps,  there  is  phenelsg  between  the  pitch  and  thrust 
variation  of  about  cwte  period  ( T  = 

The  thrust  is  proper tionnal  to  2 


(3,2) 


e,  -  At.  ee  =  fft  -  At) f(t>  -  At  f(t) 

i 

At  4  T 


,4  -  Tha  influence  of  flapping  may  have  two  aspects  « 


1  -  Tha  flapping  velocity  changes  the  circulation 
blade  i 


of  the 


(3.3)  Pr  k[  8c-  -  Ztr  -±  j>  J 

2  -  The  flapping  modifies  the  distance  from  the  blade  to  the  wake 
under  tha  disc.  This  distance  measured  in  half  chords  is  t 

(3.4)  JiAL  -  JL .  Jl 

o i  C  V  XST 


d-  angle  of  incidence  of  the  blade 
O'  local  solidity, 

*hi»  distance  can  be  quit  great  in  two  cases  i 


a)  -  light  blades  excited  at  frequencies  between  0  and  JCI  » 

b)  -  all  bin das  at  the  resonance  frequency  SX 

The  two  previous  phenomena  are  bound  together.  We  can  observe  thee  on 
the  slide  {  2%  ),  As  the  pitch  has  reached  its  final  value  there 
remains  soma  ti»«  after  only  the  movment  of  the  blade  it  self. 

Theoretically  the  thrust  variation  is  proper tionnel  to  t  (Slid#  36) 

±i 

a  I 

(Slide  22}  In  fact  the  thrust  oscillates  with  a  period  30  %  greater 
then  T  -  period  of  revolution. 

The  deeping  of  this  oscillation  is  lower  by  about  one  half  than 
theoretical  one. 

3.5  -  Ccneluaion 

We  estimate  that  in  usual  permanent  flight  the  effect  of  the  variatior 
of  circulation  3T.  which  is  attached  to  radial  and  helix  vortices 
will  be  a  r  corrective  only.  We  suggest  the  following  form  t 

n(fc-/u)  z  ri t)  -ju  r(t> 

4! 


muisistt 

1  -  We  tried  to  define  the  woks  vortices  be  using  extensive  experimental 

results  in  non  dimensional  form. 

2  -  The  strong  interference  between  the  two  hall*  deforms  the  regular  curves 

and  we  era  obliged  to  consider  the  reel  helix.  This  phenomenum  is  cer¬ 
tainly  influenced  by  the  neighbourhood  of  the  fuselage  for  low  *  kt  ", 
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ROTOR  AIRLOADS 

R.  Piziali,  H.  Daughaday,  F.  DuWaldt 
Cornell  Aeronautical  Laboratory,  Inc. 


ROTOR  AIRLOADS* 

Raymond  Piziali,  Hamilton  Daughaday,  Frank  DuWaldt 
Cornell  Aeronautical  Laboratory,  Inc. 

INTRODUCTION 

The  problem  of  adequately  predicting  the  aeroelastic  response  of  rotating 
wings  has  been  one  of  the  most  formidable  problems  to  confront  the  aircraft 
industry.  Industry  has  had  to  rely  on  expensive  semiempirical  methods  and 
previous  experience  to  design  and  develop  these  aircraft  components.  Such 
procedures  do  not  provide  the  designer  with  a  clear-cut  rational  approach  to 
his  problem  and  they  provide  little  hope  of  optimizing  the  blade  design  to  any 
given  set  of  specifications. 

The  aeroelastic  problem  can  be  thought  of  as  consisting  of  two  closely 
related  parts:  the  aerodynamic  loading,  and  the  dynamic  response  of  the 
system.  They  are,  of  course,  interdependent  and  related  such  that  the  aero¬ 
dynamic  loading  excites  the  response  of  the  structure  and,  at  the  same  time, 
depends  on  it.  In  general,  the  dynamic  response  of  the  system  will  consist 
of  fuselage  motions  and  blade  motions  —  i.  e. ,  rigid  pitching,  flapping,  lagging, 
and  all  the  coupled  bending  (both  flapwise  and  chordwise)  and  torsion  modes  of 
the  blades  — while  the  aerodynamic  loading  will  consist  of  the  total  reaction 
of  the  air  on  the  blade  sections.  This  aerodynamic  loading  will  include  the 
unsteady  lift,  induced  drag,  pitching  moments,  and  profile  drag  which  are 
appropriate  for  attached  and  separated  flow  conditions  with  the  blade  sections 
operating  in  normal  and  reversed  flow. 

The  vortical  wake  of  a  rotating  wing  is  quite  complex;  however,  for 
steady-state  translation  of  the  rotor  and  as  viewed  from  the  wing  in  the 
rotating  coordinate  system,  it  is  periodic.  The  wake-induced  velocities  at 
the  wing  will,  therefore,  also  be  periodic  and  these  together  with  the  periodic 
variation  ir.  the  local  tangential  velocity  (dynamic  pressure)  will  cause  the 
aerodynamic  loading  to  be  periodic.  Thus,  due  to  the  time- varying  airloading. 


*  This  work  was  accomplished  under  the  sponsorship  of  the  U.  S.  Army  TRECOM. 


the  vortical  wake  will  consist  of  continuous  distributions  of  shed  and  trailing 
vorticity  which,  at  any  instant,  will  be  generally  concentrated  on  a  skewed 
and  distorted  helical  surface;  these  distortions  of  the  vortical  wake  from  a 
skewed  helical  surface  are,  of  course,  the  result  of  the  time-varying  three- 
dimensional  induced  velocity  field  of  the  wing  and  wake. 

The  above  discussion  points  out  the  general  complexity  of  the  over -all 
aeroelastic  problem  for  rotating  wings  and  the  interdependency  of  its  various 
aspects.  Any  practical  method  of  predicting  this  aeroelastic  response  will 
certainly  make  use  of  many  simplifying  assumptions.  Some  of  these 
assumptions  will  be  generally  applicable;  however,  there  will  be  others 
whose  validity  will  depend  on  the  particular  wing  and  operating  condition 
being  analyzed.  For  example,  in  general,  it  will  not  be  possible  to  neglect 
the  aerodynamic  coupling  of  the  blade  modes  but  for  a  specific  problem  this 
may  be  a  valid  assumption.  Also,  it  will  probably  not  be  possible  to  compute 
the  aerodynamic  loading  and  the  structural  response  independently;  they  will 
have  to  be  computed  simultaneously  by  satisfying  the  equations  of  motion  or, 
possibly,  by  iterating  between  separate  loading  and  response  computations. 

Any  method  of  solution  of  the  rotating  wing  aeroelastic  problem  must 
adequately  predict  the  wake -induced  velocities  at  the  wing.  However,  it  is 
very  difficult  to  represent  mathematically  (or  computationally)  the  wake  effects; 
this,  in  the  past,  has  been  the  prime  obstruction  to  an  adequate  solution.  The 
attempts  of  earlier  investigators  to  solve  this  problem  necessarily  were  based 
on  relatively  drastic  simplifications  of  the  wake  in  order  to  make  the  problem 
computationally  feasible.  Modern  high-speed  digital  computing  machines  have 
now  made  it  physically  possible  to  account  for  much  more  of  the  realistic 
wake  detail.  It  is,  therefore,  now  possible  to  formulate  a  representation  of 
the  wake  which  will  enable  accurate  computation  of  the  wake -induced  velocities. 

This  paper  reviews  our  initial  approach  to  the  rotating  wing  aeroelastic 
problem  and  some  of  the  result*  of  this  effort.  A  detailed  discussion  of  this 
initial  effort  and  the  results  are  presented  in  References  1,  2,  and  3.  This 
is  followed  by  a  discussion  of  some  of  the  remaining  problems  and  the  line  of 
attack  currently  being  followed  in  their  solution. 
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THE  INITIAL  APPROACH  TO  THE  PROBLEM 


For  the  initial  effort  on  this  problem,  it  was  felt  that  it  would  be 
advantageous  to  concentrate  on  only  the  aerodynamic  loading  aspect  of  it 
(i.  e.  ,  to  bypass,  for  the  present,  prediction  of  the  structural  response). 

This  was  facilitated  by  using  measured  dynamic  response  as  a  known  input 
with  the  exception  of  the  first  harmonic  flapping  which  was  approxima^y 
forced  to  equilibrium  by  an  iteration  with  the  flapping  equation  of  motion. 
Furthermore,  an  attempt  was  made  to  simplify  the  problem  as  much  as 
possible  for  the  initial  effort  by  eliminating  all  but  the  most  essential  features 
and  effects. 

Each  blade  of  the  rotor  was  represented  by  a  segmented  lifting-line 
(bound  vortex)  located  along  the  steady  deflected  position  of  the  quarter -chord. 
The  number  of  segments,  n  ,  and  the  length  of  each  segment  are  arbitrary; 
they  are  each  straight  and  of  constant  vortex  strength.  The  lifting  line  is 
considered  to  advance  in  a  stepwise  manner  through  "  N  "  equally  spaced 
azimuth  positions. 

In  the  wake,  the  continuous  distributions  of  shed  and  trailing  vorticity 
of  each  blade  are  represented  by  a  mesh  of  segmented  vortex  filaments;  each 
segment  is  straight  and  of  constant  vortex  strength.  The  segmented  trailing 
vortex  filaments  emanate  from  each  of  the  end  points  of  the  lifting-lir.e 
segments.  The  segmented  shed  vortex  filaments  intersect  the  trailing 
filaments  in  a  manner  such  that  the  end  points  of  both  are  coincident  (Figure  1). 

The  strengths,  dP /dt  ,  of  the  shed  elements  are  equal  to  the  change  in 
strength  of  the  bound  vortex  segments  between  successive  azimuth  stations 
and  are  deposited  in  the  flow  at  each  azimuth  station  of  the  bound  vortex.  The 
strengths,  dP/dr  ,  of  the  trailing  vortex  elements  are  equal  to  the  differences 
in  strengths  of  adjacent  bound  vortex  segments  and  are  deposited  in  the  flow 
in  a  manner  such  that  they  connect  the  bound  vortex  end  points  to  the  shed 
vortex  end  points.  The  displacement  time  history  of  the  wake  elements  is  an 
input  in  this  method;  thus,  any  physically  realistic  distortion  of  the  wake  can 
be  incorporated  into  the  computation  of  the  airloads. 


For  each  azimuth  position,  the  airloads  are  computed  at  the  mid-points 
of  the  lifting -line  segments.  Thus,  the  airloads  are  computed  at  {  n  )  (  N  ) 
points  of  the  rotor  disk.  These  points  R*  of  the  disk  are  labeled  as  shown 
in  the  example  of  Figure  2;  subscript,  ,  refers  to  position  in  the  disk. 

The  expression  for  the  lift  per  unit  span,  for  a  blade  section  at  any  point,  Pfc  , 
of  the  disk,  is  taken  to  be 

i,  r.  u> 

where 


The  quantity  in  parentheses  is  the  instantaneous  '‘effective"  angle -of -attack 
at  tli e  three -quarter -chord  point,  where  e^k  is  the  instantaneous  geometric 
angle -of -attack  with  respect  to  the  tip -path -plane,  hk  the  "effective"  plunging 
velocity  with  respect  to  the  tip -path -plane,  and  the  velocity  component 
normal  to  the  tip -path -plane  induced  by  all  the  wake  vorticity  and  the  bound 
vorticity  of  other  blades.  The  strengths,  Hr  ,  of  the  bound  vortex  elements 
are  considered  to  be  the  unknowns  for  the  method  of  solution  developed. 


The  velocity  component  in  a  given  direction  induced  at  a  point  by  an 
arbitrarily  oriented  straight  vortex  filament  segment  of  constant  strength  is 
given  by  the  Biot-Savart  law  as: 

f  *  f  Tf*  <5) 

where  is  the  constant  vortex  strength  of  the  filament,  and  the  coefficient 
■f  is  a  function  only  of  the  coordinates  of  the  point  where  the  velocity  is  being 
computed  and  the  coordinates  of  the  vortex  filament  end  points.  The  velocity, 

,  of  Equation  (2)  can  be  computed  by  summing  the  contributions  (given  by 
Equation  (3)  )  of  each  individual  vortex  element  of  the  blades  and  wake  as 
indicated  by  Equation  (4). 


U  +£  fr 

a*i*. 


(4) 


However,  each  wake  vortex  element  strength  of  Equation  (4)  is  just  the 
difference  in  strength  of  two  of  the  unknown  bound  vortex  segment  strengths, 
fj  (i.  e. ,  ^  »  where  and  i>  are  two  values  of  the  disk 

position  subscript,  n  ).  These  ^  relationships  are  depicted  in  Figure  1 
relative  to  a  point  P*  of  the  rotor  disk.  If  these  differences  are  substituted 
for  the  f  ,  Equation  (4)  can  be  expanded  and  the  unknown  grouped  and 
factored  out.  Thus,  the  urh  can  be  expressed  as 

ntf 

h  *z  <51  a  <5> 


In  Equation  (5),  the  are  the  unknown  bound  vortex  strengths  and  each  ffisi 
iB  the  sum  of  all  the  coefficients  f  of  a  common  Cc  in  Equation  (4).  By 
substituting  the  Uf^  given  by  Equation  (5)  in  Equation  (2)  and  letting  1 A  = 
—/l#, )  »  the  following  expression  for  /*  is  obtained: 


n» 


i.+  a.t>.Z  n 


t 


(6) 


A  set  of  linear  simultaneous  nonhomogeneous  equations  with  constant  coefficients 
is  obtained  when  Equation  (6)  is  written  for  each  point,  ,  of  the  rotor  disk 
at  which  the  airload  is  to  be  computed.  Thus,  the  airloads  problem  is  reduced 
to  obtaining  the  solution  to  this  set  of  algebraic  equations  for  which  the  unknowns 
are  the  (  n  )  (  N  )  bound  vortex  strengths  f"J,  . 


The  effects  of  stall  were  approximated  by  effectively  limiting  the  lift 
coefficient  to  its  value  at  the  stall  angle  for  effective  angles  of  attack  greater 
than  the  stall  angle  of  the  blade  section.  It  is  possible  within  the  framework 
of  this  initial  method  to  use  more  representative  lift  coefficient  variations 
with  angle  of  attack  than  the  one  chosen  here.  It  should  be  noted,  however, 
that  incorporating  stall  in  this  manner  assumes  that  the  stalled  section  lift  is 
all  circulatory  and  there  are,  therefore,  shed  and  trailing  vortex  elements  in 
the  wake  with  strengths  determined  by  the  stalled  values  of  the  bound  circulation. 
In  reality,  the  predominant  part  of  the  airload  on  the  stalled  section  may  be 
noncirculatory  and  due  to  the  separated  flow  (i,  e. ,  analogous  to  the  pressure 
drag  on  a  flat  plate  at  large  angles  of  attack).  Some  indication  of  the  magnitude 


5 


of  the  separated  flow  forces  experienced  may  be  obtained  by  using  the  measured 
airloads  and  dynamic  pressure  to  compute  the  corresponding  normal  force 
coefficients;  they  are  found  to  be  extremely  large  for  portions  of  the  rotor 
disk  where  the  blades  are  evidently  stalled. 

Reversed  flow  effects  are  account  -  for  in  the  computation  approximately 
in  that  the  sign  of  the  lift  force  and  the  sense  of  the  bound,  shed  and  trailing 
vorticity  are  always  consistent  with  the  local  tangential  velocity,  V  ,  and  the 
"effective"  relative  angle -of -attack  of  the  section  as  it  approaches  and  passes 
through  the  region  of  reversed  flow.  However,  the  wake  remains  attached  to 
the  bound  vortex  of  the  blade  {although  its  position  and  orientation  corresponds 
to  the  local  inflow  velocity)  and  the  bound  vortex  remains  at  u  .  quarter -chord 
while  the  blade  approaches  and  passes  through  the  region  of  rever.-cU  Lu  ..  In 
reality,  when  a  section  is  operating  unstalled  in  reversed  flow,  the  center  of 
pressure  may  be  approximately  at  the  three -quarter -chord;  thus,  the  bound 
vortex  should,  possibly,  be  placed  at  that  station  and  the  induced  velocity 
calculated  at  the  quarter -chord. 

Mach  number  and  Reynolds  number  effects  can  be  approximated  as  vari¬ 
ations  of  lift  curve  slope  of  the  local  blade  se-tions.  For  a  given  rotor  and 
operating  conditions,  the  Mach  number  and  Reynolds  number  distributions 
over  the  rotor  disk  are  fixed  (known);  thus,  these  effects  can  be  incorporated 
into  the  computations  by  introducing  the  proper  distribution  of  lift  curve  slope 
over  the  disk. 

The  sensitivity  of  the  computed  airloads  to  the  number  of  revolutions  of 
wake  retained  and  the  wake  distortion  was  investigated  computationally  for  the 
NASA  model  rotor  (Reference  4)  at  an  advance  ratio  of  0.  15. 

For  this  rotor  and  operating  condition,  three  revolutions  of  wake  accounted 
for  practically  all  of  wake-induced  velocities;  extending  the  wake  to  five  revo¬ 
lutions  only  resulted  in  a  maximum  change  of  two  percent  in  the  magnitude  of 
the  lowest  harmonics  of  the  airloads  (i.  e. .  the  steady,  first,  and  second 
harmonica).  This  was  expected  because  the  distant  wake  will  have  a  relatively 
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uniform  influence  over  the  disk  and,  thus,  affect  primarily  only  the  lowest 
harmonics  of  the  results.  The  general  conclusion  ie  that  the  first  few 
revolutions  of  wake  are  the  most  important  with  respeet  to  dynamic  airloads, 
and  the  amount  of  retained  wake  will  depend  on  the  rotor  and  operating  con¬ 
dition  being  investigated. 

Sensitivity  of  the  computed  airloads  to  distortions  of  the  wake  from  a 
regular  skewed  helical  surface  was  investigated  for  several  relatively 
arbitrary  distortions;  the  sensitivity  experienced  was  moderate  for  this  case. 
The  only  conclusion  to  be  drawn  from  these  results  is  that  the  wake  distortion 
could  possibly  be  important  for  certain  specific  rotors  and  operating  conditions 
and  that  further  investigation  should  be  undertaken  to  determine  and  properly 
represent,  at  least  to  first  order,  these  wake  distortions.  It  should  be 
remembered  that  the  wake  distortions  are  the  result  of  a  time -varying  three- 
dimensional  induced  velocity  field  and  cannot  be  properly  determined  by  only 
the  induced  velocities  at  the  blades  or  in  the  plane  of  the  rotor.  However,  it 
is  noted  that  reasonably  good  agreement  with  measured  airloads  has  been 
achieved  using  a  regular  skewed  helical  wake  which  corresponds  to  a  constant 
and  uniform  induced  velocity  field.  Thus,  for  the  cases  computed,  this  is 
evidently  a  reasonable  first  approximation. 

The  measured  and  computed  airloads  are  presented  here  for  two  different 
rotors  and  advance  ratios.  They  are:  (1)  the  NASA  model,  which  was  tested 
in  a  wind  tunnel  at  =  0.  15  (Reference  4),  and  (2)  a  full-scale  HU-1A 
helicopter  rotor  tested  in  flight  at  =  0.  26  (Reference  5).  To  emphasise 
the  fact  that  apparent  good  agreement  with  respect  to  azimuthal  distributions 
does  not  necessarily  guarantee  good  agreement  with  respect  to  the  harmonics, 
the  following  two  forms  of  presentation  are  used:  (1)  the  azimuthal  distributions 
(i.  e. ,  time  histories)  at  each  radial  station,  and  (2)  the  radial  distribution  of 
the  sine  and  cosine  components  of  each  of  the  first  eleven  harmonics. 

For  the  NASA  model  at  j*-  *  0.  15,  the  computed  and  measured  azimuthal 
airload  distributions  (time  histories)  are  presented  for  fire  radial  stations  in 
Figure  3.  The  computed  distributions  have  been  adjusted  to  have  the  same 


mean  value  as  the  measured  distributions  at  each  of  the  radial  stations,  so 
that  the  nonuniform  parts  of  the  loads  can  be  directly  compared.  In  general, 
the  agreement  appears  to  be  good.  A  more  stringent  comparison  is  made  in 
Figure  4  where  the  radial  distributions  of  the  sine  and  cosine  component  of 
each  harmonic  of  the  airloads  are  presented.  It  is,  here,  obvious  that  when 
the  computed  and  measured  radial  distributions  of  the  individual  harmonics 
are  compared,  the  correlation  does  not  appear  to  be  as  good;  it  is,  however, 
very  encouraging. 

For  the  HU-1A  at  yti  =0.  26,  the  airloads  were  computed  with  and  without 
the  measured  motions  of  the  blades;  the  computation  without  the  measured 
motions  did,  however,  have  the  steady  blade  deflections  and  the  first  harmonic 
rigid  flapping  motion.  The  measured  motions  used  included  the  flapwise 
bending,  rigid  flapping,  rigid  pitching,  and  the  torsion  modes;  only  those 
harmonics  of  these  motions  which  exceeded  an  equivalent  recording  trace 
deflection  of  1/50  of  an  inch  were  retained.  Presented  in  Figure  5  are  the 
computed  and  measured  azimuthal  distributions  of  the  airload  at  each  radial 
station;  the  computed  results  in  this  figure  are  for  the  case  with  the  elastic 
motions  included  and  have  been  adjusted  to  have  the  same  mean  value  as  the 
measured  results  at  each  radial  station.  The  agreement  here  appears  to  be 
quite  good,  although  it  is  not  as  good  as  was  obtained  for  the  NASA  model  at 
J/~  -  0. 15.  Also  indicated  on  each  airload  time  history  in  this  figure  is  the 
azimuth  range  where  the  computed  instantaneous  effective  angle  of  attack 
exceeded  14°  (i.  e. ,  the  azimuth  range  for  which  the  lift  coefficient  was  limited 
to  a  stalled  value  in  the  computations).  It  is  interesting  to  note  the  change  in 
character  of  these  airloadings  on  the  retreating  side  of  the  disk  from  the  out¬ 
board  radial  stations  tc  the  inboard  radial  stations.  This  is  evident  in  both  the 
computed  and  measured  results  and  is  apparently  the  effect  of  stall  occurring 
on  the  inboard  blade  sections.  The  radial  distributions  of  the  sine  and  cosine 
components  of  each  harmonic  of  these  measured  and  computed  airloads  are 
presented  in  Figure  6.  Here,  the  results  of  both  computations  are  presented, 
i.  *.  .  and  without  the  measured  motions.  The  most  significant  differences 

4  w 

between  the  two  computed  results  appear  in  the  first  through  the  fifth  harmonics 
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because,  with  the  exception  of  the  fourth  harmonic  motions,  only  the  first 
five  harmonics  of  the  motions  were  of  enough  significance  to  be  retained; 
the  fourth  harmonic  motions  were  insignificant  apdj  thus,  neglected. 

It  is  interesting  to  note  that  even  though  there  were  no  fourth  harmonic  motions 
used  in  the  computation,  a  significant  difference  did  appear  between  the  fourth 
harmonic  airloads  obtained  from  the  two  computations  which  can  only  be 
attributed  to  the  other  harmonics  of  the  motions.  Here,  as  for  the  NASA  model 
results,  the  correlation,  based  on  the  individual  harmonics,  does  not  appear  to 
be  as  good  as  that  based  on  the  azimuthal  distributions  (i.  e. ,  time  histories); 
the  results  are,  however,  very  encouraging  for  this  initial  effort. 

Some  additional  interesting  results  from  the  HU-1A  computation  which 
used  the  significant  motions  are  presented  in  Figures  7  through  9.  In  Figure 
7,  the  induced  velocity  distribution  is  presented  as  a  contour  plot  over  the 
rotor  disk;  these  induced  velocities  at  the  /arious  points  of  the  disk  are  the 
velocity  on  the  blade  when  it  is  at  the  point,  i.  e. ,  they  are  the  induced  veloci¬ 
ties  in  the  rotating  coordinate  system.  It  should  be  noted  that  the  contours  are 
not  of  equal  increments  of  the  induced  velocity.  The  two  most  interesting 
features  evident  in  this  figure  are  the  nonuniformity  of  the  induced  velocity 
distribution,  especially  for  an  advance  ratio  of  0.  26,  and  the  large  induced 
velocities  and  steep  gradients  on  the  retreating  half  of  the  rotor  disk;  the 
latter  tend  to  alleviats  the  so-called  retreating  blade  "tip-stall"  by  reducing 
the  instantaneous  effective  angles  of  attack.  This  is  evident  in  Figure  8  which 
presents  a  contour  plot  of  the  nondimens ionali zed,  instantaneous,  effective 
angle  of  attack  over  the  rotor  disk;  the  angle  of  attack  has  been  nondimension- 
alized  to  an  angle  of  14°  (the  estimated  stall  angle  for  these  blade  sections). 
Experimental  wind  tunnel  results  obtained  by  the  NASA  (Reference  6)  seem,  in 
general,  to  confirm  these  results.  The  tuft  patterns  presented  for  operating 

conditions  which  are  well  into  the  stall  region  appear  to  indicate  that  the  out- 

— / 

board  portion  of  the  blades  remained  unstalled.  Further,  it  was  observed  in 
the  tunnel  tests  that  the  total  lift  did  not  tend  to  fall  off  as  the  apparent  stall 
area  increased,  and  the  variation  of  longitudinal  flapping  with  advance  ratio 
did  not  exhibit  the  change  of  slone  that  might  he  evper.teH  hari  there  Keen  pro¬ 
gressive  tip  stall.  The  large  induced  velocities  and  steep  gradients  on  the 


retreating  side  of  the  disk  and  their  absence  on  the  advancing  side  are  related 
to  the  concentration  of  the  trailing  vorticity  near  the  blade  tips  and  shifting 
of  the  airloading  toward  the  tips  on  the  retreating  side.  This  shift  of  the  air- 
loading,  which  is  influenced  by  the  inboard  area  of  stall,  is  evident  in  Figure 
9  (a  contour  plot  of  the  airloading  over  the  rotor  disk). 

The  induced  drag  and  its  harmonic  analysis  are  also  a  part  of  the  output 
of  this  program.  For  this  case,  the  first  through  third  harmonics  of  the 
induced  drag  are  sizeable  and  resonance  conditions  could  be  toublesome  at 
the  higher  harmonics.  This  induced  drag  distribution,  when  used  to  compute 
the  induced  power,  yields  a  result  which  is  approximately  three  times  that 

% 

obtained  by  a  conventional  method  based  on  the  uniform  inflow  assumption. 

T!ie  conclusion  to  be  drawn  from  these  results  is  that  more  rc^dstic  induced 
velocity  distributions  will,  in  general,  yield  a  higher  induced  power  expend!- 
ture.  This  is  consistent  with  fixed-wing  theory. 


*  These  results  are  presented  and  discussed  in  more  detail  in  Reference  3. 

**  A  k..  ULT  TJf  r. 1 _ tT>  _ _  ", \ 
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GENERAL  DISCUSSION  AND  RECENT  WORK 


In  the  initial  approach  which  has  been  diucusssd,  the  dynamic  response 
aspect  of  the  problem  was  bypassed  to  facilitate  concentration  of  effort  on 
the  development  of  the  essential  features  of  the  computational  procedure  for 
determining  the  airloading.  As  mentioned  previously,  this  procedure  was 
possible  since  computed  airloads  were  compared  with  loads  measured  in 
flight  for  configurations  where  the  dynamic  response  was  also  measured. 
However,  we  must  remember  that  our  objective  is  to  predict  the  aeroelastic 
response  and  airloads  on  a  blade  when  our  only  known  inputs  are  given  flight 
conditions  (i.  e. ,  cyclic  pitch,  collective  pitch,  advance  ratio,  etc. ). 

Two  methods  have  been  considered  fo  *  solving  for  the  aeroelastic 
response.  One  consists  of  an  iterative  procedure  where  the  airloads  would 
be  computed  for  given  inputs  and  dynamic  response,  and  these  airloads  are 
then  used  to  obtain  a  second  estimate  of  the  dynamic  responses.  The  second 
method  involves  the  simultaneous  solution  of  the  aerodynamic  load  and 
dynamic  response  problems. 

It  is  clear  that  the  iterative  procedure  has  some  fundamental  difficulties 
from  a  theoretical  standpoint,  although  it  may  be  adequate  in  some  cases.  The 
airloads  computed  at  any  step  of  the  Iteration  are  based  on  the  prescribed 
inputs  and  dynamic  responses  of  the  preceding  iteration,  and  they  include  the 
corresponding  aerodynamic  damping  loads.  There  appears  to  be  no  exact 
method  available  for  separating  the  aerodynamic  damping  from  this  total  air- 
loading,  although  an  approximate  estimate  of  the  damping  can  be  made  with 
quasi -static  aerodynamics. 

The  dynamic  response  for  the  next  step  of  the  iteration  can  be  computed 
treating  the  total  aerodynamic  force  (including  aerodynamic  damping)  as  the 
exciting  force.  In  this  case,  the  responses  near  resonances  are  very  large, 
being  only  limited  by  the  comparatively  small  structural  damping,  and  it  would 
be  expected  that  the  convergence  of  the  prousdoxe  wwuld  be  peer.  Ae  improved 
convergence  can  be  obtained  by,  first,  approximately  estimating  the  aerodynamic 
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damping  and  subtr-  g  it  from  the  total  computed  airloading  to  obtain  the 
forcing  function.  Now,  the  same  approximate  treatment  of  aerodynamic 
damping  would  be  included  in  the  computation  of  the  dynamic  responses  for 
the  next  step  of  the  iteration. 

Currently,  we  are  evaluating  the  advantages  of  the  exact  treatment  of 
the  problem  which  closes  the  loop  on  the  airload  and  response  calculations. 

It  turns  out  that  the  general  technique  which  has  been  described  for  handling 
the  airloads  problem  alone  can  be  extended  to  include  the  determination  of 
the  dynamic  responses. 

In  steady  flight,  the  response  in  each  of  the  dynamic  modes  of  interest 
can  be  expressed  as  a  Fourier  series  with  the  fundamental  period  of  one  rotor 
revolution.  Thus,  for  example,  the  amplitude  velocity,  and  acceleration  in 
the  q^1  dynamic  mode  can  be  written, 

nSlt  +  bh  xw  ff:lr>£lYj(-Q<>ai*nA.t+bnc+a,nnt) 

6  n  0  «  V 

^  -  H-fl.  Y  (  c*<tnjn-t  •+■  bn  Avru  nSlt') 

However,  this  harmonic  description  of  the  response  is  incompatible  with  the 
time  domain  description  which  is  used  in  the  airload  computation  discussed 
earlier.  The  situation  is  easily  remedied  by  expressing  the  Qn  and  bn 
coefficients  in  terms  of  the  N  values  of  (j  at  the  A /  azimuth  positions  used 
in  the  airload  determination.  By  this  means,  the  velocity,  ^  ,  and  acceler¬ 

ation,  j  ,  at  a  given  azimuth  position,  can  be  written  in  terms  of  the  J  's 
at  all  Al  azimuth  positions.  It  also  proves  possible  to  express  the  structural 
damping  in  each  mode  in  terms  of  the  N  f 

When  the  dynamic  response  is  considered  simultaneously  with  :he  airloads, 
an  equation  of  motion  is  added  for  each  of  the  5  modes  considered  which  must 
be  satisfied  at  each  of  the  N  azimuth  positions.  Thus,  there  are  W  times  3 
new  variables  to  be  included  in  the  machine  Cilcjijtion,  In  addition  to  the 
equations  of  motion  obtained,  Equation  (6)  is  also  modified  when  the  dynamic 
response  is  determined  simultaneously.  Some  of  the  prescribed  inputs,  , 
become  linear  functions  of  the  new  variables.  However,  the  revised  equations 
can  still  be  arranged  in  the  same  form  and  solved  in  an  analogous  manner. 
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It  has  been  found  conceptually  convenient  to  view  the  aerodynamic  loading 
aspect  of  the  over -all  aeroelastic  problem  in  the  following  manner:  The  wing 
and  its  "immediate"  attached  wake  can  be  thought  of  as  a  "wing-wake  system" 
for  generating  a  force.  This  system  will  have  a  characteristic  transfer  function 
which  will  describe  the  variation  of  the  output  (i.  e. ,  the  airloading)  with  the 
frequency  of  the  input  (which  is  a  time -varying  angle  of  attack).  The  remainder 
of  the  wake,  which  is  the  returning  wake  of  the  wing  and  all  the  wake  of  the 
other  wings  in  the  rotor,  and  the  motions  of  the  wing  are  then  thought  of  as 
providing  the  input  to  this  system.  For  example,  when  the  classical  two- 
dimensional  oscillating  wing  problem  is  viewed  in  this  way,  the  wing  section 
and  its  entire  wake  will  make  up  the  wing -wake  system,  and  the  inputs  are 
just  due  to  the  motions  of  the  wing.  When  the  immediate  wake  of  the  rotating 
wing  is  considered  as  part  of  the  system,  it  will  influence  the  system  transfer 
function  and,  therefore,  be  effective  in  controlling  the  system  response  to  the 
inputs  rather  than  functioning  as  part  of  the  input  source.  The  point  of  division 
between  the  immediate  wake  and  the  remainder  of  the  wake  is,  here,  unim¬ 
portant  because  the  model  is  used  only  as  a  convenient  way  of  thinking  about 
the  problem  and  not  as  a  basis  for  computation.  The  problem  is,  thus,  con¬ 
ceptually  divided  into  two  parts:  the  first  part  being  the  inputs  to  the  wing -wake 
system,  and  the  second  part  the  wing-wake  system  itself. 

The  displacement  time  history  of  the  wake  will  certainly  influence  the 
inputs  to  the  wing-wake  system,  and  it  has  been  determined  that  the  airloading 
can  be  sensitive  to  this  displacement  time  history.  At  any  instant  of  time,  the 
wake  will  appear  to  be  a  skewed  and  distorted  helical  surface.  The  skewed 
helical  surface  upon  which  the  distortions  can  be  thought  to  be  superimposed 
is  determined  by  the  rotor  translational  velocity,  the  tip -path -plane  angle,  and 
the  rotational  velocity  of  the  rotor.  The  deviation  of  a  particular  wake  element 
from  this  surface  at  any  instant  of  time  depends  m  the  induced  velocities  it 
has  encountered  since  it  left  the  wing.  Because  th-^  first  revolution  of  wake  i« 
the  strongest  source  of  input  to  the  wing-wake  system,  it  should  be  positioned 


with  the  greater  precision.  Although  the  computational  procedure  of  the 
initial  effort  is  capable  of  accounting  for  any  rather  general  wake  displacement 
time  history,  there  still  remains  the  question  of  the  correct  displacement  time 
history.  To  determine  this  exactly  would  require  that  the  total  induced  veloci¬ 
ties  be  computed  throughout  the  wake  for  each  instant  of  time  and  that  the  wake 
elements  be  moved  according  to  their  local  induced  velocities  for  each  incre¬ 
ment  of  time;  such  a  procedure  becomes  computationally  impractical,  even  on 
a  high-speed  digital  computer.  A  possible  approximation  which  is  being 
attempted  would  require  repeating  the  loads  computation  at  least  once  {i.  e. , 
at  least  one  iteration.)  to  determine  the  proper  wake  distortion.  First,  the 
wake  displacement  time  history  would  be  assumed  as  it  is  now  for  the  initial 
method.  The  velocity  time  histories  of  selected  wake  elements  would  then  be 
computed  as  well  as  the  velocities  at  the  blades.  These  wake  element  veloci¬ 
ties  would  then  be  used  to  generate  the  wake  displacement  time  history  for  the 
second  computation.  It  is  anticipated  that  the  position  deviation  from  the 
assumed  displacement  time  history  will  not  significantly  alter  the  induced 
velocity  at  each  wake  element;  therefore,  the  convergence  of  the  position 
should  be  rapid.  Additional  iterations  would  be  carried  out  to  test  this  postu¬ 
lated  behavior. 

In  addition  to  having  an  adequate  representation  of  the  input  source  for 
the  wing-wake  system,  a  successful  computational  procedure  must  have  an 
adequate  representation  of  the  wing -wake  system  itself.  In  the  initial  approach, 
the  wing  was  represented  by  a  single  bound  vortex  at  the  quarter -chord, 
satisfying  boundary  conditions  at  the  three -quarter -chord  and  the  immediate 
wake  was  represented  by  a  system  of  concentrated  vortex  filaments  shed  at 
equal  increments  of  time.  The  major  approximations  involved  in  this  initial 
approach  are:  (1)  the  boundary  conditions  are  satisfied  at  only  one  point  on 
the  chord,  and  (2)  the  immediate  wake  is  represented  by  concentrated  vortices 
equally  spaced  in  time  from  the  trailing  edge.  These  approximations  of  the 
wing -wake  system  seem  relatively  drastic,  but  they  produce  relatively  good 

nr^u ^  wAaiilfe  fVtuii  far  raXfja  liainff  tVn  ft  t* sinr ^ * ftntAHon 

may  be  due  to  the  relatively  low  reduced  frequency  range  of  the  predominent 
part  of  the  airloading  and/or  may  also  be  the  result  of  compensating  errors 
arising  from  the  approximations. 
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From  a  mathematical  point  of  view,  the  boundary  conait’  >n  at  any 
instant  of  time  over  the  chord,  to  be  satisfied  by  a  distribution  of  sing’ilarities 
on  the  chord,  is  the  sum  of  all  the  velocities  normal  to  the  chord  at  thai 
instant  of  time.  This  boundary  condition,  which  is  termed  the  "instantaneous 
effective"  boundary  condition,  is  composed  of  the  normal  velocities  due  to  the 
angle  of  attack,  camber,  gusts,  etc.  (i.  e. ,  due  to  all  the  inputs  to  the  wing- 
wake  system),  plus  the  normal  velocities  induced  by  the  immediate  wake;  thus, 
it  is  through  this  influence  on  the  instantaneous  effective  boundary  condition 
that  the  immediate  wake  is  effective  in  controlling  the  wing-wa^e  system 
response.  In  the  airfoil  problem,  this  instantaneous  effective  boundary 
condition  is,  of  course,  related  to  past  history  of  the  singularity  distribution, 
i,  u ,  the  airfoil  problem  is  a  succession  of  related  boundary  value  problems. 

If  the  singularity  distribution  (i.  e.  ,  the  distribution  of  bound  vorticity) 
is  expanded  in  a  Glauert  series  and  the  instantaneous  effective  boundary  con¬ 
dition  in  a  cosine  series,  the  corresponding  coefficients  of  the  two  series  can 
be  shown  to  be  directly  proportional  to  each  other.  Because  the  lift  and  moment 
arc  functions  of  only  the  first  four  terms  of  the  Glauert  expansion,  they  therefore 
dep  md  only  on  the  first  four  terms  of  the  instantaneous  effective  boundary 
condition  expansion.  Thus,  to  predict  the  lift  and  moment,  the  wake  repre¬ 
sentation  need  only  represent  adequately  the  first  four  terms  of  the  expansion 
of  the  wake -induced  velocity  distribution  (i.  e. ,  the  low  order  part  of  the 
distribution)  and  then  the  wing  representation  should  be  capable  of  satisfying 
this  boundary  condition  over  the  entire  chord. 

The  amount  of  improvement  in  the  computed  lift  which  would  result  from 
the  use  of  such  improvements  in  the  computational  model  is  not  clear.  However, 
it  is  believed  it  would  yield  much  more  realistic  aerodynamic  moments;  thus, 
such  improvements  are  being  implemented  to  determine  their  benefits. 

The  chordwise  boundary  condition  will  be  satisfied  over  the  entire  chord 
by  directly  computing,  for  each  element  of  wake  vorticity,  its  contribution  to 
each  of  the  first  four  Glauert  coefficients  (instead  of  computing  the  induced 
velocity  at  the  three -«»»*?*»••• -f'fcord  due  to  each  element  as  in  the  initial  approach). 
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This  is  accomplished  by  the  use  cf  Glauert’ s  well-known  integral  which  is 
transformed  to  apply  to  the  problem  coordinate  system.  This  amounts  to 
determining  the  first  four  terms  of  the  cosine  e.rpansion  of  the  chordwise 
induced  velocity  distribution  due  to  an  element  of  vorticity  arbitrarily  oriented 
with  respect  to  tne  blade  chord.  It  is  clear  that  the  more  distant  a  wake 
element  i&  from  the  chord,  the  less  will  be  its  contribution  to  the  higher -order 
terms  of  its  cosine  expansion.  Thus,  most  of  the  wake  (specifically  that  part 
functioning  as  the  input  source  to  the  wing-wake  system)  will  contribute 
significantly  only  to  the  first  two  Glauert  coefficients  (those  which  correspond 
tc  the  linear  and  constant  part  of  the  instantaneous  effective  boundary  condition 
over  the  chord).  The  immediate  wake  (that  part  which  functions  as  part  of  the 
wing -wake  system)  will,  of  course,  contribute  significantly  to  all  four  of  the 
Glauert  coefficients. 

There  are  two  possibilities  for  the  representation  of  the  immediate  wake. 
First,  there  is  evidence  £r<->m  some  of  our  earlier  computations  that  by  proper 
choice  of  the  distance  of  the  first  concentrated  shed  wake  element  from  the 
trailing  edge,  it  may  be  pot-ible  to  reproduce  the  first  four  terms  of  twe 
expansion  of  the  wake -induced  velocity  distribution  with  sufficient  accuracy 
to  enable  computation  of  the  lift  and  moment  adequately  for  the  reduced  frequency 
range  of  interest.  The  second  possibility  is  a  formulation  which  has  been 
developed  wherein  the  shed  wake  is  effectively  made  a  continuous  distribution 
behind  the  wing  from  the  trailing  edge  back  as  far  as  deemed  necessary.  This 
distribution  is  not  fixed;  for  each  azimuth  position  of  the  blade,  this  distribution 
will  be  a  function  of  the  total  circulation  about  each  blade  section  for  each  of  the 
azimuth  positions.  Both  of  these  wake  representations,  together  with  the  above 
method  of  satisfy  ng  the  chordwise  boundary  condition-,  are  being  tried  on  the 
two-dimensional  oscillating  wing  problem  in  order  to  determine  their  adequacy. 

In  the  initial  approach,  the  mesh  of  shed  and  trailing  vortex  filaments 
was  maintained  for  the  entire  retained  wake.  It  is  now  believed  that  it  may  be 
possible  to  eliminate  all  the  shed  vorticitv  for  most  flight  conditions,  except 
that  in  the  immediate  part  of  the  wake;  this  could  reduce  the  computation  time 
by  approximately  one -half. 


CONCLUDING  REMARKS 


It  should  be  remembered  that  the  over -all  problem  is  an  aeroelastic 
problem  and  the  airloading  is  but  one  aspect  of  it.  Because  of  the  inter¬ 
dependency  of  the  airloading  and  structural  response,  it  may  be  necessary 
to  compute  them  simultaneously  in  order  to  obtain  satisfactory  results. 

Direct  application  of  such  a  formulation  to  aeroelastic  stability  analyses 
is  evident,  Another  interesting  area  of  application  for  this  work  is  the 
maneuvering  or  transient  flight  conditions  of  the  rotor. 

The  work  accomplished  thus  far  indicates  that  a  practical  solution  to 
the  aeroelastic  problem  is  feasible.  However,  there  are  many  possible 
approximations  and  assumptions  which  may  be  utilized,  and  those  which  are 
permissible  are  not  yet  clearly  defined.  Thus,  it  may  still  be  possible  to 
improve  the  results  while  simplifying  or,  at  least,  shortening  the  computation, 
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NOMENCLATURE 


Q.  -  Lift  curve  slope: 

b  -  Blade  semichord  at  each  radial  station, 

0 

h  -  "Effective"  plunging  velocity  of  the  blade  section;  includes  the 
component  of  the  free  stream  velocity,  V4  ,  normal  to  the  local 
blade  sections  (i.  e.  ,  the  effects  of  the  local  blade  slopes). 

k  -  Subscript  indicating  position  in  the  rotor  disk. 

^  Lift  per  unit  span. 

h  -  Number  of  radial  bound  vortex  segments  representing  blades. 

N  -  Number  of  equally  spaced  azimuth  positions  at  which  the  airloads 
are  computed;  must  be  multiple  of  the  number  of  blades  in  the 
rotor. 

P  -  The  point  of  the  rotor  disk  where  the  airloads  are  computed. 

r  -  The  radial  coordinate  of  the  rotor. 

/?  -  Blade  radius. 

V  -  Tangential  velocity  in  tip -path -plane:  V*_fl  r  +  Vj 

\(f-  Free  stream  velocity  due  to  rotor  translation. 

(xr-  The  induced  velocity  component  normal  to  the  tip -path -plane  at 
the  three-quarter-chord  of  the  section. 

o<  -  Instantaneous  geometric  angle  of  attack  with  respect  to  the  tip- 
path -plane. 

A 

d.  -  Stalling  angle  of  attack  for  the  blade  sections. 

Instantaneous  effective^  angle  of  attack  with  respect  to  the  tip- 
path -plane:  <*«. »  c 4  ur/j 

<*r-  Longitudinal  tip-path -plane  inclination  with  respect  to  the  free 
stream  velocity,  Vfyt  . 

P  -  The  unknown  strength  of  the  bound  vortex  segment  representing 
the  blade 
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A  "  Advance  ratio:  ju.  •  ty  co-*  p 

/°  -  Air  density. 

^  -  Aaimuthal  position. 

JX  -  Rotational  speed  of  the  rotor 
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Figure  2  SAMPLE  MAP  OF  COMPUTATION  CONTROL  POINTS  IN  TIP  PATH  PLANE  FOR 
FOR  n  =  4,  N  =  16 
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ROTOR  BLADE  CHORIttlSE  PRESSURE  DISTRIBUTIONS 
AND  TWO-DIMENSIONAL  AIRFOIL  CHARACTERISTICS 

By  James  Scheiman  and  Henry  L.  Kelley 
NASA  Langley  Research  Center 

SUMMARY 

A  comparison  is  made  between  airfoil  chordwise  pressure  distribu¬ 
tions  from  helicopter  jtor  flight  tests  and  static  two-dimensional 
wind-tunnel  tests.  Differences  in  actual  and  two-dimensional -airfoil 
pressure  distributions  are  shown  to  exist.  These  differences  in  air¬ 
foil  characteristics  are  expected  to  amplify  the  blade  flap wise  and 
torsional  vibratory  forces  determined  from  two-dimenBional-airfoil 
data.  Possible  reasons  for  these  airfoil  differences  are  briefly  dis¬ 
cussed.  The  point  is  made  that  in  endeavoring  to  confirm  current 
refined  theories  of  calculating  section  angle  of  attack,  it  is  essential, 
in  making  data  comparisons,  that  care  be  used  to  prevent  these  differences 
between  — uul  and  static  two-dimensional- section  data  from  obscuring 
the  effectiveness  of  the  angle-of-attack  calculations. 

INTRODUCTION 

Experience  has  shown  that  the  ability  to  perform  an  adequate 
structural  dynamic  analysis  of  the  rotor  blade  is  marginal.  This  lack 
of  ability  has  generally  been  viewed  as  attributable  to  unknown  air 

1  norla  anr!  ^  n  r\nwM  aiO  nv>lrv</-v**v\  -t  "VC"1  “  “  rC»v**C  I*  wavU*  v/C 

the  applicability  of  two-dimsnsional-airfoil  characteristics.  This 
view  has  tended  to  be  confirmed,  for  example,  by  the  results  of  rotor 
test  inflow  velocity  measurements  and  by  the  adequacy  of  predicting 


helicopter  performance  by  use  of  two-dimensional  airfoil  data.  In  any 
case,  both  inflow  velocities  and  actual  airfoil  characteristics  must 
be  known  in  order  to  perform  a  reasonably  accurate  dynamic  blade 
analysis. 

In  regard  to  the  inflow  velocities,  it  is  believed  that  the 
capability  of  theory  to  predict  these  velocities  for  triLi  level  flight 
has  significantly  improved  recently.  With  these  new  theories  the 
danger  exists  for  blaming  the  remaining  inadequacy  of  the  inflow  theory 
for  any  lack  of  correlation  between  test  and  theory,  when  the  differences 
may  actually  be  caused  by  airfoil  characteristic  discrepancies. 

The  validity  of  two-dimensional  data  has  been  given  little  detailed 
attention  because  of  a  lack  of  actual  operating  test  data.  Partly  to 
help  fill  this  gap,  the  NASA  Langley  Research  Center  has  recently 
completed  a  helicopter  flight- test  program  which  has  utilized  extensive 
blade  pressure  instrumentation.  These  data  provide  a  comp ail  son  of 
the  actual  and  two-dimensional-airfoil  chord-wise  pressure  distributions 
to  the  extent  needed  to  illustrate  that  important  airfoil  characteristic 
discrepancies  do  occur  in  the  flight  conditions  sampled. 

Portions  of  the  flight  measured  chordwise  pressure  distributions 
for  two  flight  conditions  are  discussed.  Samples  of  these  distribu¬ 
tions  are  directly  compared  with  two-dimensional  full-scale  data  (see 
ref.  l)  by  equating  the  two  normal  force  coefficients.  The  chordwise 
pressure  distribution  for  other  flight  conditions  and  the  movement  of 
the  blade  center  of  pressure  are  discussed. 
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All  reference  to  two- dimensional -airfoil  characteristics  in  this 
paper  refers  to  static  two-cL'msnsional  characteristics  in  distinction 
to  oscillating  unsteady  two-dimensional  characteristics. 

SYMBOLS 


c 

CK 

Ap 

r 

E 

VF 

x 

x 

h 

t 


airfoil  chord 
normal-force  coefficient 

differential  pressure  measured  on  the  airfoil 
dynamic  pressure 

radial  distance  to  blade  element  measured  from  center  of 
rotation 

blade  radius  measured  from  center  of  rotation 
forward  speed 

chordwise  distance  measured  from  blade  leading  edge 
center  of  pressure  of  airfoil  section  measured  from  leading 


edge 

nondimensional  tip-speed  ratio, 


blade  nominal  azimuth  angle,  measured  from  downwind  position 


in  the  direction  of  rotation  and  disregarding  blade  lag 


motion 


n 


rotor  angular  velocity 
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DISCUSSION 

Normal-Force  Coefficient  for  Flight 
With  Blade  Section  Stall 

A  plot  of  the  local  normal-force  coefficients  C^  along  the  blade 
radius  for  different  azimuth  positions  is  shown  in  figure  1.  The 
flight  condition  is  for  a  trim,  level-flight-cruise  forward  speed  and 
a  reduced  rotor  rotational  speed  and  is  thus  for  a  flight  condition 
ejected  to  produce  local  blade-section  stalling.  Further  details  of 
this  flight  condition  are  available  in  table  IV  of  reference  2. 

Notice  the  high  values  of  normal-force  coefficient  in  the  area  of 
r/R  *  0.55,  *  ■  210°  to  2^0°  and  for  r/R  =  0.75,  *  »  250°. 

These  coefficients  correspond  to  dynamic  pressures  of  approximately 
50  pounds  per  square  foot  at  r/R  =0.55  and  100  pounds  per  square 
foot  at  r/R  =*  0.75 •  The  normal-force  coefficient  values,  in  these 
areas  of  the  rotor,  correspond  to  values  above  the  maximum  static 
two-dimensional  values  of  Just  prior  to  these  high  normal- 

force  coefficients  a  rapid  rate  of  change  in  the  normal-force 
coefficient  is  noted.  This  change  in  Cjj  can  be  directly  related 
to  a  two-dimensional-airfoil  angle-of-attack  change  and  the  corresponding 
high  rates  of  angle-of-attack  change  can  be  explained,  for  example,  by 
the  rapid  changes  in  local  Inflow  velocities  through  the  rotor.  In 
this  instance  an  estimate  based  on  successive  normal-force  coefficients, 
in  the  previously  mentioned  high  angle-of-attack  area  of  the  rotor, 
indicates  a  rate  of  roughly  100°  per  second  or  1°  per  2~l/2  blade- 
chord  lengths.  This  rapid  angle-of-attack  increase  will 
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provide  a  partial  explanation  for  the  lack  of  ehordwise-pressure- 
dlstrlbution  correlation  explored  further  in  this  paper. 

The  circles  on  this  figure  are  points  where  chordwise  pressure 
distributions  are  discussed  in  figures  2  to  *)•>  the  solid  circles  are 
points  vhere  two-dimensional  and  flight  data  are  compared. 

Chordwise  Pressure  Distributions  for  Plight 
With  Blade  Section  Stall 

Figures  2  to  4  are  for  the  same  flight  condition  as  figure  1, 
which  was  selected  with  the  expectation  of  producing  local  blade- 
section  stalling.  A  plot  of  the  chordwise  pressure  coefficient  distri¬ 
bution  for  r/R  =  0.55  arid  *  =  165°,  195°,  210°,  225°,  and  255°  is 
shown  in  figure  2.  The  blade-azimuth  position,  Integrated  normal-force 
coefficient,  and  the  centers  of  pressure  are  as  indicated.  At  i  *  165°, 
the  flight-test  distribution  agrees  with  the  two-dimensional  dataj  the 
center  of  pressure  is  close  to  (slightly  aft  of)  the  quarter  chord. 

The  normal-force  coefficient  is  below  the  two-dimensional  stall  point. 

At  i  =*  195°  the  normal-force  coefficient  of  1.3  is  above  the  two- 
dimensional  stall  value  but  the  pressure  distribution  appears  unstalled. 
At  the  remaining  azimuth  locations  the  normal-force  coefficient  is  above 
the  airfoil  section  two-dimensional  stall  point  and  no  two-dimensional 
data  are  available  for  comp  aid  son.  The  pressure  distribution  is  such  as 
to  correspond  to  some  separation  and,  therefore,  the  section  can  be 
viewed  as  exhibiting  stall  characteristics  although  the  details  of  the 
distribution  have  no  counterpart  in  two-dimensional  data.  Based  on 
examination  of  the  contours  of  figure  1,  this  increased  maximum  Cjj 
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is  expected  to  increase  the  vibratory  amplitude  of  the  actual  airfoil 
blade  loads  as  compared  with  the  predicted  two-dimensional  air  loads. 

It  may  be  of  interest  that  some  variations  in  pressure  occurred 
between  rotor  revolutions  for  these  last  three  plots  (average  values 
are  shown) j  however,  the  distribution  shapes  are  believed  representative. 
These  variations  in  themselves  suggest  a  stalled  type  of  flow. 

Sample  pressure  distributions  for  r/R  =  0.75  are  shown  in 
figure  3;  these  distributions  are  for  the  same  flight  condition  shown 
in  figures  1  and  2.  For  t  =  165°  195°  the  distribution  shows 

good  agreement  with  two-dimensional  data.  As  the  azimuth  angle  increases 
from  t  =*  195°  to  2^0°,  the  normal-force  coefficients  again  increase  to 
values  above  the  two-dimensional  stall  point,  although  the  actual  air¬ 
foil  retains  the  unstalled  two-dimensional  pressure  distribution. 

The  pressure  distribution  for  r/R  =  0.95  is  shown  in  figure 
The  normal-force  coefficients  are  all  below  the  static  two-dimensional 
stall  point  and  therefore  good  pressure-distribution  correlation  would 
be  expected.  For  f  =  ^5°  and  75°  the  agreement  between  the  flight  and 
two-dimensional  data  is  indeed  reasonable,  but  at  ^  »  90°  and  120°  the 
correlation  is  not  so  good.  Thus,  while  a  large  part  of  the  rotor  does 
behave  in  accordance  with  two-dimensional  data,  figures  2  to  4  show 
that  poor  correlation  can  occur  to  a  degree  which  would  be  expected  to 
have  a  major  effect  on  periodic  blade  loads. 
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Chordwise  Pres sure  Distributions  for  Flight 
With  High  Blade -Tip  Mach  Numbers 

The  next  trim  level  flight  condition  discussed  for  flight  at  a 
high  tip  Mach  number,  for  which  the  maximum  blade- tip  Mach  number 
was  0.76.  The  chordwise  pressure  distribution  for  the  O.95  blade 
station  is  shown  in  figure  5.  The  normal- force  coefficients  are  all 
below  the  two-dimensional- airfoil  stall  point  and  good  correlation 
would  be  expected.  The  centers  of  pressure,  however,  are  all  farther 
forward  than  would  be  expected,  even  for  high  Mach  number  operation  of 
a  two-dimensional  airfoil.  For  t  =  30°  there  is  reasonable  agreement 
between  the  flight  and  two-dimensional  distributions,  although  for 
V  31  75°,  90°,  and  105°  the  flight  data  depart  from  the  two-dimensional 
data. 

The  0.75-radius  station  shown  in  figure  6  is  for  the  previously 
described  high  tip  Mach  number  flight.  The  flight-measured  chordwise 
pressure  distributions,  the  centers  of  press-are,  and  normal-force 
coefficients  are  typical  of  unstalled  two-dimensional  data.  -The 
correlation  shown  is  good. 

The  high  tip  Mach  number  test  for  the  0.55-blade-span  station  is 
shown  in  figure  7.  Again  the  normal-force  coefficients,  centers  of 
pressure,  and  the  distribution  are  typical  of  two-dimensional  data. 

The  correlation  with  two-dimensional  data  is  again  good. 

In  {summary,  figures  5,  6,  and  7  for  the  high  tip  Mach  number  flight 
indicate  that  a  large  percentage  of  the  actual  chordwise  pressure  distri¬ 
butions  are  in  agreement  with  two-dimensional  airfoil  data.  Only  a 
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small  (though  important)  percentage  of  the  pressure  distributions  are 
not  in  agreement.  Since  the  disagreement  in  this  case  is  primarily 
in  the  high  Mach  number  regions,  it  appears  that  with  careful  selection 
of  the  flight  test  conditions,  it  will  be  possible  to  find  cases  that 
warrant  comparison  with  theories  using  the  two-dimensional  data  to 
study  the  adequacy  of  new  angle-of-attack  prediction  theories. 

Other  Flight  Conditions 

A  small  portion  of  the  chordwise  pressure  distributions  for  a 
number  of  other  trim  level-flight  conditions  have  been  reviewed  and 
the  results  were  similar  to  the  results  of  the  previous  two  flight 
conditions  discussed  in  detail  in  this  paper;  namely,  that  portions  of 
the  actual  operating  helicopter  blade  do  not  behave  In  accordance  with 
two-dimensional  airfoil  data.  Because  there  are  these  cases  where 
important  differences  do  arise,  an  exact  knowledge  of  the  rotor  inflow 
velocities  is  not  necessarily  sufficient  to  describe  the  exact  rotor 
blade  loading.  Caution  should  therefore  be  exercised  in  interpreting 
the  correlation  of  flight  measured  and  theoretical  rotor-blade  span- 
wise  loadings. 


Measured  Center-of- Pressure  Movement 
In  an  attempt  to  generalize  the  actual  airfoil  center-of-pressure 
movement,  a  plot  was  made  of  the  center  of  pressure  as  a  function  of 
the  tmiuiutu  eoigle  iur  nnree  airrerent  flight  conditions,  and 

this  plot  is  shown  in  figuie  8.  Note  tne  forward  shift  in  center  of 
pressure  on  the  advancing  aide  of  the  rotor  for  all  three  flight 
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conditions  and  the  rearward  center  of  pressure  on  the  retreating  side 
of  the  rotor  for  the  first  flight  condition  (flight  with  blade  section 
stall  discussed  previously) .  The  forward  shift  in  center  of  pressure 
could  not  be  explained  by  Mach  number  effects. 

The  flight- test  blade  was  a  modified  KACA  0012  airfoil  (ref.  1) 
and  hence  had  no  camber.  The  interesting  possibility  thus  arises  that 
if  a  small  amount  of  camber,  which  would  tend  to  add  a  constant  moment 
coefficient  with  varying  angle  of  attack  below  stall  were  added,  the 
variations  in  dynamic  pressure  with  azimuth  would  then  modify-  the 
measured  centers  of  pressure  in  such  a  way  as  to  result  in  reduced 
one-per-re volution  aerodynamic  control  forces.  In  other  words,  the 
added  source  of  moment  variation  with  azimuth  would  be  expected  to  have 
a  phase  angle  such  as  to  offset  partially  the  measured  variations. 

Discussion  of  Actual  and  Two-Dimensional -Airfoil 
Pressure-Distribution  Differences 

The  reasons  for  the  differences  found  between  actual  and  two- 
dimensional  airfoil  data  are  not  conpletely  understood.  As  is  well 
known,  the  flow  conditions  on  a  rotor  are  highly  complex  and  many 
potential  contributing  explanations  have  long  been  at  hand  should 
such  problems  arise .  Since  the  problem  has  now  been  verified  m 
tangible  fora,  an  effort  is  being  made  to  sort  out  some  of  these 
possibilities. 

As  one  example,  the  fact  that  a  high  rate  of  increase  in  angle 
of  attack  can  give  higher  than  static  C^y  values  is  well  known 
(for  example,  ref.  3),  and  this  effect  has  long  been  looked  for  in 
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rotor  measurements .  In  some  early  investigations  this  effect  vas 
apparently  negligible,  that  is,  stall  vas  evidenced  roughly  where 
expected.  Apparently  the  other  complexities  of  rotor  inflow,  and  the 
specific  design  details,  prevented  any  significant  occurrence  of  higher 
than  static  values.  In  several  more  recent  investigations, 

including  the  present  one,  the  opposite  has  been  true,  Dynamically, 
this  effect  would  be  expected  to  increase  the  actual  amplitude  of  the 
oscillating  air  loads  as  compared  to  the  calculated  loads  based  on 
two-dimensional  data. 

It  should  be  noted  that  the  most  drastic  source  for  high  rates 
of  change  of  angle  of  attack  is  likely  to  be  the  striking  of  the  tip 
vortex  from  the  previous  blade.  Consequently,  the  high  rates  of  change 
and  the  Cjj  values  in  excess  of  static  two-dimensional  values  may 
occur  in  specific  cases  in  basically  mild  flight  conditions  as  well  as 
in  the  low  rotor  speed  or  high  forward  velocity  conditions  normalfy 
associated  with  blade- section  stalling. 

Time-varying  blade  yaw  angles,  spanwise  flow  on  the  blade,  and 
nonuniform  velocity  gradients  in  front  of  the  airfoil  are  other 
possible  factors  that  may  cause  disagreement  between  actual  and  two- 
dimensional  airfoil  characteristics. 

CONCLUDING  REMARKS 

It  has  been  shown  that  the  actual  helicopter  rotor  blade  does  not 
always  behave  in  accordance  with  two-dimensional  airfoil  data.  These 
airfoil- characteristic  differences  are  expected  to  amplify  both  the 
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*pvi  -  and  torsional  blade  oscillating  loads.  Possible  reasons  for 
-•-  JL  f  fere  nee  e  were  briefly  explored.  The  point  is  made  that  when 
*  ads  are  predicted  from  refined  inflow  theories  as  compared  with 
« . loadings,  caution  should  be  exercised  in  interpreting 
-t's  in  blade  loading,  since  these  may  arise  because  of  the 
up l 1  c ability  of  two-dimensional  data  rather  than  inadequacies 
-  r.flow  theory.  Thus,  before  comparisons  of  actual  and  predicted 
s»io  sire  uccd  to  dstsnni.no  validity  of  GJigls—of ~s.tts.clc  celLcuIsl** 
h  each  experimental  case  used  must  be  reviewed  for  evidencr.  of 
.-eset.ee  or  absence  oi  discrepancies  between  the  actual  section 
mzic  characteristics  as  reflected  by  chor dulse  pressure  dietri- 
.  a:.t  che  section  characteristics  being  assumed  in  the  analysis. 
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Figure  2.-  Chordwise  pressure  distributions. 
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Figure  Chordwise  pressure  distributions. 
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figure  7-~  Chordwlse  pressure  distributions. 
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A  DISCUSSION  OF  ROTOR  BLADE  HARMONIC 
AIRLOADING 


R.  H.  Miller 

Department  of  Aeronautics  and  Astronautics 
Massachusetts  Institute  of  Technology 


Introduction 


The  subject  of  this  symposium  covers  one  ot  the  most  important  and  yet  in  many  ways 
one  of  the  most  neglected  areas  of  V/STOL  technology.  Certainly  in  an  engineering, 
sense  the  subject  cannot  be  said  to  have  been  neglected  since  the  problems  of  helicopter 
vibration  and  fatigue  (both  human  and  structural)  must  be  reduced  to  manageable  propor¬ 
tions  by  the  exercise  of  ingenuity  and  much  time  and  expense  during  the  engineering 
development  stages,  and  this  has  been  the  history  of  V/STOl  development  since  the  first, 
practical  helicopter  flew  in  the  1930s. 

However  the  basic  research  which  could  lead  to  a  better  understanding  of  the  physical 
phenomena  of  VTOL  dynamics  and  hence  to  a  more  rational  solution  of  the  associated 
problems  has  been  neglected  in  the  past  and  the  result  has  been  a  containment  of  the 
problems  rather  than  their  solution.  There  are  many  valid  reasons  for  this  lag  in  research, 
the  most  obvious  being  the  extreme  complexity  of  the  system.  However  we  now  have  both 
experimental  and  analytical  techniques  available  to  us  which  could  permit  rapid  investi¬ 
gation  of  the  aerodynamics  and  dynamics  of  V/STOL  aircraft  and  which  could  lead  to  a 
better  understanding  of  the  basic  phenomena.  It  is  axiomatic  that  without  this  understanding 
control  of  these  phenomena  will  never  be  possible  except  on  a  hit  and  miss  basis. 

The  decision  to  hold  a  symposium  on  this  topic  is  therefore  most  welcome  and  timely.  In 
responding  to  the  invitation  to  present  a  discussion  at  this  meeting  on  work  conducted  at  MIT 
on  rbtor  dynamics,  a  selection  was  made  from  the  projects  conducted  during  the  past  few  years 
and  only  those  pertinent  to  the  title  of  this  paper  and  of  current  interest  were  selected.  Re¬ 
search  by  others  in  the  same  area  will  not  be  referred  to  in  the  interests  of  time  and  because 
a  fairly  complete  coverage  has  previously  been  included  in  Ref.  1. 


Early  Investigation  of  Airloads  by  Meyer  and  Falabella 

From  1947  to  1952  a  series  of  wind  tunnel  tests  were  conducted  on  the  aerodynamic  loading 
of  rigid  and  flapping  rotor  blades.  The  initial  studie:  (Ref.  2)  were  confined  to  the  measure¬ 
ments  of  bending  moments.  One  of  the  most  interesting  results  was  the  demonstration  of  the 
critical  effects  of  blade  resonance  combined  with  the  very  low  aerodynamic  damping,  clearly 
evident  from  Fig.  1.  Reasonable  overall  agreement  between  computed  and  measured  bending 
moments  was  obtained  for  both  the  flapping  and  rigid  rotor  but  it  was  concluded  that  information 
regarding  the  actual  aerodynamic  loading  or  induced  flow  was  needed  to  establish  the  higher 
harmonic  loadings. 
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Consequently  a  new  program  was  started  in  which  the  blade  airloads  were  measured 
directly  by  means  of  a  pressure  transducer.  Fig.  2,  from  Ref.  3  presents  one  of  the  results 
of  this  investigation  and  shows  the  downwash  distribution  at  the  rotor  disc  deduced  from 
the  measured  airloads.  These  results  are  useful  for  verification  of  more  recent  theoretical 
calculations  of  the  downwash  distribution  and  at  the  time  provided  experimental  evidence 
of  the  highly  non  uniform  character  of  the  rotor  inflow  which  is  now  recognized  as  being 
one  of  the  primary  sources  of  the  periodic  excitation  of  helicopters. 

It  was  during  this  program  that  the  pronounced  effect  of  flapping  hinge  offset  on  the 
airload  distribution  of  the  rotor  was  noticed.  This  led  to  the  conclusion  that  the  effect 
of  a  rigid  roior  cculd  be  achieved  without  the  need  for  providing  a  high  degree  of  struc” 
tural  stiffness;  that  is  a  hinged  rotor  could  be  designed  to  carry  large  first  harmonic  bending 
moments  at  the  hub  by  offsetting  the  flapping  hinge  appreciably  from  the  center  of  rotation. 

Of  interest  here  is  the  pronounced  influence  of  this  first  harmonic  load  variation  on  the  in 
flow  through  the  rotor,  as  shown  in  Fig.  3.  A  much  more  uniform  downwash  distribution 
was  obtained  which  would  suggest  that  an  equivalent  rigid  propeller  carrying  large  first  hafmontc 
moments  in  forward  flight  would  be  smoother  than  a  flapping  rotor.  Structural  and  design 
considerations  determine  the  degree  to  which  such  moments  could  be  tolerated  in  an  actual 
vehicle.  This  point  is  discussed  further  in  Ref.  4. 


Direct  Measurement  of  Hub  Loading 


The  investigations  described  above  clearly  indicated  that  appreciable  higher  harmonic 
excitations  were  occurring  on  the  rotor.  Simple  rotor  blade  aerodynamic  theories  using 
uniform  downwash  indicated  that  these  loads  would  decrease  rapidly  with  advance  iatio. 
However  this  clearly  did  not  fit  the  observed  facts.  Consequently  a  program  in  which  hub 
loads  were  directly  measured  was  initiated  in  order  to  explore  this  phenomenon  more  fully. 
The  program  was  conducted  in  two  parts.  In  the  first  the  rotor  was  subjectea  to  forbad  os 
ci Nations  and  the  results  substantiated  the  large  reduction  in  damping  predicted  by  Loewy  ^ 
due  to  unsteady  aerodynamic  effects.  In  the  second  part  the  rotor  was  used  as  a  mechanical 
amplifier  for  the  aerodynamic  harmonic  load  by  placing  the  blade  in  resonance  with  the 
harmonic  to  be  investigated  end  measuring  the  resultant  hub  shears.  As  a  result  of  this 
series  of  tests  it  was  obvious  that  a  mechanism  existed  which  produced  high  periodic  airloads 
on  a  rotor  blade  at  integers  of  the  rotor  speed,  that  existing  methods  for  predicting  these 
loads  were  inadequate  and  that  unsteady  aerodynamic  effects  in  the  presente  of  these  oscil¬ 
latory  airloads  were  of  prime  :mportance,  (Refs.  1 1  and  12) 

The  experimental  program  was  continued  in  order  to  investigate  more  harmonics,  parti¬ 
cularly  the  third  which  had  not  been  covered  in  the  earlier  investigation  and  in  order  to 
determine  the  effects  of  tandem  rotor  geometry  on  the  vibration  level.  Because  of  the  ob 
vious  difficulty  of  changing  geometry  for  such  purposes  in  a  flight  test  program  results  of  this 
nature  are  most  conveniently  obtained  from  wind  tunnel  tests.  This  program, conducted  by 
Duvivier,  produced  the  results  shown  in  Fig.  4,  taken  from  Ref,  5  in  which  the  pronounced 
effects  of  geometry  on  tandem  rotor  vibration  are  evident. 


During  this  program  techniques  were  developed  for  measuring  blade  airloads  directly  by 
means  of  single  pressure  pickups  located  at  various  spanwise  stations  along  the  rotor  blade. 
This  permits  the  investigation  of  many  interesting  phenomena  with  relatively  small  models 
in  the  wind  tunnel  and  forms  the  basis  of  a  continuing  program  involving  the  direct  measure¬ 
ment  of  rotor  blade  airloads  under  various  simulated  flight  conditions  and  airframe  geometries. 
Because  of  the  difficulties  of  covering  a  wide  variation  of  trim  conditions  and  geometry  on  a 
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Aeroelastic  Investigations 

Simultaneously  research  was  being  conducted  on  the  flutter  characteristics  of  rotor 
blades.  Much  of  this  work  is  not  pertinent  to  this  meeting,  however  two  items  are  of 
interest  and  will  be  briefly  discussed. 

Stall  flutter  of  propellers,  rotors  and  compressors  has  been  a  problem  for  vears  because 
of  the  high  angle  of  attack  at  which  these  operate.  On  a  rotor  blade  this  phenomenon  is 
believed  to  be  one  of  the  causes  for  the  very  high  pitching  moments  experienced  at  high 
advance  ratios.  Investigations  by  Ham  (Ref.  6)  defined  the  parameters  which  govern  this 
phenomenon  and  provided  a  valuable  criterion  by  which  its  presence  can  be  detected.  For 
example  Fig.  5  taken  from  Ref.  6  shows  the  great  dependence  of  stall  flutter  amplitude  on 
the  amount  of  stalled  area  in  the  lotor  disc  regardless  of  advance  ratio. 

Another  aeroelastic  phenomenon  which  is  of  importance  in  establishing  the  harmonic 
loading  of  rotor  blades  is  that  arising  from  the  deflections  of  the  blade  out  of  'he  plane 
of  rotation  and  out  of  the  plane  of  its  twist  axis.  The  effects  of  this  phenomenon  on 
flutter  are  discussed  in  Ref.  7.  Its  potential  significance  for  blade  loading  is  discussed 
in  Refs,  1  and  8  and  may  be  demonstrated  by  a  simple  illustration.  In  Fig.  6  a  blade  is 
shown  with  its  normal  elastic  deflection.  Primary  torsional  flexibility  for  simplicity  is 
assumed  to  occur  at  the  feathering  hinge  located  near  the  blade  root.  As  the  blade  twists 
about  the  feathering  hinge  components  of  centrifugal  force  act  to  increase  this  twist.  That 
is,  a  rotor  blade  bent  out  of  its  plane  of  rotation  is  not  in  static  equilibrium.  Since  the 
centrifugal  force  is  several  times  the  gross  weight  of  the  ship,  clearly  powerful  blade 
twisting  moments  can  be  produced  by  this  mechanism.  In  Ref.  7  it  was  shown  that  in  cer¬ 
tain  cases  these  moments  couid  be  equivalent  to  those  produced  by  a  6°/o  shift  in  chord- 
wise  CG  of  the  rotor  blades. 

Since  rotor  blades  must  bend  elastically  in  the  presence  of  aerodynamic  lift,  the  cen¬ 
trifugal  force  component  as  well  as  the  stead  state  drag  will  cause  periodic  twisting  moments 
due  to  the  periodic  changes  in  blade  bending  deflection  caused  by  the  higher  harmonic  air¬ 
loads.  In  addition,  periodic  changes  in  ind"ced  drag  due  to  the  nigher  harmonic  content 
of  the  downwash  will  combine  with  the  steady  state  elastic  deflection  of  the  blade  to  produce 
another  source  of  periodic  twisting  moments.  Twisting  moments  from  these  two  sources  will 
cause  harmonic  changes  in  angle  of  attack  and  are  therefore  of  considerable  importance  in 
establishing  the  blade  loads. 


Theoretical  Prediction  of  Blade  Airloads 


The  tests  discussed  above  clearly  indicated  the  need  for  an  analytical  tool  for  computing 
blade  downwash  velocities  which  would  take  into  account  the  individual  blade  wake  geometry 
and  also  introduce  the  effects  of  unsteady  aerodynamics.  Attempts  to  obtain  a  closed  form 
solution  to  this  problem,  or  one  based  on  tabulated  integrals,  were  not  successful  and  it  was 
evident  that  extensive  computer  facilities  would  be  required  to  explore  this  problem  and, 
hopefully,  to  provide  a  basis  for  obtaining  simplified  solutions  suitable  for  engineering  ap¬ 
plications.  In  1960  the  availability  of  an  IBM  709  computer  at  the  MIT  Computation  Center 
and  funds  from  a  Carnegie  grant  permitted  initiation  of  such  a  program.  This  is  a  continuing 
investigation  under  BuWeaps  sponsorship.  The  results  to  date  have  been  published  in  Ref.  i, 

8  and  9  and  therefore  only  the  intent  of  the  research  and  the  current  status  will  be  summarized 
here. 
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If  was  decided  that  an  attempt  would  be  made  to  construct  the  program  of  research  in 
such  a  way  as  to  lead  primarily  to  an  understanding  of  the  basic  phenomena  involved  and 
to  include  insofar  as  possible  a  rigorous  trearment  of  the  unsteady  aerodyndmics.in  order  to 
accomplish  both  these  purposes  the  following  approach,  was  followed. 

1.  The  computer  program  was  developed  in  such  a  way  that  the  rotor  inflows  could 
be  obtained  independently,  if  possible,  of  blade  motions. 

2.  The  classical  approach  in  which  the  trailing  and  shed  wakes  were  treated  inde- 
pendently  and  the  unsteady  aerodynamic  effects  determined  separately  was  followed. 

3.  The  program  could  be  modified  to  a  lifting  surface  rather  than  a  lifting  line  ap~ 
proach  when  necessary. 

4.  Since  the  program  was  to  be  used  primarily  for  research,  numerical  precision  of 
the  different  integrations  would  be  varied  over  wide  limits  and  over  the  different  ele¬ 
ments  of  the  program  as  required. 

5.  The  classical  assumptions  of  a  rigid  wake  of  propeller  vortex  theory  would  not 
be  accepted. 

It  was  hoped  that  this  approach  in  addition  to  clarifying  the  physics  of  the  problem  would 
also  lead  to  a  much  simplified  computer  program. 

Computer  solution  times  vary  widely  depending  on  the  degree  of  accuracy  desired.  Fig. 

7  shows  a  comparison  of  the  correlation  between  theory  and  experiment  as  a  function  of  the 
number  of  stations  used  in  establishing  the  trailing  vortex  pattern  for  a  harmonic  of  general 
interest.  The  large  number  of  stations  of  Fig.  7a  are  required  to  define  the  tip  loading  ac¬ 
curately  but  the  smaller  number  of  stations  of  Fig.  7b  would  probably  be  sufficiently  accurate 
for  most  engineering  purposes.  The  two  station  solution  (Fig.  7c)  is  useful  for  demonstrating 
many  of  the  phenomena  of  rotor  blade  loading  and  is  a  valuable  research  tool.  At  moderate 
speeds  the  bound  circulation  on  the  blade  does  not  vary  radically  over  the  outer  5C'/o  span, 
regardless  of  twist,  from  the  minimum  energy  condition  of  constant  circulation  and  hence  much 
qualitative  information  may  be  quickly  obtained  by  observing  behavior  of  o  hypothetical 
"ideally1’  twisted  blade,  as  was  done  in  Ref.  i,  having  a  tip  vortex  and  a  roct  vortex  located 
at  a  representative  inboard  station. 

The  degree  of  accuracy  shown  in  Fig.  7a  is  by  no  means  typical  f jr  all  rotors  and  har¬ 
monics.  For  certain  conditions  more  accurate  solutions  are  required  using  much  closer  in¬ 
tervals,  a  reiteration  for  non-rigid  wake  effects  and  a  lifting  surface  theory  rather  than  the 
lifting  line  approximations.  Such  solutions  and  the  conditions  under  which  they  are  used 
will  be  discussed  later.  The  ideal  fluid  potential  flow  model  used  in  these  solutions  has, 
however,  certain  fundamental  limitation!  and  there  is  some  question  as  to  the  degree  of  ul¬ 
timate  computation  accuracy  which  would  be  meaningful  in  view  of  both  these  limitations, 
and  the  difficulty  of  obtaining  equally  accurate  test  results  or  actual  rotors  under  known  con¬ 
ditions.  It  is  believed  that,  as  in  the  case  of  fixed  winqs,  tb  s  theory  is  of  most  value  as  a 
research  tool  in  order  to  clarify  physical  phenomena  a;'J  to  o’  in  the  interpretation  of  test 
results. 

The  decision  to  separate  the  blade  motions  from  the  computer  solutions  tor  downwash 
resulted  from  the  previously  discussed  experimental  evidence  of  the  pronounced  effect  of 
blade  elasticity  on  the  airloads.  if  it  could  be  assumed  that  the  blade  responded  to  ars. 
aerodynamic  input  which,  to  first  order,  would  not  be  modified  by  this  response  then  the 
problem  could  be  considerably  simplified  since  the  solution  for  blade  airload  and  for  blade 
response  to  that  airload  could  be  separated.  The  basis  for  this  hypothesis  came  from  the 
investigations  reported  in  Ref.  1  From  this  study  it  was  evident  that  interharmonic  coupling 
effects  could  to  first  order  be  neglected.  That  is,  for  example,  the  nth  component  of 

downwash  generated  by  an  mth  harmonic  lift  variation  was  an  order  of  magnitude  less  then 
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that  generated  by  the  nth  harmonic  lift  variation.  Furthermore  higher  harmonic  airloads 
are  an  order  of  magnitude  less  than  the  steady  state  load  acting  or.  the  ship.  Consequently 
their  contribution  to  the  downwash  distribution  at  the  rotor  disc  will  be  1  ’*o  orders  of  mag" 
nitude  less  than  that  contributed  by  the  spiral  wake  generated  by  steady  state  blade  Jift  dnd 
this  downwash  is  readily  computed  for  any  gross  weight  and  mean  inflow. 

Having  established  the  nature  of  the  downwash  variation  around  the  rotor  disc  the  next 
step  is  to  compute  the  airload  variation  corresponding  to  this  downwash  distribution.  Quasi¬ 
static  lifting  line  solutions  give  reasonable  estimates  of  these  loads  but  introduce  appreciable 
errors  both  in  magnitude  and  phase.  It  is  therefore  necessary  to  introduce  the  classical  lift 
deficiency  functions  of  unsteady  aerodynamics.  In  Ref.  1  a  closed  form  solution  was  ob¬ 
tained  for  this  lift  deficiency  function,  C(k)  =  F  +  iG  (to  use  standard  nomenclature)  for 
hovering  flight  and  for  an  infinite  number  of  blades  and  was  shown  that  F  and  G  reduced 
to  the  following  simple  expressions,  independent  of  frequency. 


F  = 


G  =  0 


cri t 

JX 


B 


Si 
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This  result  leads  to  the  hypothesis  that  over  all  approximate  values  of  the  lift  deficiency 
function  F  and  the  corresponding  phase  shift  G/F  could  also  be  obtained  for  the  forward 
flight  cose.  Consequently  typical  values  were  computed  some  of  which  ore  shown  in  Fig. 

8.  A  more  marked  variation  with  frequency  was  apparent  than  For  the  hovering  flight 
case;also  the  lift  deficiency  function  gradually  increases  with  p  and  eventually  approaches 
the  two  dimensional  classical  value.  Similarly  G  increases  from  a  value  very  close  to  zero 
ar  d  again  approaches  the  two  dimensional  va lue away  from  the  blade  tip.  G,  which  deter" 
mines  the  phase  shift  cannot  be  predicted  by  using  lifting  line  theory  because  of  the  singularities 
involved  as  the  shed  wake  approaches  the  trailing  edge  of  the  surface,  and  lifting  surface 
theory  must  rherefore  be  used  in  obtaining  these  results.  An  approximate  means  of  de¬ 
termining  I'his  while  still  retaining  the  convenience  of  lifting  line  theory  in  the  computer 
solutions  was  sujM^sted  in  Ref.  8  in  which  an  equivalent  upper  limit  for  the  shed  wake  in¬ 
tegrations  was  established.  Lift  deficiency  functions  obtained  by  these  means  are  highly 
accurate  in  comparison  with  the  forma!  approach  using  lifting  surface  theory  given  in  Ref. 


Computation  oi  the  rotor  blade  airload  thus  proceeds  in  the  folk  wing  steps: 

Compute  the  downwash  distribution  at  the  rotor  disc  due  to  a 
trailing  vortex  system  of  constant  unit  strength  originating  from  as 
many  spanwise  stations  as  desired  along  the  blade.  The  use  of  unit 
vertex  strength  permits  tabulating  this  downwash  distribution  for 
various  values  of  mean  inflow  ,  advance  ratio,p,  and  the  number 
of  blades.  This  tabulation  is  currently  being  started. 

2,  The  actual  downwash  is  then  computed  from  the  known  blade 
geometry  and  aircraft  trim  by  inversion  of  the  matrix  obtained  from 
Eq,  17  of  Ref.  8  in  the  manner  described  in  that  reference. 

3.  Knowing  this  downwash  distribution  quasi  static  airloads  may  be 
obtained  from  Eq.  14  of  Ref,  8  and  the  actual  airloads  obtained  by 
suitably  modifying  the  quasi  static  airloads  by  the  appropriate  F  and 
G  function*,  similar  to  those  aiven  in  Fig.  a. 
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Sjeps  2  -nd  3  are  trivial  solutions  for  the  computer,  the  bulk  of  time  being  required  for 
the  integration  of  step  I,  hence  the  atte.tr  ,o  tabulate  these  results. 

>he  limitations  of  the  lifting  line  theory  will  now  be  discussed  using  for  illustration 
the  simple  two  dimensional  case.  Fig.  9  v,  a  sketch  of  an  airfoil  in  two  dimensions  in 
the  presence  of  on  element  of  vorticity  /T  J '3  at  a  distance  £  from  the  airfoil.  The 
norma!  velocity  w  indue-. j  at  this  airfoil  at  a  chordwise  station  x  due  to  Y  d  j 


w(x)  - 


r d<; 

“*) 


When  f  >>  b,  i.  e. ,  when  the  element  of  vorticity  is  relatively  far  from  the  blade 
then  f  »  x  and  w  is  substantially  constant  with  x.  A  mean  value  of  w  may  then  be 
used  and  this  is  the  basis  for  the  well  known  lifting  line  theory.  When  £  is  not  large 
compared  to  h  lifting  surface  theory  must  be  used.  This  requires  that  w  be  expressed 
as  a  Fourier  es  in  terms  of  the  blade  chord.  The  first  two  components  of  the  series 
and  theii  tii  /ariation,>  define  the  blade  lift.  These  steps  have  been  described  in  de~ 

*[!  *  "  ^ '  Evidently  computer  time  for  such  a  solution  would  be  several  times  more 
than  that  required  for  the  iifting  line  solution,  depending  on  the  number  of  chordwise 
stanons  required  to  define  w(x)  adequately.  Fortunately  there  are  only  a  few  instances 
where  such  an  approach  would  be  required  and,  as  is  evident  from  the  above  discussion, 
the  multiple  chord  station  computations  need  be  oerformed  only  over  a  short  segment  of 
the  wake  in  the  vicinity  of  the  airfoil,  using  eqccMons  4  end  5  of  Ref.  1.  Two  cases 
requiring  this  more  rigorous  lifting  surface  theory  will  now  be  givan  as  illustration. 

The  first  is  in  the  computation  of  the  lift  deficiency  function  C(k)  =  F  +  iG  referred 
to  above  As  the  shed  wake  approaches  the  airfoil  or'g  -*  b,  a  singularity  at  the  trailing 
L*r?  •  develops  which  can  not  be  ignored  if  G,  which  determines  the  phase 

snitt,  is  to  be  correctly  evaluated.  This  requires  the  use  of  lifting  surface  theo.  v.  For 
reasons  discussed  in  some  detail  in  Ref.  8,  F  can  be  closely  approximated  for  the  low  re¬ 
duced  frequencies  of  interest  in  rotor  dynamics  by  ignoring  this  singularity  and  using  mean 
v°  u®s  of"W«!n  lifting  line  theory.  However  G  is  always  critically  dependent  on  the 
chordwise  variation  of  w(x).  This  was  the  reason  for  introducing  the  concept  of  a  "near" 
and  a  far^  wake  in  Ref.  1  and  using  lifting  surface  theory  in  the  near  wake.  In  hovering 
tlignt,  p  -  0,  the  near  and  far  wakes  contribute  about  equally  to  the  final  value  of  G  but 
as  p  increases  the  effect  of  the  far  wake  gradually  diminishes  until  at  a  p  of  about  ,  3  only 

the  effect  of  the  near  wake  on  G  is  significant.  Again  this  follows  from  the  above  discus¬ 
sion. 

In  Fig.  10  the  build  up  of  the  lift  deficiency  function  is  shown  fora  typical  case.  Sub¬ 
script  I  refers  to  the  lift  deficiency  function  generated  by  the  near  shed  wake  where  the 
near  wake  is  defined  as  being  that  portion  extending  from  the  rotor  blade  aft  over  the  first 
quarter  azimuth.  This  portion  of  the  lift  deficiency  function  generally  is  close  to  the  two 
dimensional  ciassical  va lue.  Subscript  2  refers  to  the  lift  deficiency  function  obtained  from 
both  the  near  and  far  shed  wakes.  Subscript  3  refers  to  the  lift  deficiency  function  as  obtained 
trom  both  the  far  shed  and  near  shed  wakes  and  the  cc-nplete  trailing  wa(ce  system.  It  is  evi¬ 
dent  that  as  the  advance  ratio  increases  the  near  shed  wake  dominates.  At  the  blade  tip, 
however  and  for  loadings  concentrated  over  a  short  segment  of  the  span,  the  three  dimensional 
value  of  become?  appreciably  smaller  than  the  two  dimensional  value.  Frequently  Its 
e  ec*  at  least  for  blade  loading  computations,  will  not  be  of  great  significance  and  lifting 
line  solutions  may  therefore  preside  good  approximations  to  this  loading. 

Another  case  of  interest  in  which  iifting  surface  theory  must  be  used  is  found  when  the 
,rr - ...  T<?wrs<D.  nils  ou ii  uttui  wnen  non t i g ta  wakes  are  snrro- 


duced  into  the  computations  at  low  advance  ratios  as  discussed  in  Ref.  10.  It  may  also 
occur  at  any  advance  ratio  in  the  case  of  a  rotor  with  a  large  number  of  blades  operating 
at  low  thrust  coefficients  and  hence  having  low  mean  inflow  velocities.  The  intense  tip 
vortex,  characteristic  of  rotors,  generated  by  any  one  blade  will  come  close  to  the  fol¬ 
lowing  blade  and  lifting  line  theory  then  predicts  erratic  and  unreliable  lift  variations. 

An  interesting  conclusion  drawn  from  an  examination  of  the  time  history  of  airloads 
both  experimental  and  analytical  is  the  very  pronounced  higher  harmonic  content  of  this 
load  close  to  the  90  azimuth  value.  Figure  16  of  Ref.  1  showed  the  very  rapid  change  in 
downwash  occurring  near  this  point  as  the  blade  passes  over  the  leading  edge  of  the  vortex 
spiral  and  the  resulting  lift  is  almost  in  the  nature  of  an  impulse.  It  was  found  that  even 
nine  harmonics  were  insufficient  to  describe  this  rapid  change  of  lift  and  good  agreement 
v/ith  test  was  only  obtained  in  this  regime  when  the  computer  solution  as  obtained  prior  to 
harmonic  analysis  was  used.  The  significance  of  this  result  lies  in  the  fact  ihat  a  vibratory 
input  from  the  rotor  of  frequency  equal  to  the  number  of  blades  times  the  rotor  speed  will 
always  occur  regardless  of  the  number  of  b laces  used.  Since  this  type  of  vibratory  input 
arising  from  an  impulse  is  particularly  unpleasant  any  rotor  modification  which  reduces 
the  impulsive  nature  of  the  lift  would  be  desirable. 

ir  may  well  be  asked  at  this  point  what  modifications  are  indicated  which  would 
improve  the  vibration  characteristics  of  the  rotor.  Certainly  twist  appears  tp  have  a  pro¬ 
found  effect  on  the  airload  distribution  with  he  non-twisted  blade  having  a  higher  downwash 
at  the  tip  than  the  twisted  blade.  This  would  tend  to  move  the  wake  further  away  from  the 
rotor  as  each  subsequent  blade  passes  over  it.  Another  factor  is  the  effect  of  blade  modifi¬ 
cations  on  the  higher  harmonic  spanwise  load  distribution.  The  response  of  a  rotor  blade  to 
a  given  input  depends  on  the  generalized  force  and,  if  the  spanwise  airload  distribution  is 
known,  this  response  con  be  modified  by  suitably  varying  the  spanwise  mass  distribution. 

The  value  of  a  better  understanding  of  rotor  aerodynamics  in  forward  flight  will  become 
more  evident  once  a  working  theory  has  been  developed  and  extensively  used. 


Summary 

It  has  been  shown  how  experimental  investigations  over  a  period  of  years  have  led 
to  the  development  of  analytical  tools  for  the  determination  or  'otor  blade  dynamic  loading 
and  how  these  tests  hove  helped  in  defining  objectives  for  the  analytical  solutions.  This 
work  represents  a  small  part  of  a  larger  effort  over  a  period  of  years  by  many  organizations. 
However  this  total  effort  is  still  small  compared  to  the  magnitude  of  the  problem  and  it  is 
hoped  that  the  discussions  in  this  paper  and  the  results  reported  by  other  investigators  will 
stimulate  greater  effort  in  this  important  area. 

In  particular,  refinements  in  the  theory  are  needed  to  permit  rigorous  treatment  of 
special  oases.  The  effects  of  blade  twist  on  the  harmonic  load  need  clarification.  Simpli¬ 
fied  engineering  solutions  must  be  developed,  possibly  the  effect  of  a  viscous  vortex  core 
introduced  and  the  effects  of  rotor/rotor  interference  determined  analytically.  Stall  flutter 
and  other  potential  instabilities  at  high  advance  ratios  due  to  the  elastic  deflection  of  the 
blade  out  of  the  plane  oi  rotation  must  be  furrher  investigated  preparatory  to  the  proposed 
increases  in  helicopter  speeds  during  the  next  decade. 

Such  problems  are  among  the  most  interesting  and  challenging  ones  yet  presented 
to  the  aeroelastician.  Progress  to  date  indicates  that  rational  solutions  to  mese  problems 
will  eventually  be  found  and  in  a  time  frame  depending  only  on  the  level  of  effort  expended. 
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FIG, 4  REAR  ROTOR  THIRD  HARMONIC 
VERTICAL  SHEAR 


FIG.  5  EFFECT  OF  DEGREE  OF  STALL  ON  MAXIMUM  FLUTTER  AMPLITUDE 


AXIS  OF  ROTATION 


FIG.  6a  BLADE  DEFLECTED  OUT  OF  PLANE 
OF  ROTATION 


FIG.  6b  PLAN  VIEW  OF  BLADE  DEFLECTED  OUT  OF 
PLANE  OF  ROTATION  AND  TWISTED  ABOUT 
FEATHERING  AXIS  SHOWING  COMPONENT 
OF  CENTRIFUGAL  FORCE  NORMAL  TO 
FEATHERING  AXIS  AND  CAUSING  FURTHER 
TWISTING 
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FIQ.7  CORRELATION  BETWEEN  THEORY  AND  TEST 
AS  A  FUNCTION  OF  THE  NUMBER  OF  8RANWJSE 
STATIONS  USED  IN  SOLUTION.  THIRD  HARMONIC 
SPANWISE  DISTRIBUTION  2  BLADES  p  *  .2 
3  SPIRALS  IS*  INTERVALS 

— — -  EXPERIMENTAL  RESULTS  FROM  NACA 
RM  i  R6!07 

-  THEORY 

(a)'  10  STATIONS* SOLUTION  TIME  17 MINS 
(i>)  6  STATIONS- SOLUTION  TIME  10 MINS 
io)  2  STATIONS -SOLUTION  TIME  4 MINS 
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FIG.  8 

THREE  DIMENSIONAL  LIFT  DEFICIENCY  F,  AND 
PHASE  SHIFT,  G/F  AS  A  FUNCTION  OF  HARMONIC 
n  OF  ROTATIONAL  SPEED.  THREE  BLADES, ASPECT 
RATIO  OF  BLADES  10.  INFLOW  RATIO  X  »  0.05 


FIG.  9  TWO  DIMENSIONAL  INDUCED  FLOW  ON 

AIRFOIL  DUE  TO  ELEMENT  OF  VORTICiTY 
IN  WAKE 
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1.  Introduction 

Recent  theoretical  work  at  M.  I.  T.  on  rotor  blade  harmonic  air  loading.  Refs.  1  and  2, 
has  indicated  the  possibility  cf  a  portion  of  a  blade  tip  vortex  being  drawn  up  into  the 
leading  edge  of  the  rotor  disc  under  certain  flight  conditions  with  consequent  major  effects 
on  the  blade  local  air  loading.  It  is  believed  that  these  effects  contribute  significantly  to 
the  vibration  ievel  of  the  rotor  during  the  associated  flight  conditions.  However  no  experi¬ 
mental  confirmation  of  the  existence  of  the  phenomenon  has  been  available. 

is  the  purpose  of  the  present  paper  to  first  discuss  the  theoretical  results  indicating 
the  possible  existence  of  this  phenomenon,  then  indicate  how  these  results  can  be  utilized 
in  establishing  those  test  conditions  that  will  permit  observation  of  the  phenomenon,  and 
finally  to  present  the  experimental  evidence  to  date  of  the  occurrence  of  the  phenomenon 
in  terms  or  blade  measured  air  loading. 


2.  Theoretical  Background 

Figures  1  and  2  illustrate  the  path  of  the  blade  tip  vortices  for  two  selected  advance 
ratios  corresponding  to  the  so-called  transition  regime  of  a  typical  helicopter.  Note  the 
several  intersections  of  the  forward  blade  by  the  tip  vortices  of  all  three  blades.  In  most 
flight  conditions  the  rotor  downwash  is  sufficiently  large  to  carry  these  tip  vortices  well 
below  the  rotor  disc.  However  under  certain  transition  conditions  a  pronounced  upwash 
occurs  at  the  leading  edge  of  the  disc  which  may  carry  the  intersecting  tip  vortices  up  through 
the  disc  and  back  down  again  farther  downstream  as  suggested  in  Ref.  3.  A  typical  example 
is  shown  in  Figure  3. 

At  any  point  (x,y)  on  the  trajectory  of  that  segment  of  the  tip  vortex  at  =  180°, 


Then  assuming  that  near  the  disc  v(x)  does  not  vary  significantly  with  y. 
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and  the  condition  for  vortex  intersection  of  the  disc  is  given  by  the  integral  equation 


A(*)  dx  =  0 

0 


where  x.  is  the  value  of  the  x-coordinate  of  the  trajectory  when  intersection  occurs, 
measured  along  the  disc  fore^and-aft  axis  from  the  leading  edge.  The  trajectory  is  readily 
obtained  graphically  as  shown  in  Fig.  3.  Therefore  theoretical  plots  of  rotor  downwash  at 
l/j  «  180can  be  used  to  predict  possible  intersections  of  the  rotor  disc  by  a  tip  vortex  at 
y/  =  180°.  Note  that  disc  incidence  .  is  of  prime  significance  in  determining  the  inter¬ 
section  point. 

Downwash  distributions  for  advance  ratios  of  .  05  and  .  10,  for  a  three  bladed  rotor  with 
twisted  blades  are  shown  in  Figures  4  and  5.  These  distributions  were  calculated  using  the 
theory  of  Refs.  1  and  2. 

Two  types  of  rotor  wake  geometry  were  considered  in  the  inflow  analysis.  The  first  type,(1), 
that  of  the  so-called  "norfr-rigid"  wake,  is  based  on  the  assumption  that  the  tip  vortex  moves 
with  the  local  resultant  velocity  if  experiences  at  the  instant  of  shedding.  The  second  type, (2),  also 
that  of  the  "non-rigid"  wake  is  rosed  on  the  assumption  that  the  tip  vortex  moves  with  the  mean 
local  resultant  velocity  it  experiences overan  interval.  Inflows  resulting  from  both  assumptions 
are  shown  in  the  figures.  Note  the  effect  of  a  non-rigid  wake  on  the  distributions.  Under  the 
given  conditions  the  non-rigid  wake  geometry  leads  to  a  large  predicted  upwash  near  the  lead¬ 
ing  edge  of  the  rotor  disc. 

Application  of  the  above  inflow  distributions  to  calculation  of  the  trajectories  of  tip  vortices 
at  ijt  =  180°  indicates  that  vortex  intersections  should  occur  near  the  80yo  radius  of  the  blade, 
as  shown  in  Figures  6  and  7. 

3.  Experimental  Technique 

The  experimental  equipment  used  in  the  present  study  was  that  described  in  Refs.  4  and  5, 
modified  in  the  present  case  to  permit  remote  control  of  the  angle  of  shaft  tilt  while  the  rotor 
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was  operating.  A  genera!  view  is  shown  in  Fig.  8.  Available  blades  had  pressure  trans¬ 
ducers  mounted  at  49^0,  83^0,  and  93°/o  of  biade  radius.  The  experimental  technique 
involved  the  recording  of  a  timing  signal  and  the  pressure  signal  from  the  pressure  trans" 
ducer  at  the  83^o  station  on  an  Offner  pen  recorder  for  a  rest  condition  having  parameters 
corresponding  to  those  for  a  theoretically  predicted  upwash  at  the  leading  edge  of  the  disc, 
while  varying  shaft  angle  in  the  range  indicated  by  the  theory  of  Section  2  to  lead  to 
vortex  intersections  of  the  forward  blade  near  the  83°/o  pressure  station.  When  the  pres¬ 
sure  trace  showed  peaks  indicating  the  occurrence  of  the  predicted  intersection  of  the  tip 
vortex  with  the  blade  near  the  83/o  radius  point  (see  Fig.  8),  a  complete  record  was 
taken  on  a  Heiland  oscillograph.  This  record  shov/ed  the  signals  from  all  three  pressure 
transducers  as  wall  as  a  timing  signal,  and  a  blade  flapping  signal.  Typ;aal  Offner 
pressure  record',  are  shown  in  Figs.  9  and  10. 


4.  Results  and  Discussion 


The  experimental  results  of  Figure  9  indicate  the  occurrence  of  a  large  local  pressure 
fluctuation  at  the  83°/o  blade  pressure  station  for  rotors  having  three  twisted  or  untwisted 
blades  in  simulated  transition  flight  at  p  =  .  05.  The  pressure  peak  occurs  at  an  azimuth 
position  of  approximately  ^0°,  Since  the  upwash  conditions  theoretically  predicted  for 
the  given  flight  condition  and  Iwistea  blades  indicated  a  tip  vortex  intersection  with  the 
rotor  disc  in  the  neighborhood  of  ip  —  180°,  it  is  concluded  that  the  observed  pressure 
fluctuation  is  due  to  the  tip  vortex  shed  by  the  second  preceding  blade  passing  near  the 
83°/o  pressure  station  of  the  instrumented  blade.  Note  that  the  effect  has  approximately 
the  same  intensity  for  both  twisted  and  untwisted  blades,  though  the  shaft  angles  differ 
for  the  two  cases. 

The  results  presented  in  Figure  10  for  p  =  .  10  are  less  clear,  but  are  believed  to  show 
the  effects  of  the  presence  of  a  tip  vortex  near  the  03^/a  pressure  station.  It  is  suspected 
that  in  this  case  the  tip  vortex  from  the  immediately  preceding  blade  passes  close  to  the 
83^0  pressure  station  over  a  much  larger  portion  of  the  disc  than  in  the  u  =  .05  case. 

Hence  the  effect  appears  as  a  small  dip  at  ip  =  130°  approximately  and  as  a  general  in¬ 
crease  in  pressure  loading  in  the  <^=150°  to  220°  region.  Again  the  results  are  approxi¬ 
mately  the  same  for  both  the  twisted  and  untwisted  blade. 

The  following  explanation  is  advanced  for  the  distributed  nature  of  the  effect.  Exami¬ 
nation  of  the  calculated  downwash  field  for  the  rotor  with  twisted  blades  in  the  p  =  .  10 
flight  condition  revealed  the  presence  of  an  upwash  occurring  at  the  i|/  =  127.5  azimuth 
position,  as  shown  in  Figure  11.  It  is  suspected  that  a  region  of  upwash  exists  between 
this  azimuth  position  and  that  at  (p  =  180°  and  between  blade  stations  at  80°/o  and  95^0 
radius.  Unfortunately  the  chosen  computation  points  for  the  theoretical  downwash  field 
calculation  did  not  include  this  region.  However,  existing  information  permits  the  con¬ 
struction  of  Figure  12.  The  two  plotted  points  represent  possible  disc  intersections  by  the 
tip  vortex  shed  from  the  immediately  preceding  blade.  It  is  assumed,  but  not  shown,  that 
a  locus  c'  intersections  exists  between  these  two  points.  The  sol'd  curve  represents  the 
instantaneous  position  of  the  tip  vortex  shed  by  the  preceding  blade  when  the  instrumented 
blade  is  at  the  corresponding  azimuth  position.  Finally  the  path  of  the  83^0  blade  pressure 
station  is  shown  by  the  dash-dotted  curve.  It  is  seen  that  the  tip  vortex  is  close  to,  and 
ahead  of,  the  instrumented  pressure  station  in  the  region  ip  =  130°  to  150°,  and  close  to, 
and  behind,  the  pressure  station  in  the  region  ip  =  150  to  180°.  The  observed  harmonic 
pressure  fluctuations  in  the  same  azimuth  regions  and  at  the  given  pressure  station  are  of  a 
form  consistenr  with  the  above  hypothesis,  i.  e. ,  first  a  down  loading  due  to  the  downwash 
of  an  arriving  vortex  and  then  an  up  loading  due  to  the  upwash  of  a  departing  vortex,,  The 
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above  discussion  is  speculative  only  and  awaits  quantitative  comparison  between  theoreti- 
cal  and  experimental  harmonic  airloading  for  the  givon  flight  conditions.  It  is  hoped  ro 
obtain  sucn  comparisons  in  the  near  future. 


5.  Conclusions 

The  existence  of  rotor  blade  Impulsive  airloading  due  to  the  impingement  of  tip  vor¬ 
tices  from  preceding  blades  is  now  believed  to  be  confirmed.  The  practical  significance 
of  the  tffect  includes  the  increase  of  rotor  operating  vibration  levels,  particularly  if  the 
Impulse  is  applied  to  the  anti-node  of  a  blade  bending  mode.  Flight  conditions  where 
the  effect  is  experienced  by  helicopters  include  transition  flight,  and  hovering  in  ambient 
winds,  for  example  while  performing  plane  guard  or  sonar  dipping  duties. 

Fouible  means  of  eliminating  such  impulsive  loading  include  rapid  acceleration  through 
transition,  since  large  disc  Incidences  ana  rotor  inflows  prevent  the  occurrence  of  the  up- 
wash  field  at  the  leading  edge  of  the  disc.  Another  approach  is  distribution  of  the  tip  vortex 
strength  by  appropriate  tip  geometry  or  by  tip  blowing,  as  suggested  in  Ref.  3. 

Finally,  the  tip  vortex  interference  effects  demonstrated  by  the  model  rotor  of  the 
present  paper  are  likely  to  influence  the  airloads  currently  being  measured  on  full  scale 
rotors  in  corresponding  flight  conditions  and  should  be  considered  when  evaluating  such 
airload  measurements. 
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Model  Rotor  Characteristics 


Blade  Radius 
Number  of  Blades 
Blade  Chord 
Airfoi1  Section 
Solidity 


R  ~  4  ft. 
3 

5  in. 

NACA  0012 
0.1 


Symbols 

Standard  NASA  nomenclature  Is  used  throughout  except  for: 
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FIO.  I  TIP  VORTEX  PATTERN  FOR  THREE 
BLADES'  n  *.06 
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TIP  VORTEX  TRAJECTORY 
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(b)  UNTWISTED  BLADE  0,*OO  a=-7° 


FIG.  9  EXPERIMENTAL  AIRLOAD  VARIATION  WITH 
AZIMUTH  AT  THE  83%  BLADE  STATION 
A  *.°5  867r*  10° 
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FIG.  10  EXPERIMENTAL  AIRLOAD  VARIATION  WITH 
AZIMUTH  AT  THE  83%  BLADE  STATION 
fi=.  10  S™*  10°  a  *  -  7° 
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FIG.  II  INFLOW  DISTRIBUTION  AT  *  » 127.8* 
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— *—  INSTANTANEOUS  TIP  VORTEX  POSITION 

- - 83%  BLADE  STATION 

o  PONT  ON  LOCUS  OF  INTERSECTIONS 
OF  DISK  BY  TIP  VORTEX 
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